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SYNOPSIS 


M.M. CHOUDHARY 
Ph.D. 

Department of Electrical Engineering 
Indian institute of Technology, Kanpur 

June, 1986 

MULTI-PHASE (n-phase) SYSTEMS : DYNAMICAL 
MODELLING AND ANALYSIS 

The local area demand of electric energy was 
previously met by a locally constructed generating station 
known as urban, system. However., gradually the power 
stations with higher capacity were installed which have 
been supplying the electric energy covering a wider area. 
Steadily with the installation of bigger plants, such as 
superhydel and superthermal plants near the source of 
energy for economic considerations, the electric energy is 
now being transmitted to the load centres far off from the 
generating stations. Thus the problem of electric utilities 
has been to transmit large amount of power economically and 
efficiently. Although the solution to this problem has 
been achieved to some extent by increasing the transmission 
voltage which ensures greater economy, reduced losses and 
improved regulation, yet the three-phase transmission 
system at the increased transmission voltage, for exanple, 
at EHV/UHV levels, presents certain drawbacks, such as diffi- 



culty in acquiring new rights-of-way, high surface voltage 
gradient on conductors, strong electric field on the ground 
surface with possible biological effects, increased audible 
noise and radio interference, switching surge overvoltages 
which cause problem to air gap insulation, problem in insu- 
lation co-ordination, increased short-circuit currents and 
possibility of ferrovresonance conditions. Because of these 
factors, increasing the. transmission voltage seems to have 
also reached a point of saturation. 

Further, HVDC transmission emerged as an alterna- 
tive which has some advantages, such as increased stability 
limit, lower losses and costs than the equivalent a.c. lines; 
it has the advantage of its use as an asynchronous connection 
which does not raise the fault level appreciably, and it can 
also be used as a link between a.c. systems having different 
frequencies. But, it has some diadvantages too, namely 

A 

higher terminal costs, introduction of harmonics resulting 
in greater losses and communication interference, and also 
difficulties in providing adequate protection scheme. 

The other alternative may be multi-phase i.e. high- 
phase order transmission (HPOT) system. The concept of the 
multi-phase system was initiated during the year 1972 and 
various studies have been conducted and reported since then* 



The multi-phase transmission system has the advantage of 
acquiring new rights-of-way efficiently and effectively. 

The general feasibility of multi-phase lines has been inves- 
tigated and the six-phase line appears to be more suitable 
than the other multi-phase lines. The feasibility studies 
of six-phase line obtained by converting double-circuit 
three-phase lines are also reported to have been carried out. 
The installation and the field testing of some experimental 
six-phase and twelve-phase lines are reported to have been 
undertaken. The overvoltages and insulation requirements 
for six-phase lines have also been studied. 

The various elements of six-phase systems for the 
purpose of steady-state analysis have been mathematically 
represented. The model of six-phase synchronous machines in 
phaser co-ordinates is reported for the analysis under balanced 
as well as under unbalanced condition. 

The mathematical model of three-phase/six-phase 
transformers with different primary and secondary connections 
with and without leakage impedances are available in the 
literature/ and based upon these models, the three-phase and 
the single-phase equivalents of a six-phase transmission line 
integrated with two three-phase/six-phase transformers have 
been developed. 



The load flow study is reported to have been 


carried out on, mixed 3-phase and six-phase, and completely 
six-phase systems under balanced as well as under unbalanced 
conditions. The impact of converting double -circuit 3-phase 
lines to single-circuit six-phase lines at the same line to 
line (i.e. adjacent phase) voltage on the load flow perfor- 
mance has been investigated. It has been concluded from 
these studies that, the mixed three-phase and six-phase 
system, and a completely six-phase system have the advantage 
of better voltage magnitude, improved phase angle and increa- 
sed efficiency, even for the higher system loadings than that 
of the conventional three-phase systems. 

The short circuit analysis, an important study for 

the purpose of designing protection scheme, is reported to 

have been carried out in order to investigate the effect of 
* « * 

six-phase line in comparision to double-circuit three-phase 
lines with equal line to line (i.e. adjacent phase) voltage, 
by deriving appropriate six-phase transformations, such as 
symmetrical component transformation and Clarke* s component 
transformation. 

It is not only the multi-phase transmission (high 
phase order transmission i.e, HPOT) system, but also the 
six-phase (Multi-phase) synchronous machines have been examined 



for higher rating* employing different types of stator wind- 
ing arrangements since a long time. The steady-state analysis 
of a six-phase machine is reported to have been studied with 
the help of an orthogonal transformation 1 * conceptualizing it 
as a combination of two sets of Park* s transformation. 

Although a considerable amount of work on the multi- 

I 

phase transmission and generation system has been reported 
in the literature* yet the other aspects of investigations* 
such as transient and dynamic stability studies including 
modelling for the purpose of such studies have not been under- 
taken in detail so far. These aspects have been dealt in 
detail in the present work. 

The salient features of the work reported in the 
thesis are as follows: 

1. to develop a detailed dynamical modelling of a six-phase 
synchronous machine; 

2. to develop mathematical representations for various multi- 
phase (n-phase} elements* such fes n-phase synchronous 
machine* 3-phasa/n— phase transformers'* n-phase transmission 
lines* of power system networks for steady-state condition; 

3. to investigate the dynamic stability of a six-phase syn- 
chronous machine and its corrparision with that of three- 
phase machine; and 



4. io investigate in detail the load flow, transient stabi- 
lity and dynamic stability with a view to examine the 
performances of mixed-phase (3-phase and 6-phase) systems 
and completely six-phase systems, with that of conven- 
tional three-phase system. 

A brief account of the work reported in the thesis 
is presented in the following paragraphs. 

The dynamical equations of a six-phase synchronous 
machine have been developed with the help of an orthogonal 
transformation starting with the inductance calcul actions of 
the machine right from the fundamentals. Further, a linearized 
model in terms of current variables has been developed, and 
employing the eigenvalue technique, the dynamic stability of 
a six-phase machine has been investigated and compared to 
that of a three-phase synchronous machine. 

Generalized mathematical representations of the 
various components of multi-phase (n-phase) power system 
elements for steady— state condition have been developed. A 
mathematical model of a multi-phase (n-phase) synchronous 
machine has also been developed in phasor co-ordinates from 
which the model of a six-phase machine has been derived. The 
mathematical models of 3-phase/n-phase transformers for various 
connections', viz. 3-phase star/n-phase star, 3-phase delta/ 



n-phase star* have been developed with tappings on the pri- 
mary as well as on secondary side taking leakage impedance 
of the transformer into considerations. The three-phase 
equivalent of a n-phase transmission line has been developed 
in terms of impedance and ABCD parameters for it representa- 
tion. This representation is helpful in analysing the system 
on completely three-phase basis. The equivalent single-phase 
representation has also been obtained with the help of the 
three-phase equivalent in terms of ABCD parameters represent- 
ing Tt -circuit to analyse a fully balanced system. The n-phase 
equivalent of a three-phase line connected with a 3-phase/ 
n-phase transformer at each end of the line* has been obtained. 
This n-phase equivalent representation will be needed to carry- 
out the analysis of the mixed system completely on n-phase 
basis (i.e. if the interest of investigation lies on the n- 
phase part of the system) . 

The load flow study has been carried out to investi- 
gate the impact of converting double circuit three-phase lines 
to single circuit six-phase or twelve phase lines under two 
voltage conditions* namely 1) when the phase to ground voltages 
of the multi-phase line and the double circuit three-phase 
lines are equal and 2) when the phase to ground voltage of the 
multi-phase line is V 3 times that of the double circuit three- 
phase lines. Two types of cases have been taken into consi- 



derations: 1) for a given slack bus voltage; and 2) for a 
given reactive generation at the slack bus. The single-phase 
representation of a n— phase transmission line connected via 
two 3-phase/n-phase transformers has been obtained on the 
same base quantities as three-phase transmission line. 

The direct and quadrature axis transient and sub- 
transient inductances of a six-phase (multi-phase) synchro- 
nous machine have been calculated from the fundamental. The 
various time constants of a six-phase synchronous machine 
have also been obtained. The transient stability investiga- 
tion has been carried out on a sample network for the various 
system configurations and voltage conditions. In this analy- 
sis, the synchronous machine is represented as a constant 
voltage source behind a direct-axis transient reactance, and 
the loads as constant impedances. The mixed 3-phase and six- 
phase system, and a completely six-phase system have been 
investigated and cornpa red with conventional three-phase sys- 
tem from the transient stability point of view. 

The dynamic stability of a sample network for its 
different system configurations and voltage conditions has 
been investigated employing eigenvalue technique based on a 
linearized model. The two-axis model and the classical model 
of the synchronous machine have been used. The performance of 



a mixed 3-phase and six-phase system with the different system 
loadings has been compared to that of a completely three-phase 
system. Also, the dynamic stability of a completely six-phase 
system has been examined with the same system loadings as the 
conventional three-phase system and also with the increased 


system loadings 



CHAPTER 1 


INTRODUCTION 

1.1 INTRODUCTORY REMARKS 

Few people could have really thought about the 
electricity at the time of its invention in the last century 
for its presently holding such a position of prime irrportance 
as a vital means of serving people, especially towards their 
prosperity because of its manifold applications in the deve- 
lopment. The electricity has now been exploited for its use 
In diversity of applications including domestic* agricultural, 
commercial and industrial fields ranging from smaller to 
bigger ones which essentially enhance the industrial and 
commercial growth with which the individual economy cannot 
remain unaffected. Although it was initially available for 
the use of limited localities, and that too in the specific 
fields, it is now meeting the growing demands with Increasing 
generation. Moreover, electricity is now gradually becoming 
a kind of necessity. Then, it is but natural to arise a sort 
of inquisitiveness for its coming to such a stage, and therefore 
it is necessary to have a look on the growth of electric power, 
which is mentioned below in brief. 

The electricity, an English word, was first used by 
Sir Thomas Browne in. 1646, but William Gilbert Coined in his 
publication * De Magnete 1 in 1600, the word 'electric©* for the 
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substances which behave like amber (elektron in Greek) which 
when rubbed with a dry cloth produces static electricity — 
an idea already enunciated by Thales of Miletus as back as 
6th century B.C. Moreover/ Gilbert pointed out the difference 
between the attraction of artber and that of load stone/ latter 
only attracts iron. But the continuous source of electricity 
was not available till Volta invented Battery in 1800 # and it 
remained the only source till the spate of inventions and 
developments including the discovery of the relation between 
electric current and magnetic field in 18 31, the electro- 
magnetic laws of induction and the construction of an electric 
motor by Faraday in the same year, the first generator and 
commutator by Pixii in 1832 and 1833 respectively, paved the 
way for their commercial exploitation by Gramme in 1870. 

Near about 1851, there was a great interest in apply— 
ing electricity as the source of illumination for the light 
houses and this gave a considerable impetus to the design of 
magneto-electric machines including that of Saxton, Sturgeon and 
Stohrer. In 1857 two dc machines with 90 rpm were installed 
in south fourland lighthouse and the electric light was thrown 
on the sea for the first time in 1858. In 1873 Gramme installed 
his first machine for arc lighting to replace the old generator 
in the french lighthouse, and at the same time, he build several 
generators to replace batteries. In 1878 Edison devised a 
platinum wire filaments, and by 1882 he installed over 150 plants 
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furnishing power for 3000 incandescent lamps from Edison 1 s 
Jumbo dc dynamos driven by steam engines. 

AC systems were being developed in Europe during 
1880s / one of the most successful was that of Lucien Gauland 
and J.D. Gibbs of Paris who first demonstrated their system in 
London in 1881 , and in the same year Savoy theatres Grossnover 
Gallery with 2-kV ac machine was equipped. In USA, the first 
experimental ac distribution system was installed in 1885 which 
supplied 150 lamps in the town of Great Barrington, Massachusettes 
The first single-phase ac transmission line was put into opera- 
tion in 1890 to transmit electric energy generated by water power 

to a distance of 13 miles from Willamette falls to Portland 

1 

(oregon), USA, operating at 3.3 kV with a frequency of 33-^ c/s. 

In 1891 the three-phase long distance ac transmission 
line from Laufen to the Frankfurt electrotechnical exhibition 
in Germany at 30-kV for a length of 109 miles was a convincing 
and successful demonstration. The two-phase ac distribution 
system was demonstrated to the public at the Columbian Exposition 
in Chicago in 1893 where the first Large Westinghouse alternator 
at a frequency of 60 c/s was shown and thereafter, the three- 
phase ac transmission gradually replaced dc system. In 1894., 
there were five polyphase generating plants in USA of which one 
was two-phase and the remaining four three-phase. In Great 
Britain the earliest three-phase generation and transmission was 
installed in 1900 at Wood Lane Station transmitting power at 
6.6-kV. 



During 1906-1911 the transmission voltage in USA 
mounted from 13-kV to 150-kV. In 1921 the first 220-kV came 
into operation. In 1936 a transmission line with 287 *-kV from 
Hoover Dam to Los Angles was installed. The first 345-kV line 
in USA was commissioned in 1953; 500-kV line in 1965; 765— kV 

line in 1969 and the transmission lines- upto 1500-kV is consi- 
dered for future, though feasibility studies have been con- 
ducted upto 1100— kV only by BPA (Oregon) , USA. 

Electric systems were virtually operated till 1917 
as individual units because they started as isolated systems 
and spread out only gradually to cover the wider area. The 
demand for large blocks of power and increased reliability 
suggested the interconnection of neighbouring systems. 

As the function of a power system is to generate 
electrical energy economically and without the ecological dis- 
turbance, and to transmit this energy over transmission lines 
and distribution networks with the maximum efficiency and 
reliability for delivery to consumers located at far peaces at 
virtually fixed voltage and frequency, the research regarding 
this went through various stages to achieve this objective with 
the investigations examining mainly different aspects concern- 
ing the subnetworks, viz. generation, transmission and distri- 
bution subnetworks. But this is not very easy because of the 
rising number of users on one hand, and the technical and the 
economic limitations on the other* However, the efforts have 



gone to increase the generating capacity and simultaneously 
transmit the power to the point (i.e., place) of utilization. 
The 3-phase generation and transmission is already established 
through the research and development and is working satisfac- 
torily. This is because of various advantages offered by a 
three-phase over a single-phase ac system, such as a single 
three-phase synchronous generator occupies less floor areas, 
costs less and also needs much less maintenance in cornparision 
to three single-phase synchronous generators, and three-phase 
transmission is certainly more economical, efficient and reli- 
able compared to single-phase transmission system. However, 
with the growing size and complexity of the modern power 
systems, the main problem for the electric utility has been to 
increase the power transfer capability of the transmission 
lines so as to cater to the needs at the load centres far off 
from the generating stations. The electric utilities have 
found its solution by increasing the transmission voltage of 
the three-phase transmission line. But the uprating of the 
three-phase transmission voltage also seems to have reached a 
point of saturation because of its inherent drawbacks, such as 
increased audible noise level, increased voltage gradient on 
the conductor surface and the strong electric field on the 
ground with possible bilogical effects, when it is viewed as a 
means of transmitting more power at EHV and UHV levels and also 
it has difficulties in acquiring new rights of way, esp ecially 
in urban and conjested areas even at HV and EHV levels. 



Some of the common concern of electric utilities 


have been to economise the use of land and to increase the 
acceptability of overhead transmission lines from the points 
of visibility and appearance. These problems have been tried 
to be solved through improved design methods, namely compac- 
tion, aesthetic designs, etc., aiming at substantial reduction 
of sizes of three-phase high-voltage lines. But at the same 
time, the limitations in acquiring new rights-of-way attracted 
the ■ attention of power system engineers towards searching 
its viable alternative that ultimately resulted in the idea of 
a multi -phase transmission system. This topic of multi-phase 
i.e. high-phase order (HPO) transmission system has been taken 
up for the investigations, and the results arising out of 
several case studies have revealed its considerable advantages 
over the conventional three-phase transmission system, when 
the two systems, namely three-phase and six-phase transmission 
systems are considered for increasing the power transfer capa- 
bility. Although the feasibility studies of multi-phase trans- 
mission system have been carried out and the research are in 
progress; impact of the conversion of a double-circuit three- 
phase line to a single— circuit &lx-— phase line has been studied, 
However, it requires the in-depth investigations of transient 
and dynamic stability analyses of the power system consisting 
of multi-phase transmission lines* 
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From the generation point of view, the efforts 
towards the larger rating of a synchronous machine have been 
made as it is justified from the security and economic consi- 
derations', for it reduces the operational cost per unit of 
output. To generate more power from a single unit, a multi- 
phase stator winding has been envisaged, as it occupies less 
floor area and needs less maintenance compared to a traditional 
three-phase synchronous machine. 

The advantages of a six-phase transmission system 
over the conventional three-phase transmission system reported 
so far, coupled with the idea of six-phase transmission lines 
having the same number of conductors as those of double-circuit 
three-phase lines, and an inquisitiveness towards the querry 
whether a six-phase transmission line will really serve the 
purpose of increasing the transmission capability of overhead 
lines in a better way than a conventional three-phase line along 
with an intention to study the performance of a six-phase synch- 
ronous machine with respect to that of a conventional three-phase 
machine have indeed provided us with a motive to find out the 
answers to the quarries, and as a result, this topic of multi- 
phase system has been taken for its transient and dynamic 
stability investigations, which are relatively; untouched, along 
with appropriate mathematical modelling. 

The power transfer capability of a transmission line, 
if it is viewed with respect to a double-circuit three-phase 



line, can be increased either by increasing the transmission 
voltage of a three-phase line or by converting the double- 
circuit three-phase line to a single -circuit six-phase line 
with equal line to line (i.e., adjacent phase) voltage which 
will be better out of the two options can be adjudged if the 
analysis is carried out on both the systems. The steady- 
state analysis has revealed several advantages of the latter 
option (i.e. single-circuit six-phase line) with respect to 
the former. It is important to observe their performances 
from the stability point of view as well. 

Indeed, the stability is a single entity which 
expresses the ability of the machines to keep the relative 
angular positions of the rotor constant, if any sort of dis- 
turbances is experienced by the system. In power systems, 
two types of stability are encountered, namely the dynamic 
stability and the transient stability, depending upon the 
magnitude and the duration of the disturbance. The power 
system which is practically a very large and complex networks 
always experiences small disturbances; the stability study 
caused by such disturbances is termed as dynamic or operating 
point stability study. The other, transient stability study, 
is carried out when the system experiences large disturbances 
such as faults on a high voltage transmission line, etc., for 


a short duration, 
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For the investigation of any system, it is essen- 
tial to develop the suitable mathematical model as a realistic 
analytical tool to evaluate the system performance; so is true 
with the multi-phase systems, and the present work is an 
attempt in this direction by developing suitable mathematical 
models to carry out the dynamic and transient stability investi- 
gation of multi— phase systems in order to observe its perfor- 
mance in comparision to that of three-phase systems. 

Before the work carried out in this thesis is pre- 
sented, a brief account of the various aspects of investiga- 
tions on the multi-phase systems reported so far in the lite- 
rature has been highlighted under the following section. 

1 . 2 LITERATURE SURVEY 

An idea of multi-phase (more than three phases) 
transmission system i.e., high phase order transmission (HPOT) 
system was first presented byBarmes and Barthold[l] by in- 
vestigating various aspects, such as voltage gradient on con- 
ductor surface which is important from the viewpoints of corona 
inception, electric field intensity on the ground level from 
environmental viewpoints which is helpful in designing EHV and 
UHV lines, space energy flow distribution around a conductor 
which is essential in studying the efficiency of spade utiliz- 
ation for power transmission systems, on the different networks 
of high-phase order ranging from three to thirty six phases 
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with conductors symmetrically arranged in a circular configu- 
ration in order to carry out a comparative study of the diffe- 
rent phase order transmission! system. It is found that the 
voltage gradient on the conductor surface for the six— phase 
system is less than that for the three-phase system while the 
field intensity on the ground level of the six-phase line is 
higher. As indicated, the energy flow distribution around 
conductors is an important index for the effective utilization 
of rights-of-way, it has been observed that the high-phase 
order transmission (HPOT) is better than three-phase transmiss- 
ion system in the effective and efficient utilization of rights- 
of-way; also, it can be used to transmit more power at the 
increased transmission voltage with the same rights-of-way. 

It thus attracted the attention of researchers and utility 
engineers and consequently, the investigations are being carr- 
ied out since then in order to study further the high-phase 
system as an alternative to a conventional three-phase system* 

A significant amount of research efforts has been 
applied in investigating the feasibility of multi-phase trans- 
mission systems since then, and the conclusions derived from 
the various feasibility studies have been reported in the 
literature [ 1-12 ] . a systematic and an elaborate study on 
multi-phase system was carried out under the joint efforts of 
Allegheny power system and West Virginia University (APS-WVU) , 
and later on, under a project sponsored by the US Department of 
Energy. 
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Venkata et al [ 2 ] have outlined the concept of multi- 
phase transmission systems in terms of six-phase lines and 
have carried out the study of load flow, stochastic reliability 
and stability analyses by making a preliminary investigation on 
a 138-kV system under two configurations , namely 1) a completely 
three-phase system having one double-circuit three-phase line and 
2) a combined three-phase and six-phase system by converting the 
double-circuit three-phase line to a single -circuit six-phase 
line with 138-kV (phase to ground).. It has been shown from the 
load flow analysis that the more demands can be met from confi- 
gurationi 2 with improved voltage regulation and higher trans- 
mission efficiency without overloading the line. Further, two 
ways of increasing the power transfer capability of a double- 
circuit three-phase line have been considered: a) increasing the 
transmission line voltage of the double-circuit three-phase 
line by)/"3j b) converting the double-circuit three-phase line 
to a single-circuit six-phase line, the latter having y 3 times 
the phase to ground voltage than the former. Although both of 
these alternatives are electrically ideal, the six-phase trans- 
mission line has been shown to require less transmission 
corroidor[ 2 ] . 

Bhat et al [ 3 ] have derived symmetrical component 
transformation for a six-phase system as a straight forward 
extension of the three-phase system, and have carried out the 
fault analysis to determine the fault currents and voltages on 
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a six-phase transmission system. The feasibility study on the 
multi-phase system has been carried out by Guyker et al [ 4] and 
various aspects have been taken for investigation* and especia- 
lly, an idea of EHV six-phase line as a viable alternative to 
UHV three-phase line has been presented. 

Stewart and Wilson in their companion papers [5,6] 
have systematically examined the multi— phase transmission sys- 
tems in detail. In their first paper [5] , the feasibility 
analysis has been carried out under steady-state condition on 
two transmission voltages, for example, 80-kV and 289-kV (phase 
to ground) for six-phase and twelve phase transmission lines in 
order to compare them with 138-kV and 500-kV (phase to phase) 
double-circuit three-phase lines with a view to multi-phase 
system voltage nomenclature, the effect of transmission line 
impedance on the steady-state operation, the line and the gene- 
rator current unbalances and the electrical environmental aspects, 
such as radio noise, audible noise and the electric fields at 
the ground; while in their second paper [6] , an extensive 
switching surge studies arising out of fault clearing on the 
six-phase lines have been taken for investigation to analyse 
the fault overvoltages, overvoltage due to interphase coupling 
and lightning performance to study several aspects of line 
design, such as the insulation level, the conductor clearances, 
etc. From the study, it has been shown that the phase to ground 
voltage is appropriate for system voltage definition for multi— 
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phase system; the thermal line loading increases proportional 
to the number of phases while surge impedance loading (SIL) 
increases at the lesser rate; the generator current and the 
line current unbalances due to six-phase transmission line 
with or without transposition are comparable to, or better 
than those of double-circuit three-phase line; radio and audi- 
ble noise level are reduced with the increase in number of 
phases; the electric field at the conductor surface decreases 
with increasing phase order while the field at the ground level 
increases; fault overvoltages for a six-phase line can be 
slightly higher than those for a comparable three-phase line 
while for the phase order higher than six, the fault overvol- 
tage magnitudes are similar to those of a three-phase line. 

The interphase coupling has been shown to cause high overvoltages 
on the open phase when a single-phase switching is applied to 
high-phase order (HPO) systems with shunt compensation. Further, 
the phase to phase switching surges have been shown to assume 
primary importance in the design of high-phase order systems 
thereby requiring a surge control especially for compact phase 
spacing. The rate of rise of recovery voltage (RRRV) during 
fault clearance is less for high-phase order systems than that for 
three-phase systems with comparable short-circuit MVA. The 
lightning performance has been shown to be comparable to the 
three-phase system. 
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The terminal insulation level for a multi-phase 
system has been shown to be somewhat higher. However/ the 
typical insulation levels, even for a six-phase system does 
not exceed to those presently applied on three-phase systems 
of the same phase to ground voltage. 

The fault analysis of a six-phase transmission 
system has been carried out by Onogi and Okumoto [ 7 ] with the 
help of the combined use of two-phase and three-phase symmet- 
rical component method, and further a very interesting method . — 
connecting two single -phase transformers or one single-phase 
three-legged transformer for each phase — has been proposed 
for suppressing the fault current. The fault current has been 
shown to be suppressed by impedances connected to the tertiary 
winding, and if the tertiary impedances are infinitely large, 
the line to ground voltage of the sound phase has been observed 
to be twice as high as the prefault line to ground voltage and 
the fault current (either ground fault current or short circuit 
current) is as small as the exciting current. Two schemes have 
also been proposed to suppress the overvoltages. 

Willems [8] has developed the mathematical model 
for a six— phase line by ABCD parameters similar to that derived 
for a three-phase system by obtaining the mathematical descrip- 
tion of the various three-phase/six-phase transformers in the 
phasor co-ordinates without considering their leakage impdeances. 
Employing these models, the three-phase equivalent represents- 
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tion of a six-phase line both in impedance and admittance 
forms and its equivalent ABCD parameters have been derived to 
carry out the analysis of a composite three-phase and six- 
phase systems on a three-phase basis. In order to analyse 
the system completely on six-phase basis, the equivalent six- 
phase of a three-phase line has been obtained. Further, the 
equivalent single-phase of the composite system has been 
derived to carry out the analysis on single-phase basis. 

A six-phase transmission line simulator has been 
developed by Chinnarao and Venkata [9] to study the multi- 
phase transmission system. Peeran et al [lo] have carried out 
the fault analysis on a six-phase system with the help of Alpha- 
Beta -Zero component transformation. 

Willems [ 11 ] has derived the symmetrical component 
transformation of a six-phase system by considering it as the 
two coupled three-phase systems, and this concept can be exten- 
ded to develop transformation for the multi-phase systems that 
are multiple of three. By this derived component transformation, 
the impedance matrix of a six-phase transmission* line can be 
diagonalized. 

Grant et al [ 1 2 3 have constructed experimental six- 
and twelve-phase lines achieving a considerable compaction, and 
choosing a voltage of 80-kV phase to ground (corresponding to 
138-kV phase to pha.se for a three-phase circuit) the electrical 
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and mechanical tests have been performed demonstrating practi- 
cal insulator and hardware designs to ascertain the correct- 
ness of analytical procedures developed for a high-phase order. 

The use of symmetries inherent in power system net- 
works and its exploitation for the purpose of simplifying 
their analyses was examined by Singh et al [ 13 ] > and based 
upon the symmetry considerations alone and using group the- 
oretic techniques, the symmetrical and Clarke's component 
transformations for the multi-phase systems have been cons- 
tructed from the fundamentals [l4 ] . These transformations 
were further extended for the generalized N-Port network [ 15] . 

The fault analysis on a mixed three-phase and six- 
phase system using symmetrical component transformation has 
been carried out by Nanda et al [l6j . Eleven types of faults 
on the six-phase line have been considered for investigations 
after obtaining voltage and current relation as well as sequence 
network connections corresponding to each type of fault. The 
six-phase and the double -circuit three-phase transmission lines 
have been compared under two conditions, namely the most severe 
all phases to ground fault and the most common single-phase to 
ground fault. It has been found that the fault current invol- 
ving all phases in the six-phase line is slightly higher than 
that for the three-phase line, whereas the fault current in 
case of a single-line to ground fault of the six-phase line is 
nearly half to that of the three-phase line. 
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Kallaur and Stewart [ 1 7 ] have examined the uprating 
to 230-kV (230-kV phase to phase for three-phase operation and 
2 30/ ]f3 - kV phase to ground for a six-phase operation) by- 
selecting five double-circuit three-phase lines. The study has 
revealed that, among the five specific line designs, one is 
technically suitable for 230-kV three-phase operation and the 
three are for six-phase operation while the remaining one is 
unsuitable for either case. Further, it has been concluded 
that the phase conversion required for six-phase results in 
higher terminal cost. However, the benefits of the reduced 
losses and higher power handling capability may offset more than 
the cost penality. 

The steady-state characteristics viz., voltage gra- 
dient on the conductor surface, ground level electric field 
intensity, space distribution of energy flow density of a six- 
phase transmission line has been investigated by Takasaki et 
al [18] . A single-circuit six-phase line with sixteen phase 
arrangements has been compared to a double circuit three-phase 
line with two phase arrangements. The study reveals that the 
field intensity on the ground level in case of a six-phase 
line is higher than that of a double-circuit three-phase line, 
while the electric intensity on the conductor surface is lower. 

It is very interesting to observe that the phase arrangement 
which gives the smaller voltage gradient on the conductor surface 
gives a higher field intensity on the ground level. 
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Stewart and Grant [ 19 ] have carried out the studies on 
the six— phase and twelve— phase test lines, and have confirmed 
the analytical predictions of electrical and mechanical beha- 
viours demonstrating a simple design and attractive structures 
built with standard hardware. 

Venkata et al '[ 2o] have carried out the fault ana- 
lysis on a 138-kV six-phase line by discussing the various 
types of fault and their different 23 types of significant 
combinations likely to occur on a six-phase line, and have 
developed the analytical expressions for all the 23 types of 
fault. 

Tiwari et al [21] have developed the mathematical 
model of a six-phase synchronous machine on the same line as 
that of the three-phase synchronous machine [22] . The mathe- 

1 

matical models for the various connections of three-phase/six- 
phase transformers taking leakage impedances into considerat- 
ions have been developed in the phasor co-ordinates. Further, 
the models of a six— phase line integrated with a three-phase 
wye/six-phase star connected transformer at the each end of 
the line have been derived in impedance and admittance forms. 
Moreover, the three-phase equivalent and its ABCD parameters 
have been developed and from this, the single-phase equivalent 
has been systematically derived [21] . The six-phase equivalent 
of a three-phase transmission line has also been obtained [23] . 
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Although the fault analysis and the load flow ana- 
lysis for a three-phase system have been reported [ 24-29] in 
the phase frame of reference, the load flow investigation on 
a six-phase line has been carried out by Tiwari et [ 30, 31 ] 
for the balanced system on the single-phase basis [ 30] and 
for the unbalanced system in the phasor co-ordinates [31] . 

In the investigation, taking the phase voltage of the six- 
phase line equal to V"3 times that the three-phase line, 
the same conclusions have been derived as observed in earlier 
paper [ 2] , when a double-circuit three-phase line is conver- 
ted to a single-circuit six-phase line. 

Tiwari et al [ 31 ] have developed six-phase equivalent 
of a three-phase line to facilitate the analysis on the six- 
phase basis and have carried out the fault analysis by utili- 
zing the representation of three-phase elements in the phasor 
co-ordinates. 

The impact of converting a do ubl e-circuit three- 
phase line to a single circuit- six -phase line has been 
studied while carrying out the fault analysis on the combined 
three-phase end six-phase systems. 

Grant and Stewart [32] have examined an alternative 
to UHV three-phase transmission in which 462_kV six and twelve 
phase lines and also 317-kV twelve-phase lines have been 
compared with 1200— kV UHV three-phase transmission of compa- 
rable power transfer capability. It has been found that HPO 
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lines are smaller and require less rights-of-way than their 
three-phase counterparts. Stewart et al [33] have demons- 
trated that, for lines longer than a few miles, saving in the 
line offset the terminal equipment costs and makes HPO an 
economic alternative. Terminal equipment costs for HPO sys- 
tems are usually slightly higher than those for UHV three- 
phase system. 

Wilson et al [34] have compared a 1200— kV, three- 
phase single-circuit line with, a 462-kV six-phase and twelve- 
phase single— circuit lines and have observed that the insulation 
design complexities increase with the number of phases. For 
a circular conductor array, the most critical phase to phase 
spacing due to switching surges has been shown to be normally 
between adjacent phases. Switching surges on six-phase line 
are in the same range, or slightly higher than those on a 
three-phase line. If the switching resistors are incorporated 
with the circuit breaker, surges on the six-phase line are 
less than those on the three-phase line. Switching surge 
magnitudes on six-phase systems for practical combinations of 
variables, result in transmission line insulation and clearance 
requirements comparable to those on three-phase lines. 

Further, the transients associated with six-phase 
capacitor switching have been studied by Ramaswami et al [ 35 ] 
on a 138— kV six-phase system by first deriving the analytical 
expressions and then by getting the results verified by computer 
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programming. The impact of restrike in a six-phase system 
has been found to be less than that in a three-phase system. 

Guyker and Sh ankle [36] have presented the study 
of the terminal equipment design and the results therefrom 
for the costs of 138-k.V double -circuit three-phase line to 
a 138— kV six-phase or 230-kV double-circuit three-phase 
operations, A comparision of the total costs shows that the 
cost of the 138-kV six-phase uprating option is higher than 
that of the 230-kV three-phase uprating option due to the 
transformer cost. However, 138-kV six-phase operation does 
offer some reduced radio interference and audible noise 
compared with uprating to 230-kV double-circuit three-phase 
operation. 

Gross and Thompson [ 37 ] have presented a method for 
calculating phase and sequence voltage and currents through- 
out a power system of mixed phase order subj ected to an arbi- 
trary fault type at an arbitrary bus in the system, and have 
described the procedure on a power system network having twelve 
phase transmission line. 

The fault analysis on a twelve-phase transmission 
line has been carried out by Pal and Singh [38] , and 23 types 

of faults are analysed with the help of symmetrical component 
transformation for the twelve-phase system. From the study of 
the 12 -phase to ground fault and the single-line to ground 
fault, it has been shown that the fault level for a twelve- 
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phase system is less than that for a three-phase system. It 
is observed that, the faults other than the twelve-phase to. 
ground, are dual in pair, for example, single-phase to ground 
and eleven-phase to ground, as the sequence networks are 
connected in parallel fob eleven— phase to ground fault while 
they are connected in series in the case of single-phase to 
ground fault, similarly ten-phase to ground fault and two-phase 
to ground fault and so on. 

Swarup et al [39] have carried out the load flow 
investigation on a combined three-phase and twelve-phase 
transmission lines under balanced as well as under unbalanced 
condition by deriving the models of power system elements 
required for such studies. 

Tiwari and Singh [ 40 ] have carried out the transient 
stability analysis on a mixed three-phase and six-phase system, 
and have investigated the impact of converting the double- 
circuit three-phase lines to single-circuit six-phase lines 
with the same adjacent phase voltage, and therefrom the mixed 
transmission system has been shown to be better than the 
completely three-phase system from stability point of view. 

Saleh and Laughton [ 41 ] have developed mathematical 
representations for the various three-phase/six-phase trans- 
formers in the phasor co-ordinates. 
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The efforts towards the design and fabrication for 
the larger rating of a synchronous machine have .been made 
since a long time [ 42-44 ] But the attempt was initially 
hampered by certain limitations in the circuit breaker inte- 
rrupting capacity and the enormous reactors needed to limit 
the fault current. To alleviate the difficulty, a double- 
winding generator was suggested [43] in which each of the 
two separate stator windings is connected to a different 
sections of the station bus; the generated voltage in the two 
windings are equal in magnitude and are in phase. Several 
types of windings have been considered [44] . These machines 
have been found adequate until the fault current problem is 
alleviated by connecting each generator to its own step-up 
transformer. However, the problem of enlarging the generator 
rating accrues from the limitations arising out of stator 
voltage, mechanical forces on the stator coils, stray-surface 
losses and the physical strength of the rotor forging and the 
major coil retaining rings. 

The problem of the coil forces has been solved with 
the ordinary double-winding [45] . However, the unsymmetri- 
cal winding causes additional losses and create difficulties 
in the end winding arrangement. To avoid these problems, two 
separate stator windings displaced at 30 electrical degrees 
have been suggested [46] . 



p, Robert et al [ 46] have shown that the stray losses 
due to phase harmonics can be materially reduced by means of 
two separate stator windings displaced at 30 electrical degrees. 
The methods of protection of multi-phase generator against 
faults have also been dealt with [47] , 

Holley et al [ 48] have shown that electromagnetic 
forces on the stator windings of large turbine generators can 
be substantially reduced by using the winding arrangements (two 
three-phase stator windings displaced at 30 electrical degrees) 
that permit an increased number of stator slots. This genera- 
tor has been shown to have a very low harmonic contents. It 
has been found to be particularly attractive for feeding power 
to a dc transmission grid because of their reduced sensitivity 
to harmonics in the line current that might be introduced by 
the rectifier system. 

Notwithstanding the great strides being made, two- 
pole generators appear to be reaching a plateau, at about 
1500 MVA [49] . Fuchs and Rosenberg [ 50 ] have carried out the 
steady-state analysis of a six-phase synchronous generator by 
conceptualizing it as the combination of two three-phase stator 
windings with 30 degree displacement by means of an orthogonal 
transformation which has been taken as a set of two three-phase 
Park 1 s transformation in dqo variables. 

The use of a six-phase synchronous generator along 
with its associated six-phase/ three-phase transformers for 
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high power applications have been discussed by Hanna et al 
[ 51 ]. 

The six-phase synchronous machine has been a alysed 
by Schiferl and Ong in their companion papers [ 52,53 ] . In 
their first paper [52] a harmonic phasor method for steady- 
state analysis of a six-phase machine with mixed AC-DC stator 
connection has been examined, and in their second paper [53] , 
a novel application of the six-phase machine in an uninterrupt- 
ible power supply scheme has been proposed. 

Willems [54] has indicated that the six-phase gene- 
rator and its associated transformers can be treated for fault 
analysis employing symmetrical component transformation. 
Further, a three-phase equivalent representation has been 
obtained by combining the generator voltage and the impedances 
with the equivalent circuit of the transformer. 

After this brief literature survey, a summary of the 
work carried out in this thesis and its chapterwise discription 
are presented. 

1.3 OBJECTIVE AND SCOPE OF WORK REPORTED IN THE THESIS 

The objectives of this thesis are set forth as 
follows : 

1. io develop a dynamic model of a six-phase synchronous 

machine in order to carry out the dynamic stability ana- 
lysis of the synchronous machine and also to compare the 



stability performance of a six-phase machine to that of a 
conventional three-phase synchronous machine; 

2. to develop suitable mathematical models of various multi- 
phase elements of power system networks, viz., synchronous 
machines, transformers and transmission lines for a gene- 
ral n— phase (order) system under steady-state conditions; 

3. to carryout the transient stability investigations employ- 
ing the mathematical models in order to study the multi- 
phase system especially the impact of converting double- 
circuit three-phase lines to single circuit multi-phase 
lines * 

4.. to study the dynamic stability on a combined three-phase 
and multi-phase transmission lines with state-space 
formulation j and 

5.. to study the performance and feasibility of six-phase 
systems with the aid of several case studies. 

1.4 CHAPTERWISE DESCRIPTIONS 

A chapterwise description of the work carried out 
in the thesis is presented as follows. 

Chapter 2 starts with the inductance calculations 
of the six— phase synchronous machine from the fundamentals. 

The dynamical equations of the machine have been developed 
with the help of a generalized park' s transformation. Further, 


a linearized model of the machine has been obtained in terms 
of current variables. The six-phase (multi-phase) synchro- 
nous machine has been analysed for the different conditions, 
and the performances of the machine has been compared to those 
of the three-phase machine from the point of dynamic stability. 

Chapter 3 presents generalized mathematical models 
of the various components of n-phase power system elements 
under steady-state condition. A mathematical model of multi- 
phase synchronous machine in generalized form has been derived, 
and from which that of the six-phase machine has been obtained. 
A three-phase/n-phase transformer will be required when n-phase 
line is connected to a three-phase line and therefore, it 
requires a model for such transformers. Therefore, the models 
for three-phase/n-phase transformers with tappings for various 
connections, such as three-phase star/n -phase star, three-phase 
delta/n-phase star, etc., have been developed.. The three-phase 
equivalent of a combined n-phase transmission line and three- 
phase/n-phase transformers have been developed in terms of ABCD 
parameters representing U -circuit. The n-phase equivalent of 
a three-phase transmission line with an appropriate transformer 
at each end of the line has been obtained. These equivalent 
representations are required to carry out the analysis of the 
overall system entirely either on 3-phase or on n-phase basis. 
Further, a single-phase equivalent of the n-phase transmission 
line with an appropriate transformer at the each end of the 
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line has been derived from its three-phase equivalent in terms 
of ABCD parameters representing n -circuit. This representation 
is adequate to analyse the entire system, if it is balanced 
and symmetric, on a single-phase basis.. 

In chapter 4, the load flow investigation has been 
carried out on a mixed three-phase and six-phase lines and 
also on a mixed three-phase twelve-phase lines to study the 
performance of the mixed power system in corrparision to that 
of a completely three-phase system for the various configura- 
tion. To study the effect of conversion of double -circuit 
three-phase lines to single circuit multi-phase (six-phase or 
twelve-phase) lines, two voltage conditions have been assumed, 
namely, the phase to ground voltage of a multi-phase line is 
equal or s/3 times that of a double circuit three-phase line. 
Two illustrative examples have been presented for the sake of 
comparision* 

Chapter 5 deals with the transient stability investi- 
gation of the multi— phase systems. The direct axis and the 
quadrature axis transient and sub transient inductances of a 
six-phase synchronous machine have been obtained. The various 
time constants pertaining to the synchronous machine have also 
been determined. Further, the transient stability investiga- 
tion on a sample network has been carried out for the different 
system configurations, for example, a completely three-phase 
system with a few double circuit three-phase line, a combined 



three-phase and six-phase lines where these multi— phase lines 
are obtained by converting double circuit three-phase lines, 
and a completely six-phase (multi-phase) systems. The criti- 
cal clearing lines have been calculated for the different 
fault locations. The performance of multi-phase lines with 
respect to double-circuit three-phase lines have been examined. 

Chapter 6 presents the dynamic stability analysis 
of tile power system with six-phase transmission lines. The 
dynamic model of a power system has been developed in terms of 
voltage and current variables and the dynamic stability has 
been investigated with the help of eigenvalue techniques for 
the different system configurations on a sample network. The 
dynamic stability performance due to the uprating of double- 
circuit three-phase lines to single circuit multi-phase lines 
have been investigated. 

Finally, chapter 7 concludes with the overview of 
the work reported in this thesis and further, outlines briefly 
the future scope of work in this area. 

1.5 CONCLUDING REMARKS 

This chapter presents a brief account of motivation 
for taking up the present work and further, depicts the work 
previously carried out in the area of multi-phase systems. The 
objective and scope of the thesis have been explained briefly 
and the chapterwise work reported in this thesis has been summ- 


arized 



CHAPTER 2 


DYNAMICAL MODELLING AND ANALYSIS OF A MULTIPHASE 
(SIX-PHASE) SYNCHRONOUS MACHINE 

2.1 INTRODUCTION 

In the previous chapter, the work reported in the 
thesis has been summarized along with the literature survey 
on multi-phase power systems. This chapter presents a de- 
tailed dynamic model of a six-phase (multi-phase) synchro- 
nous machine, and also its analysis. As it is intended to 
carry out the dynamic stability investigation of a six-phase 
synchronous machine, a mathematical model of the machine 
under dynamic condition is required to be developed. Such 
models for three-phase synchronous machines are available in 
literature [ 55-59 ] , but for six-phase machines this is al- 
most untouched. 

This chapter initiates in section 2.2.1 with the 
calculation of the inductances of the six-phase (multi-phase) 
machine from the fundamentals. Further, dynamical equations 
of the machine have been derived in section 2.2.3 with the 
help of on orthogonal power invariant transformation presen- 
ted in section 2.2.2. 

In order to investigate the dynamic stability of 
the six-phase (multi-phase) machine, a linearized model in 
state space form has been derived in section 2.3. Further, in 
section 2.4 a case study has been carried out on a sample sys- 
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tem to investigate the dynamic stability of the six-phase 
machine employing eigenvalue technique. In section 2.4.1, 
methods for finding the initial conditions of a six-phase 
synchronous machine have been provided for two boundary 
conditions. Then, with the help of linearized model deve- 
loped in section 2.3 and the initial conditions obtained for 
two boundary conditions from the machine parameters, the 
dynamic stability of a six-phase machine with and without 
darrper windings has been investigated in section 2.4.2 and 

2. 4. . 3 respectively. The results of the investigation have 
been presented and discussed in section 2 ..4.4 and finally, 
the chapter concludes in section 2.5 with the finding of the 
in ve s tiga tion . 

2. . 2 DYNAMICAL MODELLING OP A SIX-PHASE (MULTI-PHASE) SYNCHRO- 

NOUS MACHINE 

In order to develop a detailed dynamic model, a 
six-phase synchronous machine has been considered to consist 
of six-windings (two sets of balanced three-phase winding dis- 
placed at s electrical degrees) on the stator, one field 
winding F on the rotor d-axis, and two fictitious short cir- 
cuited windings on the rotor d and q-axes to represent damping 
effect as shown in Fig. 2,1. One additional (short circuited) 
winding 1 G 1 on the rotor q-axis (not shown in the Fig. 2.1) has 
been included to take eddy current effects into consideration. 
In this section, inductances of the machine windings have been 
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?.2.1 Schematic winding arrangement of a Multi-phase 
(six-phase) synchronous machine . 



calculated from the fundamentals/ and further the dynamic 
model has been derived with the help of an orthogonal trans- 
formation . 

2.2.1 Inductances of a Six— Phase (multi-phase) Synchronous 
Machine 

The machine has been considered to consist of six 
windings namely a,b,c,d,e and f on the stator, one main 
field winding F on the rotor d-axis, the fictitious G winding 
on the q-axis, and two fictitious short-circuited D and Q 
windings on the rotor d- and q-axis respectively. 

Stator Self-inductances : 

Let us consider the flux linkages with any phase i 
of the stator winding. If only phase i is excited, the peak 
amplitude of the mmf wave due to the current in phase i is 

Pi = Vi t 2 - 1 ’ 

where is the effective turns/phase and i^ is the current 
in phase i. 

If F ^ and F ^ are the d-axis and the q-axis compo- 
nents of the mmf, their arrplitudes can be given as 

p di - F i co - £ (e - e i’ (2 - 2) 

and F . = F. sin (© - 0.) 
qi i i 


(2.3) 
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where © is the angle between the reference phase and the 
rotor d-axis; 

between the reference phase and the phase i in con- 
sideration. 


If and P be the permeances along the d-axis and the 

q-axis respectively, the fundamental air gap flux per pole 
along the two axes can be given as 

# ^ = Corr P onent the a i r 9 a P flux along the d-axis 


= F 


di 


p gd = p i V C ° S ( ® " ®i > 


(2.4) 


and &gqi = component of the air gap flux along the q-axis 


= F 


. P = F. P sin (0 - 0.) 
qi gq i gq i' 


(2.5) 


If gf ^ is the flux linking with the winding of the phase i 
when only phase i is excited, then 


S6 • ■ = $ . cos (® - ® . ) + jzf . sin (© - 0. ) 
gn gdi l gqi l 

( 2 . 6 ) 

From equations (2.4) - (2.6), and using trigonometric identi- 
ties 

P - + P P - — P 

jrf . , = N . i. [ -2X— — 22 + -22_- . 2 2 cos 2 (© - ©. ) ] (2.7) 

Cfll JL m Z Z X 



( 2 . 8 ) 


Hence, 


j . . = N. j6 . ./i. 
gn i gn i 


L ^ 
aaO 


+ L o cos 2 
aa2 . 



where 


L aaO 



CP gd +P gq )/2 ' L aa2 ~ N i ( ^d^gq 5 /2 


(2.9) 


If the flux which does not cross the air gap (i.e. leakage 

flux) is also considered, . the self inductance L. . can be 

n 

expressed as 


L . ■ — L-. . + L . . — L _+ L 0 COS 2 

ii li gii aaO aa2 


(© 


® ± ) 


( 2 . 10 ) 


Equation (2.1o) thus gives in general the self -inductance 
of the stator winding i (for i = a,b,c,d,e, f ) . 


If the six windings on the stator are considered as two sets 
of balanced three-phase winding displaced at e electrical 
degrees (where e = 60° for six-phase machine), the self in-* 
ductances of the winding can be computed from equation (2.10), 
and are as follows (assuming the axis of phase a to be the 


reference 

axis) 








L 

aa 

L aaO 

+ 

L aa2 

cos2© 


(2 , 

,11) 

^bb = 

L aaO 

+ 

L aa2 

COS 

(2# 

_ 2E ) 

(2. 

,12) 

Hi ^ 

CG 

L 

aaO 

■f 

L aa2 

cos 

(2© 

- 240°) 

(2, 

.13) 

L dd 

L 

aaO 

+ 

L aa2 

cos 

(2®-240°-2e ) 

(2. 

.14) 
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L = L _ + L o cos (2© + 240 ) 
ee aaO aa2 


L ff = L aaO + L aa2 cos (2s + 240 ~ 2 C > 


Stator windings mutual inductances i 


(2.15) 

(2.16) 


Let the mutual inductance between the phase i and 
the phase j of the stator windings be L. .. Let . . be the 
air gap flux linking with phase j when only phase i is exci- 
ted. The air gap flux is given by 


0 . . ± si .. cos (® - ®.) + 0 . sin (® - ®f) (2.17) 

gxj gdi j gqi J 

where ©j is the angle between the flux linking phase j and 
the reference phase. 

Substituting the values of j# ^ and ^gq± ^ rom equations 
(2.4) and (2.5) in equation (2.17), we have, 

P'-fP P , — P 

s6 , . = N . I • [ ~ 3 9^ ‘- 2 2 cos (©.-©.)+ g - q cos(2®-s -® )] 

g;Lj xx z -LJ z 1 J 

( 2 . 18 ) 

Then, the mutual inductance between the phases i 
and j due to the air gap flux can be given by - 

L ij - L abo cos <®i- J V +L aa2 008(20-^-®,) C 2 - 1 ^ 

i 

where L abo = L aaO in ec 3> iation (2.9). 

The equation (2.19) gives the mutual inductance between the 
phases i and j and therefore, the mutual inductances between 
the stator windings can be obtained and are given below by 
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using axis of phase a as the reference. 


L ab 

= 

L , 
abo 

cos e + L aa 2 cos e ) 

(2 . 20a) 

L ac 

= 

"^abo 

cos 12o° + L aa 2 cos (2©— 120 °) 

(2.20b) 

L ad 

= 

L abo 

cos (-12o°- e ) +L aa 2 cos (20-12o°- e ) 

(2.20c) 

L 

ae 


L abo 

cos 120°+L aa2 cos (2©+12o°) 

(2.20d) 

L 4= 
af 

= 

L abo 

cos(120°- e )+b aa2 cos (2©+12o°- e ) 

(2. 20e) 

^c 

= 

L . 
abo 

cos (-12o°+ e ) +I * aa 2 cos (2®-12o°- e ) 

(2.2 0 f) 

hod 

= 

L abo 

cos (-120°) +L aa2 cos ( 2 ®_ 120 °— 2e) 

(2.20g> 

hoe 

ss 

L abo 

co s (120 + £ ) +L a 2 cos ( 2®+12o — £ ) 

(2.2oh) 

hof 

= 

L abo 

GO s 1 2 0 -HL ry co s ( 2 © +1 2 0 — 2 E ) 
aa^ 

(2.20i) 

L cd 

= 

L . 
abo 

cos £+ L aa 2 cos (2®_240°- £ ) 

( 2 . 2 0 j ) 

L _ 
ce 

== 

L abo 

cos 240°+L aa 2 GOS 2 ® 

(2.20k) 

L ^ 
cf 

= 

L V. 

abo 

cos(24o°-e )+L aa 2 cos (2®— e) 

(2.201) 

L de 

= 

L abo 

cos (24o°+ e ) +L aa2 cos(2®_ e ) 

(2.20m) 

L df 


L abo 

cos 240° + L aa2 cos (2 ©-2 £ ) 

(2.20n) 

L ef 


L , 
abo 

cos e + L a2 cos (2®+240°- £ ) 

(2 ,20p) 

Stator 

to 

Rotor Mutual Inductances % 



These inductances vary with the position of the 
rotor, since it has salient pole construction. 

The mutual inductance between any phase i of the stator wind- 
ing and the rotor winding can be expressed as 

L Fi “ L i F ' L mF cos < 8 - < 2 - 21 > 


L Di ~ L iD 


L ^ cos (© — © . ) 
iriD x 1 


( 2 . 22 ) 
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J Qi 

L iQ ~ L mQ 

sin 

(0 - ©.) 

1 

(2.23) 

'Gi 

L iG = ~ L mG 

sin 

O 

I 

© 

(2.24) 


where i = a,b,c,d,e,f. 


Rotor Inductances : 

The self inductances of the rotor windings, namely 
Lp F / ^dd' L qq' L qg are Gonstan ' t * as these values do not depend 
on the rotor position because the armature has cylindrical 
structure (only rotor has salient pole construction) . 


The mutual inductances between the rotor windings 

are also constant. and L are non-zero while L_ 

FD QG DQ FQ 

L qd and are zero as the corresponding windings (referred 
to by the capital letter subscripts) are at right angles to 
each other.- 


2.2.2 Generalized Orthogonal Transformation 

In order to develop the dynamic model an orthogonal 
transformation for a n-phase machine is required, as the self 
and the mutual inductances of the machine derived in section 
2.2.1 are time dependent. Such type of orthogonal transfor- 
mation known as Park 1 s transformation for three-phase system 
is available in the literature [55-56] . For six-phase system, 
Fuchs and Rosenberg [ 5o] have taken orthogonal transformation 
as the two sets of three-phase transformation. Here, the 
orthogonal transformation for the n-phase system is presented 
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and which is as follows 


T = V- 
pn n 


cos© 

sin® 

1 

V 2 


cos (®-m) 

sin(®-m) 

1 
V" 2 


cos. (©-2m) 
sin(®-2m) 
1 

V*2 


cos[©-(n-l)m] 

sin[Q-(n-l)m] 

1 

V 2 


( 2 .: 


where n is the number of phases and m = 2n/n. 


This transformation, which can be considered as a 
generalized Park' s transformation for n-phase system, trans- 
forms the phasor time varying quantities into component time- 
invariant quantities and thus simplifies the analysis. 


It can be seen that this generalized transformation 
is linear and power-invariant. This is in general a rectan- 
gul ar matrix and its Pseudo-inverse is given by 



T^ [ T x T fc ] - 1 

pn L pn pn J 


(2.26) 


where, t denotes the transpose of the matrix. 


Thus, the Pseudo-inverse of T _ is 

pn 


T 


-1 


pn 


- r 


n 


cos© 


sin© 


cos(Q-m) sin(®-m) 


1/ V2 
1/ f2 


cos[©-(n-l)m]sin[©-(n-l)m] 1/ V" 2 


(2.27) 


If the six-phase winding of the stator is constructed as the 
two sets of balanced three-phase windings displaced at e 
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electrical degree* the corresponding transformation can be 
given by 




1 cos© cos(Q- e) cos (Q-12o°) cos (©-120°*- e ) 
sin© sin(©- e ) sin(0-12o°) sin (0-12o°-e ) 

l/f 2 l/f 2 l/f 2 l/f 2 

%- . _ 

cos (®+12o°- e ) 

sin(®-12o°~ e ) 
l/f 2 


cos (©+120°) 
sin (©4-1 20 ) 

1/ f2 


( 2 . 28 ) 


By substituting corresponding e , a linear orthogonal power 
invariant transformation can be obtained for the purpose of 
simplifying the analysis. 


As the rotor quantities are not affected with this 

I 

transformation, the augmented transformation can be 

obtained by including also rotor windings and the transfor- 

i 


mat ion T_, 
P6 


thus becomes. 


I 





(2.29) 


where I is the identity matrix of order (4x4) 



41 


2.2.3 Dynamical Equations of a Six-Phase Machine 


The voltage equations for the six-phase synchro- 
nous machine in the phasor form is given by 

d 


V. 


r . a w 

] __p " (Jt -4P 


where. 


= [ V V, . . . V _ V„ V V V ] 

— p L a b rPcno-J 


f v __ v V _ V _ 

F G D Q 


ip = ^ ••• ^ ^ i G i D ^ 3 

3? = [* ^ ... ^ ^ f 0 ] 


"a b 


f F G .D Q 


(2.30) 

(2.31a) 

(2.31b) 

(2.31c) 


and [ r] is a diagonal matrix whose diagonal elements are the 
resistances of the windings, and V G = = v = 0 as these 

windings are short-circuited. 

The flux linkages of the windings in the phasor form is 

it =[L a.. Q Jip (2 - 32 > 

where [ L n lis the inductance coefficient matrix and is given 

a- • 
as 


[ L 


a. 


L aa L ab L ac ** L af L aF L aG L aD L aQ 
L ba ^b L bc ** L bF ^toG ^bQ 


(2.33) 


L L , L 
Qa Qb Qc 


L Qf L QF l qg l qd l qq 
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The current, flux and the voltage in the transformed compo- 
nents are, . _ r T ' 1 i (2 34) 

-park ~ L t P 6J ip 

-Park = ^- T p 6^ ^p (2.3 5) 

(2.36) 


V = Ct* ] v 
—Park P6 — P 


Multiplying (2.32) by [t g] and substituting the value from 
equations (2.34) and (2.35), in following expression is 
obtained. 

t 


¥ 

—Park 


“ [T p6^ 


-Park 


(2.37) 


Substituting the values of [ L _] from equations (2.11)- 

Sl* *Q 

(2.16), (2.20a-2 .20p) and (2.21- 2.24) in equation (2*37), 
we have. 


X 

.. I 


c 

*0 

c 

% 

=) 

% 

' 

' 

■c 

V 

] 

\ 

c 

i :L 


J d 


0 

L c 

0 


mF 


0 

0 

L c 

0 


V3 L 


mF 


0 


V" 3 


mD 


0 


V" 3L. 


mD 


-V\ G ° 
0 0 


-^ 3L m Q 0 


0 

0 

l ff 

0 

l fd 

0 


— V*3 L 


mG 


0 

0 


GG 


0 


-L, 


GQ 


0 

0 

l fd 

0 

0 


-V 3 L 


mQ 


0 

0 

- L 


GQ 


QQ 


q 

i 

o 


D 


Q 


( 2 . 


where , 


L d L aaO 


+ 2l , + 3L 0 

abo aa2 


L = L + 2l , -3L 0 

q aaO abo aa2 

L o — L aaO L abo 


(2.39a) 

(2.39b) 

(2.39c) 
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Further, the equation (2.30} is manipulated with the help of 
equations (2.34)-(2.36) to obtain the voltage equations in 
the transformed variables and these are given in the matrix 
form as below. 


rill* — "—I 

V 

d 


r 

a 

0 

0 

0 

0 

0 

r 

o 


’X] 




h ^ ~ 
1 

V 

’ 

0 

r 

0 

0 

0 

0 

0 


i 




! 

q 



Q 







q , 


q 


*d 

V 


0 ' 

0 

r 

0 

0 

0 

0 


i 




0 

O 

=- 



a 






o 


O 





0 

0 

0 

r_ 

0 

0 

0 




p¥ 

1 

0 

F 





F 



F 


F ■ 


V G =° 


0 

0 

0 

0 

r G 

0 

0 


H 


Pl ?G 


0 

o 

I! 


0 

0 

0 

G 

0 

rD 

0 


i i) 


P^ 


0 

l v Q=0„ 



0 

0 

0 

0 

0 

r Q ** 


lA. 


2^' 

- 

0 


where, P = d/ dt and w = dG/dt. 

Thus, equations (2.40) in matrix form gives the voltage equa- 
tions in d,q,o component for a six-phase machine 


Equation of Motion 

The dynamic equations corresponding to the electro- 
mechanical motion of the machine are 


and M -f£- = [ T m - T, - Kjj (o> - o o ) ] 

where M is inertia constant in sec/radian/sec. 
T m is the mechanical torque 
T e is the electrical torque 
Kp is the damping coefficient. 


(2.41) 

(2.42) 
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2.3 LINEARIZED MODEL OF A SIX-PHASE (MULTI-PHASE) SYNCHRO- 
NOUS MACHINE 

A detailed linear model of a six-phase (multi-phase) 
synchronous machine in the state-space form has been deve- 
loped using the currents as state-varaibles. The machine 
has been considered to consist of six windings on the stator* 
one field winding F and two fictitious windings D and Q on 
the rotor. The performance equations to represent the six- 
phase synchronous machine have been derived on the same lines 
as those in section 2.2. All the quantities are in p.u. 
except time and angle which are in second and radiant respect- 
ively. 


*d 

= Vd +f3 h + 

X mD 

(2.43) 


= V " 3X mF Wp V + X FD 


(2.44) 

1 

D 

= f3X M) ** + X FD % +X DD ^ 

(2.45) 

¥ 

q 

= x q V-^Wa 


(2.46) 

Q 

= -V" 3 
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(2.47) 

V d 

= -r 

a 

1 

ld “ ^ 

(0 

P *d - - 

o 

*q 

(2.48) 
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F 
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(2.49) 
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(2.50) 
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G> 
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*d 

(2.51) 
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Q 
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- r_ i 
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i 
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“o Q 

(2.52) 

P6 
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w o 



(2.53) 



(2.54) 


po 


T. 
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(2.55) 


The set of equations (2.43) -(2.55) have been mani- 
pulated to obtain the state space model as follows. 


where 


EX = f (X) + DV 


X [ 6 , to , ijj* ig , i-p » i D i 1 q] 


V = T 


m 


(2.56) 

(2.57a) 

(2.57b) 


and 


E = 


1 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

x d 
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^3X _ 
mF 

V~3X 

mD 

0 

0 

0 

% 

0 

X 

q 

0 

0 
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mQ 

0 

0 
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0 
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D =[°“ 0 0 0 0 0] (2.57e) 

Using Taylor's series expansion of equation (2.56) 
about an operating point and retaining the first order terms 
in the expansion, the linearized model is obtained in the form 


E 5m " F 2m + G 


(2.58) 


where. 


X m =[a 6, &« , A± d * Ai q / Ai F / Ai D # Ai J ( 2 . 5 8a) 

t 

— — r 

i 


G = 


0 ^ 0 0 0 0 0 
M 


-to r 

0 0 0 0 — o 0 


X 


mF 


(2.58b) 


Hjn C AT m AE FD^ 


Q . . 58c) 
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where 


a 

= < X d 

- X ) i 
q qo 

+ V" 3 X_ 
mQ 

i 

QO 

(2.58e) 

b 

= < X d 

- X ) i _ 
q do 

+ ^3 X _ 
mF 

± I0 + X mD *D0 

(2 , 58f ) 

c 

= X_ I 

- V3 X _ i « 


(2„58g) 


q 

qo 

rrm q.u 


d 

= x d 1 

do + V”3 X mF i FQ + 

X mD "'"DO 

(2.58^5 


The equation (2.58) reduces to 

X = [a ] X +[B 1u (2.59) 

-m L rtv -m L m-* — m 

where [ AjJ = E -1 F and[l^] = E~^G. (2.60) 


The equation (2.59) gives a detailed linearized model in state 
space form of a six-phase synchronous machine. 


2.4 A CASE STUDY 

The dynamic stability of a six-phase synchronous 
machine (Fig. 2.2) has been investigated first by taking damper 
windings D and Q into consideration and then by considering the 
machine without damper windings under two boundary conditions 
namely' (a) the rated power and 0.85 lagging p.f. with generator 
terminal voltage as 1 pu and (b) the infinite bus voltage is 
1.0 pu and the machine loading remains the same as in (a). 

The machine is connected by a transmission line of 0.02+j0. 4 pu 
impedance to an infinite bus. The initial condition of the 
machines have been determined. The parameters of the six-phase 
synchronous machine have been obtained by assuming its indue- 


0 


10.2.2 


TRANS. LINE 


Re a *02 P.U. I 

X £ a * 4 P.U* INF BUS 

RF.s0.B5 


A six-phase (Multi-phase ) machine 

connected to an Infinite bus through 
a transmission line. 
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tive coefficients equal to those of a three-phase synchronous 
machines, and the stability of six-phase and three-phase 
machines have been compared. However, the analysis of three- 
phase and six— phase machines have been carried out on pu values 
based on the single-phase quantities. The eigenvalues of 
system matrix [a^] have been determined to investigate the 
dynamic stability, and further, initial value responses of 
machine angle and angular velocity versus time have also been 
obtained. 


2,4,1 Determination of Initial Conditions of Six-Phase 
Synchronous Machine 

For the stability study, it is essential to deter- 
mine the initial conditions prior to the application of the 
disturbance. Therefore, the steady state equation for six- 
phase machine is required for calculating the operating point 
conditions. 


From equations (2.43) to (2.52), all the flux 
linkages and the currents at the steady-state are constant 

© a t 


p %. =0 
x = d,q,F,D,Q 
io = i Q = e ; « 0 

Then, (2.48) and ( 2 . 45 ) give after manipulation 


where 
and also. 


co 


(2.61) 


(2.62) 


v, = -r i , - X i 
d ad q q 


(2.63) 
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and v = — r i + X, i. + V3 X m _ jL_ 
q a q d d ' mF F 


(2.64) 


From equation (2.36), we have. 


v a 6 ^ V d GOS ® + sdn ®) 


q 


(2.65) 


From equations (2 .63) -(2 .65) using © =co t + 6 +n /2 


v a =V g[-( r a i d +Xqi q ) cos (o^t+6^/2) + (~r a i q +X d i d +f 3X mF i p ) sin(co o t+6+fc/2)] 


=^§ [- (r a i d +X g i q ) cos (G> o t+6+tt./2 ) + C-r^+x^+Vax^ip ) cos (c^t+6)] 

-( 2 . 66 ) 

At open circuit, a field current ip corresponds to an e.m. f. 
of ip X^p peak and if the r.m.s. value of this peak is E, Then 

j F X mF = ^ 2E and ^" 3 % x m F = V6 E (2.67) 

when E is the stator equivalent EMF corresponding to ip. 

Introducing equation (2.67) in (2.66) by defining r.m.s. voltage 
v 

MM. Q 

phasor V = — , the following expression is obtained 

a f2 

\ ZS + J ff Z6)-J + x d iS/6 + E/6 (2.68) 

In this equation V & and E are stator r.m.s. phase voltage in 

P.u. 

If the r.m.s. equivalent d and q-axis currents are taken as 

J a - ff (2 - 69) 

The stator current i , expressed as a phasor, will have two 

cL 

rectangular components I and I 


Taking q-axis as the re few 



Ate, 


u u 

p^,j 


hence, we have. 


I = I + j I 
a q 


d 


(2.70) 


Substituting (2.69) and (2.70) in (2.68), we get, 
E 4i = V a + r a I a + jx q 1^$ - X a I d /6 

By using E = E/6, I = 1^/6 and 1^ = j 1^/ 6 
Equation (2.71) becomes 

E = V +r I + j X I + j X , I 


(2.71) 


a a a 


q q 


d d 


( 2 . 72 ) 


The equation (2.72) represents the steady state equation of 
six-phase machine and its phasor diagram is shown in Fig. 2.3. 

To obtain v^ and v^ from (2.63) and (2.64), the 
r.m.s. stator equivalent voltages can be computed as 
v. 


V = ~~ 

d If 6 


- r x , - x i 
ad q q 


v 


and V, 


= - r I + Xd _ + E 


(2.73) 

(2.74) 


q f 6 a q d d 

Case a : When the terminal voltage and the generated power of 
six-phase machine are known; 

It is desired to obtain the initial condition of 
the six-phase machine under the above condition. 

With the above boundary conditions, the terminal 
voltage v & the current and the power factor F p are known. 
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Resolving I into components with V as the reference, 
a. cL 

= J r + j I x (2.7S) 

where I is the component of I in phase with V and I is 

xt 0. ct yz 

its quadrature component, we also have, 

F = cos 0 

P p 

where jzL is the angle by which I lags V . Then, 

P a a 

I r = I a cos 2f p ; I x = - I a sin «! p (2.76) 

The phasor E in Fig. 2.4 is given by 
qa 

E = V + (r + j X ) I = V + (I + jl ) (r + jX . ) 
qa a a q a a r J x a q 


(V - x I + r I ) + j (X I + r I ) 
^ a qx ar / Jv qr a r 


(2.77) 


The angle between the q-axis and the terminal voltage is 


given by 


-1 


6 -■ (3 = tan" [ (x q I r + V x )/(V a + r a I r - x g I x ) ] (2.78) 


Then we compute, 

v d = - v a sin (6-3) ; V g = V a sin ( 6 ^ 3 ) (2.79) 

and v^ and v can then be determined from equations (2.73) and 
(2.74) . 


Further, the currents are obtained from 

I d = - I a sin (6 -3+0) *1 = I a cos (6-3+0) (2.80) 

and then, h and ± can be determined from relations in 

^ St 

equation (2.69). 



*_q.-axis 



FIG-2-3 A phasor diagram of a Six- phase 
synchronous Machine . 



FIG-2.4 A phasor diagram of a Six- phase 
Machine for the location of q.-axis 
under the boundary condition of Case a 
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From equation (2.74) , E can be obtained. 

The field current = ]f6E/( V"3 (2.81) 

Thus, the initial conditions in terms of current variables 
can be known. 


Case b : When the infinite bus voltage and the power of the 
six— phase machine at a particular power factor are known: 

Since the boundary conditions in the case are mixed 
and therefore, the initial conditions can be obtained with the 
modification in (a) as shown below. 


The machine is connected to the infinite bus through 
a transmission line having resistance and the reactance as R e 
and X e respectively. 

The power at the infinite bus has been approximated 
on some estimate of current causing the power loss in the 
connecting lines. 

The power delivered to the infinite bus = Power output of the 

2 

machine — I r (2.82) 

3 . £ 

The angle & between I and V M is given by 

0 . 

J x 

tan® = (2.83) 

r 

From Fig. 2. 5, the angle p between and V M is given by 


XI. + R J. 
e r ex 


Voo - XI 


e x 


4. R I 

+ e r 


tan p = 


(2.84) 



d -a xis 



FIG. 2.5 A phasor diagram showing the terminal 
voltage of a six- phase machine and the 
voltage at the infinite bus- 
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The power factor angle jzf at the machine terminal is given by 

jZf = |3 + © (^ . 85} 

Using the trigonometric identify, 

tan '/ = (tan J3 + tan Q) / (1 — tan j3 tan 0} (2.86) 

From equations (2.85) and (2.86), I can be computed. 

The terminal voltage V is given by 

cl 


V 

a 


(Voo * X I + R I ) 
ex e it 


+ J <Vr 


+ Vx 5 


(2.87) 


and 

^ 6 =[ «q + X e> 1 r + <VV 1 x] 7 f V “ -< X q +X e> VWU 


( 2 . 88 ) 


Then the initial conditions of the synchronous machine in 
currents variables can be determined as given in case a. 


2.4.2, Dynamic Stability of a Six-Phase Machine with Damper 
Windings . 

A 60 Hz six-phase synchronous machine has been 
examined under two boundary conditions for which the methods 
for calculating the initial conditions has been discussed in 
the previous section. The machine is connected by a trans- 
mission line of (0.02+j0.4) pu impedance to an infinite bus. 
The parameters of the synchronous machine have been provided 
in Table 2.1. The six-phase machine parameters have been 
obtained by using those of the three-phase machine. Pu para- 
meters have been provided on per phase basis. The parameters 
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and have been obtained from three-phase machine using 
equations (2 . 39a) - (2 . 39c) by assuming equal inductive co- 
efficients (L _,*L 0 , L , } of the two machines. Data if 

aaO aa2' abo 

any, not mentioned are referred to reference [ 55 ] . 


TABLE 2.1 

Synchronous Machine Parameters 


x d 

t> 

« 

II 

r 

a 

= .001096 


X 

g 

= 1.64 

■'f 

= .000742 


* 

x d 

= 3.3 

r D 

= .0131 


* 

X 

q 

= 3.18 

r Q 

= .054 


X FF 

= 1.65 

V 1 
2 

X mQ = X ' 49 


X DD 

= 1.605 

V 2 
2 

X mF = X FD 

- X 

2 mD 

X QQ 

= 1.526 

M = 

= .037719; 

= 0.0 


* indicates the quantities which give the corresponding pa- 
rameters for a six-phase synchronous machine 


The initial conditions have been obtained as 
discussed in section 2.4.1 with the help of machine parame- 
ters, and have been provided in Table 2.2. 

The six— phase machine has been considered under 
two conditions : 1) when its rated power is equal to the 
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rated power of a three-phase machine and 2) when its rated 
power is equal to twice that of the three-phase machine. 


TABLE 2 . 2 


Operating Conditions of Three-Phase and Six— Phase 
Machines Under Two Boundary Conditions 


Boun- 

dary 

condi 

tion 


3 -Phase 
Machine 
rated 
Power 

6 -Phase 

Machine 

quantity 

rated Power 
equal to 3- 
Phase Mach- 
ine 

rated Power 
equal to the 
twice of 3-Phase 
Machine 

a 

6 (in deg 
° ree) 

66,994 

51.289 

74.8296 


1 do 

-1.9253 

-1.3579 

-2.8261 


i q° 

0.6673 

0.4818 

0.56369 


i 

FO 

2.9794 

2.0626 

3.5482 

b 

6 o 

53.7394 

45.87 

62.48 


"'‘do 

-1.584 

-1,132 

-2.3695 


1 qo 

0.6984 

0.4676 

0.6205 


1 F0 

2.8194 

1.914 

3.1941 


The dynamic stability has been investigated of 
the three-phase and six-phase machine under different con- 
ditions assuming both the machines with damper windings. 
This has been investigated by obtaining the eigenvalues of 
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the system matrix [ A m ] corresponding to the different opera- 
ting conditions of three-phase and six-phase machines using 

program EIGRF of the IMSL subroutines. In this study, it is 

t 

assumed that the input vector U i,e. [a T- a E__l is 

~ m x xn FD 

zero. The eigenvalues of the system matrix[ A^^are provided 
in Table 2.3. 


TABLE 2 . 3 


Eigenvalues of the System Matrix [ AJ for Three-Phase 
and Six-Phase Machines with Damper windings for Two 




Boundary 

Conditions 

Boun- 

dary 

Condi. 

tion 

Eigen- 

values 

3 -phase synchro- 
nous machine 
(rated power) 

6 -phase synchronous machine 
equal power to double power to 
3-phase mach- 3 — phase machine 
ine 

a 

X 1' X 2 

— 13.54+j 376. 335 

-11. 738+j 376.0 -11 .74+] 376.0 


?\. 3/ x 4 

-0. 3558+j 9. 345 

-0.8177+5 12 .076 -0 .5857+j 12 .4157 


^5 '^6 

-37.71# -46.37 

-39.81 , -70. 6854 -39 .6974 , -71 . 3896 


x 7 

-0.20234 

-0.3368 -0.212527 

b 


-13.5513+5 376.335 -11.723+5376,0 -11. 742 +j 376.0 


A3/A4 

— 0.5864+j 10.8792 

-1.0675+512.5589 -0 .81 33+j 14 . 77 


^5 # ^6 

-37.6865,-45.882 

-39.8104,-70.194 -39.594,-70.96 


X? 

-0.2611 

-0. 3654 -0.286874 


The machine angle and the angul ar velocity versus 


time responses have been obtained using IMSL subroutine DVERK 
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by applying a disturbance of a 6 = .1 radian and these res- 
ponses for case a and case b are provided from Figs. 2.6- 

2.9. 

2.4.3 The Dynamic Stability of a Six-Phase Machine Without 
the Damper Windings 

For the dynamic stability investigation of three- 
phase and six-phase machines without damper windings/ the 
section 2.4.2 is repeated beeping in view that/ when the 
linearized model in equation (2.59) is considered for the 
machines without damper windings/ the rows and column corres- 
ponding to the damper winding quantities (6th and 7th rows 
and columns) are removed reducing the order of system matrix 
[A^'J to (5x5.)/ and correspondingly, matrix[B in ] reduces to the 
order of (5x2). The initial conditions obtained in Table 2.2 
can well be applied in this case for the different conditions 
of three-phase and six-phase machines, as the damper winding 
parameters do not affect the initial conditions of the 
machines. 

The eigenvalues of the system matrix [a^] corres- 
ponding to three-phase and six— phase machines under two 
boundary condition have been provided in Table 2.4. 

The machine angle and angul ar velocity response 
for three-phase and six-phase machines without damper wind— 
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ings under two boundary conditions with zero input vector 
have been shown in Figs. 2.10-2.13. 

TABLE 2.4 


Eigenvalues of the System Matrix [ A^] for Three-Phase 
and Six-Phase Machines Without Damper Windings for 
Two Boundary Conditions 


Boun- 


3-phase synchro- 

6 -phase synchronous machine 

dary 

Condi- 

tion 

Eigen- 

values 

nous machine 
(rated power) 

equal power to 
the 3 -phase 
machine 

double power 
to the 3 -phase 
machine 

a 

* ^*2 

-8.129+j 376.9612 

-6.162 3+j 376.96 

-6.1674+j376.9624 


X 3'^4 

-0.165+j 9.254 

—0.2173+j 11 .8757 

-0*277+312.22 


^5 

-0.2124 

-0. 352548 

-0.2229 

b 

\i # \ 2 

-8.129+j 376.96 

-6.154+j376.957 

-6.1658+j 376.962 


\ 3 i \4 

- 0 . 1 35 + 3 10 . 722 3 

-0.2l06+jl2.32 

-0. 239936+jl4.499 


^5 

-0.2738 

-0.382365 

-0.300374 


2.4.4 Results and Discussions 

The operating conditions of three-phase and six- 
phase synchronous machines under two boundary conditions , 
namely under case a and case b have been given in Table 2.2. 
Eigenvalues of system matrix [ A^ ] corresponding to three- 
phase and six-phase synchronous machines with damper wind- 
ings under two boundary conditions have been provided in 
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Table 2.3, and those corresponding to three-phase and six- 
phase machines without damper windings in Table 2.4. It is 
observed from Table 2.3 that the pair of eigenvalues 
and the eigenvalue of the system matrix [ ] almost deter- 
mine the stability of the system. The pair in column 

3 for three-phase machine corresponds to frequencies of 
approximately (9.345/2;j) 1.48 Hz and the oscillations are 
damped with a time constant of (1/.3558) 2.81 Sec. for case a 
and the corresponding frequency and time constant are 1.73 Hz 
and 1.7 Sec for case b. When a six-phase machine has been 
considered to deliver power equal bo that of three-phase 
machine, the corresponding values for case a and case b are 
obtained aS 1.92 Hz, 1.2 Sec and 1.99 Hz, .936 Sec. respec- 
tively. 

When the six-phase machine has been considered to 
deliver power twice that of three-phase machine, these values 
are 1.97 Hz, 1.707 Sec. and 2.35 Hz, 1.23 Sec. respecti- 

vely. 

From Table 2.4, when the machines without damper 
windings are considered, it is observed that the pair 

and are the eigenvalues which almost decide the 
dynamic stability of the system. For a three-phase machine, 
the pair corresponds to the frequencies of approxi- 

mately 1.47 Hz with time constant 6.06 Sec for case a, and 
1.7 IIz, 7.4 Sec for case b. When the six-phase machine 
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delivers power equal to that of three-phase machine, the 
corresponding values are 1.89 Hz, 4.6 Sec. and 1.96 Hz and 
4.6 Sec. for case a and case b respectively, while for the 
six— phase machine delivering power twice that of three-phase 
machine, corresponding values are 1.94 Hz, 3.61 Sec and 2.3 Hz, 
4.16 Sec. 

It is observed that X^ in Table 2,3 and X^ in 
Table 2.4 have shifted left from the origin more in case of 
six-phase synchronous machine while it delivers power equal 
to that of three-phase machine and even, when a six-phase 
machine delivers double the power, the corresponding eigen- 
values for the different cases of six-phase machine have 
better location than those of three-phase machine. 

The initial value responses corresponding to 
A6^'t and < — t have been plotted and shown in Pigs. 

2.6-2.13 for the three-phase and six -phase machines deliver- 
ing equal power for different cases with and without damper 
windings. It is observed from Figs. 2.6 and 2.8 that the 
oscillations in a 6 have damped out rapidly than those 
corresponding to three-phase machine for case a and case b. 
Prom Figs. 2.1o and 2.12, when the machines are considered 
without darrper winding, the oscillations are very slowly 
dying out; however, in the case of six-phase machine, they 
die a little faster and with lower peak. From Figs. 2.7 and 
2.9, the plot of aW t for three-phase and six-phase 
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FIG. 2.6 INITIAL VALUE RESPONSE OF A6~t FOR 3-PHASE AND 
6-PHASE MACHINES WITH DAMPER WINDINGS FOR Case 
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machines with damper windings for case a and case b, it is 
observed that oscill ations die out faster in the case of 
six-phase machine, but with a higher peak for 1 cycle; and 
this is evident from the pair of eigenvalues in 

Table 2.3, Further, it is observed from the Figs. 2,11 and 
2.13, the oscillations tend to die out faster in case of 
six-phase machines but initially with higher peak for three- 
cycles, 

2.5 CONCLUSION 

The various inductances of a six-phase (multi- 
phase) synchronous machine have been obtained from the funda- 
mentals. An orthogonal (linear and power invariant) trans- 
formation for a multi-phase system has been proposed to trans- 
form the phasor quantities into dqo component quantities in 
order to simplify the analysis of multi-phase synchronous 
machine. A linearized model has been developed for a six- 
phase machine to carry out the dynamic stability analysis of 
the six-phase synchronous machine. It is observed from the 
dynamic stability investigation that the six-phase machine 
while delivering power equal to that of a three-phase machine, 
is dynamically more stable than, tne three-phase machine, and 
also when it delivers double power, the stability of six-phase 
machine has improved as compared to that of three-phase 
machine. 



CHAPTER 3 


GENERALISED MATHEMATICAL MODELLING 

OP n-phase POWER SYSTEM ELEMENTS 

3.1 INTRODUCTION 

In the previous chapter, a detailed dynamic model 
of a six-phase (multi -phase) synchronous machine has been 
developed and the dynamic stability of a six-phase (multi- 
phase) synchronous machinelias also been investigated. In 
the current chapter, mathematical descriptions (models) of 
various n-phase (multi-phase) elements of power systems have 
been obtained, for the mathematical modelling is a very 
flexible and inexpensive tool for the in-depth and the 
detailed analysis of the behaviour of systems. Usually the 
analysis of three-phase system is conducted on the single- 
phase basis assuming balanced conditions; the power system 
is normally not balanced due to the various causes, viz., 
long un transposed transmission line, bundled conductors, 
single-phase switching, unbalanced loading, etc. Therefore, 
in order to take these factors into consideration, it is 
necessary to develop the mathematical model in the phase frame 
of reference. 

A mathematical description of a n-phase (multi- 
phase) synchronous machine has been developed in section 3.2 
incorporating imbalances of the machine. If the power is' 
transmitted through a n-phase transmission line, it will 
require a three-phase/n-phase transformer at each end of the 
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line requiring thereby a mathematical model for such trans- 
formers. 

In section 3.3* three-phase/n-phase transformer 
models have been developed for different types of winding 
connections. This is done by first obtaining the voltage and 
current relations of the two sides of the windings «f three- 
phase/n-phase transformers in section 3.3.1, and using these 
relations, the nodal representations of the transformers 
have been obtained in section 3.3.2. 

Further, mathematical models for multi-phase trans- 
mission line have been given in section 3.4 in terms of phase 
impedance matrix, ja-circuit representation and ABCD parameters. 
If the three-phase transmission line is unbalanced and the 
interest of investigation lies on three-phase side of the line, 
it is expedient to replace the balanced n-phase line associ- 
ated with the transformers by their equivalent three-phase 
parameters to carryout the analysis. The equivalent three- 
phase representation of such n-phase (multi -phase) line has 
been derived in section 3.5. The equivalent three-phase 
impedance has been derived in section 3.5.1 and the equivalent 
three-phase ABCD parameter representing r -circuit in section 
3.5.2. Further, for the composite system to be corrpletely 
balanced, only single-phase parameters are required. The 
equivalent single-phase n circuit parameters have been derived 
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from the equivalent three-phase ABCD parameters in section 

3.6'. 

When a three-phase transmission line is connected 
between the two n-phase (multi-phase) transmission lines 
through n-phase/thre e-phase transformers, it is necessary to 
•btain the equivalent n-phase representation of the three- 
phase transmission line associated with two transformers at 
both the ends. These parameters are desirable when the 
analysis of the composite system is carried out on n-phase 
basis. The equivalent n-phase parameters in impedance form 
have been derived in section 3.7.1 and those in ABCD para- 
meters representing n circuit in section 3.7.2. 

Finally, the chapter concludes in section 3.8. 

3.2 MULTI-PHASE SYNCHRONOUS MACHINE MODEL 

This section presents a mathematical description 
of a multi-phase synchronous machine in phasor co-ordinates 
for steady-state analysis under balanced as well as under 
unbalanced conditions. In order to develop the model, a 
schematic representation of a star connected multi-phase 
(n-phase) synchronous machine, as shown in Fig. 3.1, has 
been considered. The generated (internal) voltages of the 
machine are assumed to be balanced, and the terminal voltages 
unbalanced as the network connected to its terminal (i.e. 
load) are unbalanced. 
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Let E^, E 2 , . .., be a balanced set of internal 
voltages of a n-phase synchronous machine, and nodes n+1, 
n+2,...,2n be joined together bo form a star (i.e., neutral) 
point; a n-phase synchronous machine can thus be represented 
by a balanced set of internal voltages behind synchronous phase 
impedances. 


where 


ind 


Then, from Fig. 

3.1, 

the voltage 

form are. 
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(3.1) 


n 


= (nxn) impedance i. e. . coefficient matrix of the n-phase 


synchronous machine (= [ Y? ] ) 


For the voltage and current directions as shown 
in Fig. 3.1, the nodal, representation of the machine can be 
written as. 
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where 


I P £ J 1 I 2 ••• x n ] 


and 

*£ = [ Zp] 


= Admittance matrix of the n -phase synchronous 
machine of order (nxn) w 


Row, for a balanced set of internal voltages, 
,n 
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(3.3) 


where 


= e^ 2n//n - cos (2jt/n) +j sin(2rt/n) = l/2rt/n = nV" 1 


.n 


The voltage vector E p in eqn. (3.3) can also be written as 

(£-1)* * o* (§-l)* t 
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(3.4b) 


when n is an odd number 
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If the nodes n+1 to 2n be joined together to form the star 
point, then we have. 


and 


I N ~ I n+1 + X n+2 + + I 2n 


V N = V n-i-l * V n+2 = ~ * V 2n 


(3.5) 


(3.6) 


The internal voltages can then be given by 
E P =[e 1 +V n E 2 +V n ••• E n + V 


(3.7) 


n 


and the currents I , by 


' la wT ' l2a (e 2 + v 


‘ I na = 
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(E +V ) 
n N 




(3.8) 


where S 2 , . ..,S n are complex powers for n phases (1 to n) 

respectively, as these are unequal. 

Further the (nxn) admittance matrix in equation (3.2) is 
given by 
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If the data are available in symmetrical component form, its 
value in the phasor form can be obtained by the relation 
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(3.10) 
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y 0 #y 1 # •••/ Y n ^ are n-phase sequence admittances, 
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[_ T^J is the generalized symmetrical component transfor- 


mation matrix for n-phase system which can however be developed 
by exploiting the cyclic symmetries inherent in the network 
using the group theoretic technique, and it can be -written as 
follows. 
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where the operator <x = nY" 1 = l / 2 rc/n = e^' ^ rZ// ' n 
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Using equations (3*10) through (3.11), the various elements 
of the coefficient matrix in equation (3.9) can be obtained 
which are as follows. 


y s = I c y c + , s , yj 

x=l 

= I [y c + , s , a2i y±l 

i=l 

i n-1 o . 

y m 3 = s' + . ^ a y± 3 


(3.12a) . 
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Substituting the values of Ep, Ip, and from equations 
(3.7)— (3.9) respectively in equation (3.2), the following 
representation of the n-phase synchronous machine in matrix 


form is obtained, 
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(3.13) 


By adding the last n equations of (3.13), an additional 
equation for node N (neutral node i.e. star point) is obtained 
as 

I N =[-y 0 -y Q ... ny o O] [v» V N E x ] ( 3 . 14 ) 

From equation (3.3), - 

= [ 1 cc a 2 ... a^ n “ 15 ] E x (3.15) 

cross multiplication of the last n equations of (3.13), and 
their addition and substitution of different element of 
admittance matrix in equation (3.2) gives the following ex- 
pression after simplification. 
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In equation (3.16), (n+1) terms are V times I and the 
Injected neutral current I = 0. Taking this point into 
account, the equation (3.16) can further be simplified as 


r 1 ,y -Y 2 rr n “ 3 ,n-l -I 

L —l — *cc — *cc • * • -a. — (X *~cc n J 


! g 


(3.17) 


combining the first n equations of (3.13) to (3.17) and taking 
I N = 0, the machine equations in the matrix form become, 
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Thus equation (3.19) represents a generalized mathematical 
model for multi-phase synchronous machine. It calculates the 
terminal phase voltages, either balanced or unbalanced of 
synchronous machine. It makes use of the total power of the 
machine and not the individual power of each phase. Therefore, 
in order to determine the left hand vector in equation (3.19), 
individual phase power will not be required. In the right hand 
side, the current vector 1^ is the injected current, so, if 
there is no local load on the synchronous machine bus, the 



injected currents are zero, if, however/ there are local 
loads on the machine, injected current cannot be taken as 
zero, rather one nth of the total local load can be assigned 
to each of the phases. Due to the imbalances in the terminal 
voltages because of unbalanced terminal currents , each phase 
power is not equal to one nth of the total power thereby 
precluding the calculation of phase currents, 

3.3 THREE-PHASE/MULTI-PHASE (n-PHASE) TRANSFORMER MODEL 

The power system networks will essentially require 
three-phase/n-phase transformers for multi-phase transmission 
line in order to be compatible with the three-phase synchro- 
nous generator and the three-phase loads. Therefore, mathe- 
matical models for such transformers will be needed to analyse 
the over all system. Therefore, this section is devoted to the 
nodal representation of three-phase/n-phase transformers by 
taking Wye or delta windings for three-phase side and star 
connection for n-phase side of the transformer by presenting 
the voltage and current relations of the two sides windings 
of a transformer with the tappings on both the sides, and 
therefrom various particular cases can be derived. 

3.3.1 Voltages and Current Relations of Three-Phase/n-Phase ! 

Transformers 

The voltage and the current relations for primary 
and secondary sides of the three-phase and six-phase transfor- ■ 
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mers are available in the literature [8,22,23 ] . Although 
a three-phase/ six— phase transformer is sirrple in construction, 
yet it is not so simple to construct a three-phase/twelve- 
phase transformers. If, however, it is realized, it is 
required to be especially built* Moreover, the constructional 
complexity will be more involved for a three-phase/n-phase 
transformer, where n is larger than twelve. In order to 
develop the voltage and current relations for the two sides 
of three-phase/n-phase transformers, an equivalent n-phase 
transformation has been obtained from a given three-phase 
supply based upon the following principle. 

' The desired n-phase system can be obtained with 
the help of an appropriate number of three-phase/six-phase 
transformers. The primary voltages are obtained for the 
different transformers from the given three-phase supply by 
getting it (three-phase supply) shifted through proper angles 
by phase shifting circuits 1 . 

The number of required transformers will be n/6 
and the angle of shift for the different primaries will be 
0,m, ..., - m) where m = 2rr/n. 

Based upon this principle, three-phase (Wye) /n— phase 
(star) transformation has been obtained as shown in Fig. 3.2. 
For this case, the induced voltages at the secondary side can 

be represented as 




*S.C.-».PHA»E SMI F T IN© CIRCUIT 


FIG. 3.2 Schematic diagram for 3-phas*/ n-phas? 
transformer 
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Equation (3.20) can be written in the conpact form as 

v£, « NVp (3.21) 

t 

where N is the coefficient matrix of [ V c ] in equa- 

tion (3.20). 

The different voltages of the secondary side and the given 
supply voltage are shown in Fig. 3.3. 

Thus, the voltage and the current relations at 
the terminals of an ideal transformer are 



( 3 . 22 ) 



/ / 



FIG. 3.3 Voltage of the given supply and the 

secondary sides at the different nodes 
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and 

Ip = Cn ]lp ; Ip = [ N ] ^ i£ ( 3 . 23 ) 

inverse 

Since [n] is a rectangular matrix of order (nx3), its/ cannot 
be determined. Therefore, by taking its Pseudo-inverse [60], 
equations (3.22) and (3.23) become 



=[n]V^ ; 

t [N] 


(3.24) 

and 

Ip = [ N ] Ip ; Ip = - 

t CN] 

*4 

(3.25) 

where. 


n 

6 




KL = 2 s /-2e. ; m. = 
i=l 1 A 




and is the angle difference between the primary of ith 
three-phase/six-phase transformer and the given three-phase 
supply. 


For a particular case of three-phase Wye/six-phase star 
connected transformer, we have 


[N] 


10 0 
0 0 1 
0-10 


-10 0 
0 0-1 
0 10 


and = 2 (3.26) 


If the off-nominal turn ratio" a in the primary 
side of the transformer is also taken into account as shown 


in the fig. 3.4, the equations (3.24) and (3.25) become 
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v p I v p = I^U]* t 3 - 2 ?) 

and 

~ a 3t [ N ] ^-p ' - l -p = ~cc k“ £ N ] Ip (3.28) 

3t m 


If, however, the tappings are provided on both sides of the 
transformer i.e. three-phase and n-phase sides, the voltage 
and current relations corresponding to equations (3,27) and 

(3.2 8) will be 

-tr 1 = JT~ W V p > V p = -jr^- -4- • [N] fc v£, (3.29) 
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(3.30) 


Further, a three-phase delta/n-phase star connected trans- 
former of turn ratio a :1 as shown in Fig. 3.5 is considered. 

3t 

Assuming nominal turn ratio, the following voltage and 
current relations on three-phase and n-phase sides are as 
follows. 
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(3.34) 
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3.4 A 3-phase wye / n-phase star transformer 
with off nominal turn ratio 0(3-^ : oC^t . 
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FIG. 3.5 A 3- phase Delta/n-phase star transfor 
with off- nominal turn ratio c* 3t : 
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where connecting matrix[P] = 



-1 


1 


0 


-1 



0 


1 


(3.35) 


From equations (3.31)— (3.34) , the voltage and current rela- 
tions are 


V?, 

P’ 


y.i[N][P] = 1 ^- [p]* [n]* 


(3.36) 


ip = Va [ [ p ] ip ; ip = ~ 


V 3 


Km 


[p]* [»f i 


(3.37) 


If the off-nominal turn ratio a^, in the three-phase 
primary side is also ta)cen into consideration, the equation 
( 3 . 31 ) — ( 3 . 34) become 


v n 
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■ \ 1 n]u p 


(3.38) 
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■ rf [ p] v > 


(3.39) 

4 
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The corresponding voltage and currents on the two sides are 
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3a, 


v£, * r^{N][P] [P] *[ N] fc vJ, (3.42) 


3<X 


3t 


m 


d - V3 a 3t [ N][p] I 3 i 4 = — j— [p]TNf l“ (3.43) 

p 3 ' Vm 



95 


Now, if the tappings on both the sides of the transformer 
are considered, the equations (3.38) -(3. 4i) become 



The voltage ani current relationship on the three- 
phase and n— phase sides of the transformer can be cast into a 
general format from the pairs of equations (3.29) , (3. 30) and 


48), 

(3.49) 

as 
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(3.50) 

(3.51) 
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where, 

t 

a 3t = a 3t ' [ H ] = [ N ] / for a three-phase Wye/n-phase 

star 

and 

a 3 t = V"3 f[ H] =[N][p] for a three-phase delta/ 

n-phase star ; and cl^ represents the tapping on n-phase side 
of the transformer. Thus equations (3*50) and (3.51) represent 
in general the voltage and current relations of the two sides 
of the three-phase/n-phase transformer with the off-nominal 
turn ratios an< ^ a n -t' w '* lere fi rst subscripts 3 and n 

denote the phase order. From these equations, such relations 
can be obtained for different cases under consideration. 

The pairs of equations (3.24)-(3.25) , (3.27)-(3.28) , (3.29)- 
(3.30), (3.36)-(3.37) , (3.42)-(3.43) , (3.48) -(3.49) are the 
particular cases with the different tappings and types of 
winding connections of the two sides of the transformers. For 
the different values of n, for exarrple n = 6,12 , the relations 
can be obtained by substituting = 2, 3-j V 3 in equations 
(3*50) and (3.51). 

3.3.2 Nodal Representation of the Transformer 

For the purpose of mathematical modelling, the 
transformer has been assumed to be represented as an ideal 
one in series with the equivalent leakage impedance/admittance 
of the winding* The terminal voltages of the transformer in 
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terms of the induced voltage v£, can be expressed in genera- 
lized form from figures 3.4 and 3.5 as 

^ = 'p' - f z £] 4 (3 - 525 


where 

[ Z^] is a diagonal matrix of order (nxn) having impedances of 
the windings as the diagonal elements = [ 


Prom equations (3.50) and (3.52), 

C ^ C H ] v | - [<] ^ (3.53) 

3fc 

and from equations (3.51) and (3.53) 
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From equations (3.53) and (3,34), the nodal representation of 
the transformer is obtained as 
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(3.55) 

Equation (3.55) describes the nodal representation of a three - 
phase/n— phase transformer in generalized form rom which a 
particular case can be conveniently obtained by substituting 
the corresponding values. 
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3.4 MULTI-PHASE TRANSMISSION LINE MO EEL 


A three-phase transmission line is represented by 
phase impedance or admittance matrix* equivalent n -circuit 
representation and ABCD parameters [ 61 ] . The multi-phase 
transmission line in general can similary be described on the 
basis of three-phase line in terms of phase impedance matrix as* 
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nn 


(3.56) 


where n is the number of phases, 


If it is assumed that the multi-phase transmission line is 
completely transposed* (though it is very difficult to realize 
In practice) the impedance matrix [ Zp ] can be made diagonal 
matrix with the help of symmetrical and Clarke's conponent 
transformations, which can be developed on the basis of higher 
order transformation available in the literature [13-15*54] , 
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Whereas the short lines can be represented by 
equation ( 3 , 56 ) t the representation of the medium multi- 
phase lines can be obtained by a straight forward extension 
of an equivalent tt -- circuit model of 3-phase line [61] and 
it is as follows. 



(3.57) 


where the second subscripts S and R denote the sending and 


the receiving ends of a multi -phase line respectively, 

[ £]- [ z pJ 


(3.58) 
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In the present formulation, every phase (1-n) has been 
assigned a node number and (n+l)th node refers to ground. 

The shunt admittances for multi-phase line is 
shown in Fig. 3.6. 

Another useful representation of the transmi- 
ssion line in terms of ABCD parameters may be expressed as. 



where A n , B n , C n and D n are the square matrices of the order 
(nxn) . 


Forti -circuit representation, 

A n = D n = diagonal [l + z£ Y^/2 . . l+ZpYgh/ 2 ] 

B = diag. [ Zp # . Zp ] 

and 

C n = diag. [ Yg h (1+2^ Yg h /4)..Ys h (l+ZpYs h / 4 ^ 

For a long transmission line, 

A n = D n = co shy 1 
B n = [ Z ] sinhyl 
C n =[ Zq]” 1 sinhyl 


(3.62) 

(3.63) 

(3.64) 

(3.65) 

(3.66) 

(3.67) 
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FIG. 3. 6 Equivalent capacitive susceptances at 
the terminating bus bars. 
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where , 

[Z ] * the characteristic impedance matrix =Vz/Y and y , pro- 
o 

pagation constant =TyZ ; Z and y are the series impedance 
and the shunt admittance per unit length of line respecti- 
vely. - 

Thus, equations ( 3 . 56 ) , ( 3 . 57 ) and ( 3 . 61 ) represent 
the complete mathematical description of multi-phase trans- 
mission line in terms of phase impedance matrix, n -circuit re- 
presentation and ABCD parameters respectively. 

3.5 EQUIVALENT THREE-PHASE REPRESENTATION OF n-PHASE TRANS- 
LINE 

In a composite system if the interest of investi- 
gation lies on the three-phase side of the system, it is 
sufficient to represent the n-phase transmission line along 
with a three-phase/n-phase transformer at the each end of the 
line by its equivalent three-phase parameters to carry out the 
analysis of the composite system. Equivalent three-phase 
representation of the n-phase transmission line has been 
derived in terms of phase impedances and also in terms of 
ABCD parameters under the following sections. 

3.5.1 Equivalent Three-Phase Impedances of a Multi-Phase 
(n-Phase) Line 

In order to derive the equivalent three-phase 
impedance, the multi-phase line has been assumed to be 
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connected between three-phase buses S and R through three- 
phase/n-phase transformers T ± and T 2 as shown in Fig. 3.7. 

Two ideal transformers and T 2 uiare connected between S and 

* I 

S t and between R and R with their leakage impedances 
connected in series between S' and S" , and between R« and R" 
respectively. Let the n-phase transmission line connected 
between n-phase buses S" and R" be represented by 

V S" - V R" = Z P ^p" 0.68) 


Now the three-phase voltage at buses S and R can 
be obtained in terms of n-phase voltages and currents V^,, , 
V», and Ip M respectively by introducing the voltage relation 
of equations (3,50) in equation (3.52) corresponding to the 
two transformers connected at each end of transmission line, 
and which are as follows. 

* i 


3 

v s = 

3tl 

1 

[H] t 

«n „ 3tl 


a ntl 

*m 

S ” a ntl 

and 

3 

v° = 
R 

a 3t2 

ct nt2 

1 

K 

[Hf 

1 

Ji a 3t2 

a nt2 



(3.6?) 


(3.70) 


where 

[ and [z£ 2 ] are the diagonal matrices of order (ran) 

having the leakage impedances of the transformers and T 2 
respectively as their diagonal elements. 





connected to 3-phase buses 5 an 
transformers. 
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Subtracting equation (3.70) from (3.69) and using the voltage 
relation of multi-phase transmission line in equation (3.68), 
the following equation is obtained. 


a ntl 

C 3 tl 

3 

V - 

S 

nt2 

__ — 

3t2 

where Ip , , = 

-11 

x s 11 

: j n 

~R" 


Vi = 


J^H] t [;[z"]-l[z£j +[z£ 2 ]]i£, ( (3. 


If the transformers T^_ and T 2 have nominal turn- 
ratios, equation (3.71) yields the three-phase equivalent 

impedance as follows 


=X [H]t[Z P +Z “ +Z?2][H] C3 ‘ 72) 

If the leakage impedances of the transformers are neglected, 
the equation (3.72) reduces to 

Z P,3q [H ^ [Z P ] [H] (3 ’ 73) 

Equation (3.72) represents the three-phase equivalent of n- 
phase transmission line associated with two transformers at 
the ends with nominal turn-ratios. 

If the n-phase transmission line is assumed to be 
completely transposed i.e. Zp is an impedance matrix of 
order (nxn) whose diagonal terms are equal to Z g and the off. 
diagonal terms are equal to Z m , then from equation (3.72), 


71) 
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J P,eq 


l * 
-=~ (Z K 
v S m 


0 


Z m V 


0 




o 


0 


~“(Z K -Z K ) 
K S m mm 


(3.74) 


whore, z s “ z s + z l + z 2 
and 

= 2 C 1 + 1/Wi + 1 / —2m + l/-(Tl/3-m) ] 


(3.75) 


Thus, the equivalent three-phase sequence impedances are 

(3.76) 


7 (o) _ z (l) _ „(2) . 

P,oq P,eq " P,eq 


1 (Z S K m - Z mV 


Particular cases i 


For a six-phase transmission line associated with the trans- 
formers at each end having nominal turn ratios, the equation 
(3.76) reduces to 


(o) _ z(^) _ Z (^) as Z — Z 

P,eq ~ P,eq P,eq S m 


(3.77) 


For a twelve— phase transmission line with the transformers 
having nominal turn-ratios, the sequence impedances are 
(o) _ (2) ^ 


where. 


■ J P,eq 

I 

K 

m 

Km 


Z (D 

P,eq 


z (2> = z* - 

P,eq S 




Z 


m 


|[ (3+2V3) + j V3J 


(3.78) 

(3.79) 


For a six— phase transmission line integrated with a three- 
phase Wye/six— phase star connected transformers, equations 
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(3.72) and (3.73) reduce' to 


and 


J P/eq 


1 

2 




+ z 


'tl 


+ z 


t2 


]N 


J P # eq 


1 , T t 6 v _ 

2 N Zp N 


(3.80) 


(3.81) 


For a twelve-phase transmission line integrated with 
a three-phase wye/ twelve- phase star connected transformer at 
its each end, the equivalent three-phase impedance can be 
obtained by substituting = (3 - j y3) in equations (3.72) 
and (3.73). 

For a three-phase delta/six-phase star connected 
system the equation (3.72) and (3,73) reduce to 

Z P,oq - I [*]*[ «fl 4 * Z tl + (3.82) 

and 

3 

Z P,eq = §[ pf [ [ z£j[N][p] (3.83) 


It is observed from the equation (3.71) that, if 

a ntl //QL 3tl ^ a nt2 //a 3t2 / the t * ire ®-P hase equivalent parameters 
cannot be determined* However, if the two ratios are equal 
(the nominal turn-ratios of the transformers is a particular 
case), the equivalent three-phase Parameters (phasor as well 
as sequence) can be determined. Moreover, it is observed 
the connecting transformers increase the magnitude of para- 
meters and from equation (3. 76), it is further observed that 
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that the equivalent sequence impedances i. e, zero# positive 
and negative components are equal. 


3.5.2 Equivalent Three-Phase ABCD Parameters of a Multi- 
Phase (n-Phase) Line 


The n-phase transmission line forming a part of 
the three-phase power system connected at buses S" and R" 
via three-phase wye or Delta/n-phase star connected trans- 
formers T^ and T 2 may be modelled by an equivalent ABCD 

parameters. 


By taking the tappings on both the sides of the 
transformers T^ are T 2 into considerations and using the 
equations (3.50) and (3.51), the voltages and currents at the 
sending ends (S,S* and S") and the receiving ends (R,R , # R") 
can be related with the help of equations (3.52) and (3.68). 
From Fig. 3.7, the equation (5.52) can be expressed in the 
following form for sending ends, 

= v£„ + [Z^] I§„ (3.84) 

From equations (3.50) and (3.84) 


a 


v; 


3tl 


a„ 


CL 


htl 




tl J 

(3.85) 


and from the equation (3.51), we get 


a 


L S 


a 


htl 

r 


3tl 




[H'f I" 


s «. 


(3.86) 
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Similary, the voltage and current relations at the receiving 
ends R # R‘ and R'* can be expressed as 


V" 

R 


a 


ht2 


3t2 


[ H ] 


a 


v. 


3t2 


R 


a, 


nt2 


[hJEz^ 


R 


and 


(3.87) 


T n _ 3t2 _3 

I R" " a ' Lll J X r (3.88) 

ntz 

From the pairs of equations (3.85), (3.86) and (3.87), (3.88), 
the following representations in the matrix form are obtained. 


1 

I 


mwmm. | 

a 3tl 

-P 

i 1 

• 

CC 3tl 

• z ti ] 


v 11 " 

S r ' 


^ntl 

a ntl 


3 

x s 



0 

a ntl 

a 3tl 

X [H]t 


I n 
X S" 
— — 


I (3.89) 

and 


1 . 

l 


~ CL 

5?^ Eh] 

a 3t2 

rftH][ z y 

nt2 

: 

OJ 

J 

n 

^R" 


0 

a 3t2 

a nt2 


I 3 

R 


(3.90) 


f 
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From equations (3.61), (3.89)* (3,90), the equivalent three- 
phase representation in ABCD parameters is 


where, 

3 




A 2 

eq 

B 3 

eq 

f 

4 

4 


c 3 

eq 

D 3 

eq 


4 

■*— 


“ * 

— - 


i— — 


(3.91) 


A 


eq 


jx 

tl 


[H] fc [A n + Z^ C n ] [H] 


(3.92) 


g- _ a ' a " 1 r.r^r / „n . r „n ^n „ „n .„n. „n 


3 

J eq 




[H] [ ( A n +[Z^J C n ) Z^+B n + Z^ D“][H] 

(3.93) 


D 


ec * ct 

3 

'eq 


I I I rr 

a nm 


1 [ H] C n [ H], 


a i b~ [ H J fc [D n -H2 n Z^ 2 ] [H ] 


Ct 




(3.94) 


(3.95) 


and 


a = 


a ^- and a”= “ 3t2 


a 


ntl 


a 


(3.96) 


nt2 


When the leakage impedances of the transformers and T 2 are 
neglected, the following equivalent three-phase parameters are 

obtained 


a L ■ 

eq 

-At- ~ [H] t A n [h] 

ct % 

(3.97) 

3 

B = 

eq 

ex' a" -L- [H] B n [H] 

(3.98) 



Ill 


eq 


i ii 

d a 


T£[Hf c n [H] 


(3.99) 


D_ 


a. 


[H]" D n [H] 


eq 


a 




(3.100) 


Equations (3.92)-(3.95) give the equivalent three-phase ABCD 
parameters in general. It is important to observe that the 


equivalent three-phase ABCD parameters can also be obtained by 
considering off-nominal turn ratios of the transformers, ex- 
pressing the versality of this representation. Further, it 

3 3 

is seen that A and are not equal due to the off-nominal 

eq eq ^ 

turn ratio and the different leakage impedances of the trans- 
formers. 


If the tapping on the primary side of transformer 


ii 


is only considered, d = & 3tl and a =1 can be substitu- 
ted in the equations (3.92) -(3.96) to obtain the parameters. 


3.6 EQUIVALENT SINGLE-PHASE REPRESENTATION OF A MULTI-PHASE 
(n-PHASE) TRANSMISSION LINE 

The single-phase equivalent of a n-phase line for 

the analysis of a balanced system can be conveniently derived 

from the equivalent three-phase parameters in equations (3.92)- 

(3.95). Here, A^ # , cJL and D 3 are the square diagonal 

eq eq eq 

matrices with A*, B \ , C* and as the diagonal entries 


respectively • 
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Let Z p #G g and be the equivalent single- 

phase series impedance and shunt admittance of a n-phase line 

1 1 

respectively, and Z p ,Y gh be the series impedance and shunt 
admittance per phase of n-phase line, then the 71 -representation 
holds the following relations. 


A eq " 1 + I Z P, eq Y sh,eg 


and B 1 = zt 

eq P,eq 


(3.101) 

(3.102) 


1 1 

The values of A,,, and B_ a are evaluated from equations (3.92) 

eq eq 

and (3.93) and are as follows. 


,1 _ a. 


oq 


a 


t i+ § z l ^ +z A^ lw - 25 zl Ah> i ' 3 - io3 > 


B 


eq 


a , a"[z 3 -KZ 1+ z 2 )Cl44^ h ) + z 1 z 2 ^ h (l-K3.254^ h ) ] 

(3.104) 


1 1 

From equations (3.101) - (3.104) , Z p and Y sh can be related to 

zi and Y 1 , by the following relations. 

P,eq sh,eq J 


Z p,aq = aV [2p + ^l +Z 2>< 1 4 Z P^h )+Z l Z 2 Y sh< 1 - K5 - 25Z P Y sh , 3 


(3.105) 


and 


1 4 2 
sh,eq 


Z 


[<g-^ 


(3.106) 


P,eq 
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Thus equations (3.105) and (3.106) give in general the single- 
phase equivalent of the n— phase transmission line in terms 
of n -circuit representation. 


If the transformers and T 2 are considered with nominal 
turn-ratios, the equations (3.105) and (3.106) reduce to 

Z P,eq = Z p + < Z 1 +Z 2> ( 1+0 - 5Z P Y sh )+Z l Z 2 Y ih (1+O - 25Z P Y sh ) < 3 - 106a > 

and 


Y 


1 

sh,eq 



[ Z P Y sh +(Z l +Z 2’ Y sh (lH<, - 25Z P Y sh ) ] 


(3.106b) 


It is important to mention that the equivalent single-phase 

3 3 

tt - circuit parameters have been obtained using A eg and in 
equations (3.92) and (3.93). But there are 4 equations of 
equivalent three-phase parameters, and any two equations can 
yield the desired single-phase jr-circuit parameters from them. 

Again, using equations (3.93) and (3.94) for Tt-repre— 
sentation,the following relations are obtained. 


eq 


Y sh,eq (1+ °" 25 4,eq Y sh,eq> 


(3.106c) 


and 

B 1 = zi as in equation (3.102) (3.l06d) 

eq P,eq 

Prom equation (3.102) and (3.l06d), equation (3.106a) is 
obtained from equation (3.94), we get 

C tq ■ ^h (1+0 - 25 44> 


(3.106e) 
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From equations (3.106c) and (3.107e), the following quadratic 

equation is obtained. 


^4W4h,eq-4 1+0 - 25z p4> " 0 


(3.106f) 


which yields 


4,eq 


—1 ■+• Y* [ l -f*z 


< Z P,e</ 2) 


(3.106g) 


using Binomial expansion and neglecting the 3rd and the 

higher order terms, we get 


4.. q = [| z P.,c/sh <x+ °- 25z P%;>:i /(z P,e</ 2 > 


4 1h -°- 252 £4 


(3.196J) 


Equation (3.106J) can also be expressed as 

1 1 r ~.l 1 r-.< 1 r« . ~ r 1 r- rJ" 1 * . r-r t-> / 1 'I 2 

Y sh,eq “ ~ 

P,eq 


[ Vsh +(z i +z 2 ) y sh (1+ °- 25 Vsh 5 +z i z 2 ^sh 5 


(l-to^sz^i+fz^Zjix.s Zp l 4 2(1+0 - 25z ^ y 5h )] 

(3.106 k) 

If the u -circuit parameters are derived from equation (3.92) 
and (3.94) the single-phase equivalent parameters will be 

Zp.cq = Z p + 


(3-.1061) 
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and 


ysh, oq = jr— [ 4^h +2z iysh (1 ' K3 - 5z pysh ,+z i<ysh ,2(1+0 - 25z P^h ) 

P#eq 

+ - 25 z p y sh ( z P y sh + 2 z iysh ( 1 + 0 - 5 z pyth ,+z i ( ^h >2(1+0 - 25 


z p 4h> ) ] 


(3.106m) 


= y«h( 1 + 0 - 25 z pyL ) 


(3.l06n) 


If the co rr o sponding parameters are derived from equation (3.95) 
and (3.96), those can be obtained by substituting Z 1 = Z 2 in 
equations (3.3061) and (3.106m). 

The pairs of equations (3.106a) -(3. 106b) , (3.106a) and (3.106j) 

or (3.106k), and (3.1061) and (3.106m) or (3.l06n) give the 

equivalent* single-phase n -representation of a balanced n-phase 

transmission line with nominal turn ratios. It is observed 

from the equivalent single -phase inpedance in equations (3.106a) 

and (3.1061) that the value differs only due to difference in 

the leakage irrpedances of the two transformers. If the 

1 

leakage impedances Z^ and Z 2 are taken as equal, Z p /S q be 

the same by all the pairs of equations considered. 

Prom equations (3.106b) , (3.106k) and (3. 106m) > it 
is observed that there is a difference in the calculated 
value of y 1 , from tho different pairs of equqtio»s, but 
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from the practical point of view, the equation (3.106b) can 
bo considered without appreciable error. The equivalent 
admittance can also be calculated from equation (3.l06j) or 
(3.106n) . 

If the transformers and T 2 having equal leakage 
impedance are assumed to operate with nominal turn ratio, 
the equations (3.106a) and (3. 1 06b) yield 


Z. 


P,oq 


zJ; + 2Z(l40.5^ h )«; 2 y^ h (lH0.2 5 ^ h ) 


and 


y sh,oq 


1 

P,oq 




(3.107) 


(3.108) 


3.7 EQUIVALENT n-PHi*SE REPRESENTATION OF A THREE -PHASE 
TRANSMISSION LINE 

If the analysis is carried out on n-phase basis of 
the composite system, it is necessary to derive the equivalent 
n-phase parameters of the three-phase transmission line inte- 
grated with n-phase/ three-phase transformers at the each end 
of the line. Thus, the equivalent n-phase of the line has 
been derived in terms of phase impedance parameters and also, 
in terms of ABCD parameters. 
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3.7.1 Equivalent n-Phase Impedanes of a Three-Phase line 

In order to derive the equivalent n-phase impedance,, 
et us consider a three-phase line connected to the n-phase 
buses S and R via n-phase/ three-phase transformers and T£ 
as shown in Fig. 3.8. 

From the Fig. 3.8, the following relations are 

obtained. 


" v s> ^ z tJ I s' 

using equation (3.51), we get 


,n „ „n 


V 


a 


ntl 


a 


3tl 


[H] V gll + [Z^j] l£„ 


and 


R 'R« “ L *t2 J X R' 




a 


[«] V p„ - C Z -t2-l I R" 


ht2 

a 3t2 L ^ J R 


(3.109) 


(3.110) 


Subtracting equation (3.110) from (3.109) and using the des- 
cription of three-phase transmission line as 

V S„ - V R" =[ Z p] 4" (3.111) 

Prom equations (3.109) and (3.110), and using equation (3.111), 


wo get 

a 3tl v n 
tt ntl S 



( 3 . 112 ) 



11 
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Assuming nominal turn-ratio of the transformer T ± and ? 2 , the 

following expression is obtained. 

Z P,aq = \ + [ z Sl +Z «] 0.113) 

Equation (3.113) represents the equivalent n-phase impedance 
of the three-phase transmission line. 


3.7.2 Equivalent n-phase ABCD Parameters of a Three-Phase 

Line 

The equivalent n-phase ABCD parameters of a three- 
phase transmission line can be derived as follows. Let us 
consider o. threw -phase line connected between two n-phase buses 
S and R via two three-phase/n-phase transformers as shown in 
Fig. 3.8. 


Then# the following relations can be obtained from the Fig. 3*8. 

v s =v £' I S' 

using equations (3.50) and (3.51)# 


V 


.n 




3tl 


(3.114) 


and 


u. 

T n „ 

S ~ U 


3 tl • r r t T -r 2 

[ H j I s ,, 

htl 


(3.115) 
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From equations (3.114) and (3.115), and using equation (3.96), 


4 


-It C h ] 

“'[z£j] [h] 

| 

V S" 

4 


0 

a- [H] 


4- 


Similarly, the voltage and current relations at R,R' 


(3.116) 


and R* ' 


are 



vj? i-[ 

R L t2 


3 1 


n 

R 


or 


and 


a 


ht2 


a 


3t2 


r«] v r.. - ^ + t z t 2 ] 4 


<x. 


V. 


3t2 


R“ 


Cl 


ht2 


1 tH]" 

ntz 




(3.117) 


a 


l r* 


nt2 1 r ui 

*» [KI 


R 


(3.118) 


From equations (3.117) and (3.118) 





I R ,. 



£11 




L J 
( 3 . 119 ) 
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The three-phase line connected between S M and R u can be 

modelled as 


I 

ro to 

i > 

» 

A 3 

B 3 1 


V R " 

3 ! 


_3 

3 


3 

I s" 


C 
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i R „ 



— 






(3.120) 


From equation (3.116),. (3.119) and (3.12o), the following n- 
phase equivalent is obtained. 
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v s 
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Xe 
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'eg 
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eq 


-h 

"eq 
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eq 


v 11 
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( 3 . 121 ) 


where , 


A Sg “ A 3 [H] t + CCi« [Z^ 1 lH]C 3 tH] t ] (3.122) 


B' 


n 

"eq 


A 3 [H] t [ Z^j+a'ot" [z£j[Hlc 3 [H] t [z'j2] 

3 r tti t 

, - I I I— I Iff 

a” L4 -tl 


+ [H]B 3 [H] t + # f [zj^] [H]D 3 [H]' t ] 


,n 

"eq 


n 


ot'a" [H]c 3 [H] 

^ITl 


(3.123) 

(3.124) 


D_ 


eq 


i.[ _|L.[h] d 3 [h]* + «" Ch] c 3 [ H f [z “ 2 ]] 


(3.125) 


Thus, equations (3.122) - (3.125) give the equivalent n -phase 
ABCD parameters of the three-phase line. 
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3.8 CONCLUSION 

rh© m atnem atical model of a general n— phase (multi- 
phase) synchronous machine has been developed in phasor co- 
ordinates. that will help for its steady-state analysis both 
under balanced as well as under unbalanced conditions. Further/ 
three-phase/n-phase transformers having off-nominal turn 
ratio on both the sides have been modelled in order to analyse 
the combined three-phase and multi-phase systems by consider- 
ing two types of three-phase windings, namely, wye and delta, 
and one type of multi-phase windings i.e., star connected 
winding. Further, the multi -phase transmission line has been 
represented in terms of phase impedances, m-circuit admittance, 
and ABCD parameters that will help in analysing the short, 
medium and long lines respectively. 

The three-phase equivalent of a n-phase transmission 
line integrated with three-phase/n-phase transformers has been 
obtained in order to carry out the analysis on a composite system 
because of loss computational efforts, when the interest of 
investigation lies on the three-phase side of the system. For 
the analysis of a balanced composite system in order to carry 
out the investigation on single-phase basis, the equivalent 
single-phase parameters representing n-circuit have been 
obtained based on the three-phase equivalent ABCD parameters. 

The equivalent n-phase parameter for the three-phase line 
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has also been developed in terms of impedances and ABCD para- 
meters from which the single— phase equivalent can conveniently 
be derived. This equivalent n-phase parameters will be needed 
when the three-phase line is connected to n— phase lines by 
means of three-phase to n-phase power transformers and the 
interest of investigation lies in the n-phase line of the 
composite system. 



CHAPTER 4 


LOAD PLOW ANALYSIS OF SIX-PHASE (MULTI-PHASE) 
POWER SYSTEM NETWORKS 


4.1 INTRODUCTION 

In the previous chapter/ mathematical discrip- 
tions of different conponents of power system networks have 
been derived. This chapter presents the load flow analysis 
of a combined three-phase and multi-phase systems with the 
help of the equivalent single-phase parameters of a multi- 
phase transmission line connected to three-phase/n-phase 
power transformers assuming the power system to be completely 
balanced/ in order to compare the performance of a multi- 
phase transmission line relative to that of the conventional 
three-phase system. 

There is a vast literature available for a three- 
phase power system load flow analysis [ 25/ 62-67]. With the 
emergence of the idea of multi-phase transmission line, a 
load flow study has been carried out on mixed three-phase and 
six-phase power systems on the single-phase basis [2/2l] / and 
also In the phase frame of reference [ 30] in order to compare 
the performance of the mixed system to that of completely 
three-phase system having double circuit three-phase lines. 

In the analysis/ the phase to ground voltage of a six-phase 
line has been assumed to be equal to V*3 times that of the 
three-phase line. This chapter is also an attempt in this 
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direction by presenting the load flow investigation for the 
different situations under different voltage conditions 
namely , the phase to ground voltage of the multi-phase line 
is V3 times or equal to that of a double circuit three- 
phase line, when a double circuit three-phase line is con- 
verted to a single circuit multi-phase line. 

The conversion of a double circuit three-phase 
line to single -circuit multi -phase line will essentially 
require a throo-phase/n-phase power transformer at each end 
of the multi-phase lino. In the cases, as mentioned above, 
tiie conversion requires for three-p ha se/n -phase transformers, 
the varying voltage on the n-phase side ( f3 times or equal 
to that of the double circuit three-phase line). Therefore, 
in order to analyse such a balanced system under different 
voltages of three-phase and multi-phase lines, the equivalent 
single-phase parameters of a multi -phase line integrated with 
two three-phase/n-phase transformers at the ends incorporating 
the different voltages on n-phase side have been developed in 
section 4,2, Different types of problems have been formulated 
for the investigation of load flow analysis in section 4.3 
along with the discussion of the load flow methods required 
for the study. Different cases for the investigation have 
been taken on two sample networks in section 4,4. In section 
4.4.1, a single machine 3-Bus system has been taken to 
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study the load flow performance for the different cases in 
order to cortpare the performance of systems, namely a mixed 
three-phase and multi -phase transmission system and a comp- 
letely three-phase transmission systems; a multi-phase trans- 
mission line is obtained by replacing a double circuit three- 
phase line. 

The same problems of load flow analysis have been 
investigated in section 4.4.2 on a three-machine nine-Bus 
power system network by replacing three double circuit three- 
phase lines to single circuit multi-phase lines one by one. 
The results of the investigation on two sample networks for 
different cases have been presented and discussed in section 
4.5. Finally, the chapter ends in section 4.6 with the 
conclusions based on the discussions of the results. 

4.2 EQUIVALENT SINGLE -PHASE REPRESENTATION 

This section derives the equivalent single-phase 
p arametors for n -representation of a multi-phase transmission 
line integrated with three-phase/n-phase transformers at the 
ends of the line to facilitate the analysis on a single-phase 
basis of 'the combined three-phase and multi-phase power 
system networks. With the help of a three-phase/n-phase _ 
transformer, the power can be transmitted through multi-phase 
transmission lino at the varying voltages, and therefore, the 
expressions for equivalent single-phase series impedance and 
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shunt admittance have been derived for general case inclu- 


ding the number of phases and the different voltages on n— 

phase side of the transformer. 


A multi-phase transmission line is connected to 
three-phase buses through three-phase wye/n-phase star connected 
transformers as shown in Fig. 3.7. 


From Fig. 3.7/ the voltages at the sending end 

can bo expressed as. 




(4.1) 
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(4.3) 


Similary, the voltages at the receiving end are re- 

la tod as follows. 
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Prom equations (4. 2) -(4. 3) and (4. 4-4. 5), the following 

representations in matrix form are obtained. 



(4*6) 




where 


3 _ Ju- !^3tl . 5a£2. [uf [A n +[Z^]]^][N] 

ec 5 Km a ntl # a 3t2 


(4.8) 


(4.9a) 


3 „ 3 _1_ ^ B n + [z^D n ][N] 

eq n a ntl %. t 2 • - - ' 

r 3 n 1 a ntl [n] 

C eq ~ 3 ct 3tl * a 3 t2 


(4.9c) 
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aT 77 M t [D n +c n [< 2 ]]' ) 


(4.9d) 


3,3 3 3 

A, eq / ^eq # C eq an< ^ D eq are ' t ^ le ^i a 9°nal matrices 
111 1 

which have A e ^, B eq' ^-eq anc ^ D eq as t * le: *- r diagonal elements 
respectively. Prom equations (4,9a)-(4.9d) , the respective 

values are 


a 1 - U 3tl a nt2 r . 1 11 1 1 i 

A eq ” a ntl * u 3t2 £ 1+ 2 z p y sh +z l y sh (l40 * 25z p y sh 5 J 


(4,l0a) 


,1 _ 1 U 3tl St2 r 7 1 , W1 1 _1 1 _ 1 n L . — .I 1 . 

3 w i [ Zi 4* \Zj -hZ 0 ) ( 1 ~H=r Z y , )4~Z 1 Z 0 y _ (14-0.25Z, y , ) J 

eq n u- ntl ^ nt2 p 1 2 2 P stl 1 2 J sh' p-^sh 7 J 

(4.10b) 


n a ntl a nt2 1 M ,1 11 . 

3 * a 3t2 ^ J? y slv 


(4.loc> 


c 1 4<&ih + Vsh t 1 "- 25 %V sh > ] w.iod 


1 1 

where, Z p and y g ^ are the series impedance and the shunt 
admittance of a multi-phase transmission line, ana Z^ and Z 2 
are the diagonal elements of diagonal matrices [ z t ^ J and - 
[ ] respectively, representing the leakage impedances of 

the trims formers and T 2 . 

If the multi-phase transmission line is represented by equi- 
valent or nominal ft -circuit parameters, the equivalent single - 
ph ase series impedance and the equivalent single-phase shunt 
admittance can be known with the help of ary two equations 
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among (4.10a) - (4, lod) . 

Using equations (4.l0a) and (4. lob) , we get 

a. 


ct 

3 3tl 


3t2 


J P / eq n ct. 


ntl 


01 


ht2 


[z p + 


+Z l Z 2Vsh tl+0 - 25Z ^> , sh)] (4.U) 


and 


Y sh, 


cq 


— + 2Z iysh (1+0 * 25z p4> 3 


T 

P,oq 


(4.12) 


using equations (4.10b) and (4.10c), 

1 

2p is iound as in equation (4,11) 
and 

*ih,oq ■ -^r 1 - [z p4 +(z i +z 2>4 (1 «- 25z p4 )+z i 2 2(ysh) : 


1 r_1 

. 11 ,. ,,, V ,„1 , 1 v 2 . „„_1 1 


P#oq 


(l+0.252^ h ) + (Z 1 -HZ 2 ).5Z p (yJ h ) (l-+0.25Z^y sh )] (4.13) 


or 
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n a ntl 


sh, cq 3 a. 


3tl 


-2&L y 1 (1+. 25 Zp y*> 
a 3t2 Sh P Sh 


(4.14) 


using e qu a ti ons (4,10a) and (4.10c), 

y 1 . is found as in equation (4.14) 

J sh,oq 

,1 ..3 Stl “3t2 [ Z 1 + az^l+O^SZ^-.SZ^y^ 


P # cq n 


a ntl a nt2 


(l+0.25Zpy^ h ) ] 
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From equations (4.11) and (4.15) which give the equivalent 
single-phase series impedance from the two sets of equations, 
it is observed that the two values are different. But for 
the practical purposes, where Z ± and Z 2 do not differ much, the 
values obtained from the two equations are almost the same, 
since the third term under bracket in both the equations con- 
tribute insignificantly towards the value. It is also 

observed that the equivalent single phase series impedance 
1 

Zp OCI decreases with the increasing number of phases while the 
shunt admittance from equations (4.12), (4.13)- and (4.14) 
increases with the increasing phases. The shunt admittance, 
from equations (4.12) and (4.13), is almost the same because 
the 3rd and 4th term under the bracket in equation (4.13) 
adds little to its value. 

4.3 PROBLEM FORMULATION 

Suppose there are double-circuit three-phase lines 
in a power system network which is generally true. There is 
a problem of transmitting more powers to meet the increased 
demands at the receiving ends and let us assume that there 
are two alternatives to accomplish this task viz. 1) the ^ 
transmission voltage of the existing double circuit three- 
phase line is increased, say by 73 percent; 2) the existing 
double circuit three-phase line is converted to a single -circuit 
six-phase line making its phase to ground voltage equal to 1/3 
time that of the existing double circuit three-phase line i.e. 
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maintaining the same line to line voltage, in order to 
compare the two options from the viewpoint in two cases of 
power transfer capability* load flow investigation under two 

conditions are desirable. 

Case a : when the slack bus voltage is fixed 

The load flow analysis usually requires one slack 
bus which is a reference bus. As the investigation of the 
load flow because of converting a double circuit three-phase 
lino in either of the options can be interpreted as an equi- 
valent to changing parameters of the line from one value to 
the another and its effect on the performance thereof; so by 
assuming slack bus* the slack bus complex generation will 
change correspondingly from one parameter to other parameter . 
Since load flow is carried out for planning* it necessarily 
gives the slack bus generation, and also the other performances 
with the changing parameters. 

Case b i when the complex generation at the slack bus in case 

a is fixed 

The load-flow investigation can* however* be 
viewed in another way. For example* without affecting the 
generation of any bus connected to the generator* the effect 
in the steady-state conditions of the power system network for 

the two options can be investigated. 
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There are , however , two types of situations 
where a single-circuit multi— phase line can be considered, to 
repl ace a double circuit three-phase transmission line* 

1) if the power transmission capability of the existing double- 
circuit three-phase transmission line is required to be 
increased to meet the growing demand. This can be accoirp- 
lished by converting the double-circuit three-phase line to 
a single -air cui t six-phase line with its phase to ground 
voltage equal to V 3 times that of the double circuit three- 
phase lines i.e. maintaining the same line to line voltages 
in both the cases. 

The investigation will facilitate to examine the 
performance of multi-phase transmission lines transmitting 
73% more power in comparision to a double circuit three-phase 
lines, and also the performance of the multiphase lines when 
the same power is transmitted at V" 3 times voltage (phase to 
ground) of the double circuit three-phase lines. Since the 
voltage between two adjacent phases of multi -phase line has 
been considered to be equal to (in case of six-phase line) 
and less than (in case of phases higher than six) that of a 
double circuit three-phase line, the adjacent phase to phase 
insulation requirement will be the same for both systems, but 
the phase to ground Insulation requirement will be more in 
case of multi-phase systems than that in double circuit three- 
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phu t ^e lin- ii> % 13 u t the insulation is generally provided con- 
servatively and therefore, tlis multi— phase uprating can talcs 

pi ace without restructuring. 

2) If an existing double circuit three-phase line is rep- 
laced by a single-circuit multi -phase line for transmitting 
the same power at the same phase to ground voltage retaining 
the advantages of the multi -phase transmission line compared 
to a double circuit three-phase line as pointed out by the 
investigations 1 . Since the phase to ground voltage is the 
same! in both cases, so the phase to ground insulation require- 
ment is tile same whereas the adjacent phase to phase insula- 
tion requirements for multi-phase system will be equal to 
1/ Y 3 times in case of six-phase transmission system, and 
even less, in case of transmission line higher than six, 
than double circuit three-phase lines. This investigation 
will examine the performance of a multi-phase transmission 
line in comparision to a double circuit three-phase line at 
the same phase to ground voltage of the two transmission 
systems. 

The investigations of the system for the two 
voltage conditions have been carried out under case a and 
case b. The load flow analysis pertaining to case a is carried 
out by different methods , such as Gauss-Siedal, Newton-raphson 
method, in which one of the buses is assumed as a slack bus 
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(l«o# re for once bus) whose voltage’ magnitude and. angle are 
known, and the voltages of the other buses as well as the slack 
bus generation cure calculated. But the load flow analysis 
pertaining to case b is hardly discussed [ 68] . In case b, 
the reactive power generation instead of the voltage magni- 
tude is known at the slack bus. The angle of the slack bus 
voltage is known as in case a, and all the bus voltages are 
calcula tod. 

4.4 CASE E TO DIES 

This section investigates the load flow analysis 
of a completely three-phase systems, a combined three-phase 
and multi-phase systems to compare their performances under 
the different cases on the two sample networks. 

4.4.1 A Three-bus Single-Machine Sample Network 

The load flow analysis for the two cases has been 
carried out on a three-bus single machine sample network [ 2] 
as shown in Fig. 4.1. The necessary data are provided in 
Appendix A.l. The pu parameters of different types of trans- 
mission lines have boon calculated and have been given in 
Table A. 2. The p.u line parameter for the different confi- 
gurations, namely, a double-circuit three-phase line, a single 
circuit six-phase line; the phase to neutral voltage of the 
multi-phase line is equal to, or V " 3 times that of the double 
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F 16.4.1 A thrct-bus Single- machine Sample 
Network. 
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circuit three-phase lines, an* the corresponding parameters 
are provided in Table 4.1. Two values of each parameters 
have been given by calculating them with two sets of equations, 
namely, agnations (4.11) and (4.12) and (4.14) and (4.15) and 
the values are almost identical. The line number between 
buses 1 and 2 is converted to a single-circuit six-phase and 
twelve-phase transmission lines for two configurations, 
namely , the phase to ground voltage of a multi -phase line is 
]f 3 times that of the double-circuit three-phase line and 
the phase to ground voltage of the two systems are the same, 
under casus a and b. in Case a, the system has also been 
investigated with the increased system loading. The leakage 
reactance of the transformer at each end of the multi-phase 
line has been taken as 8%. 

4.4.2 A Nine Bus Thr e ee -Machine Sample Network 

A sample network [ 55] , shown in Pig. 4.2, has 
been taken, where line nos. 2,3, and 6 have been assumed as 
double-circuit three-phase lines, to investigate the inpact of 

conversion of those three-phase lines to multi-phase lines for 

of 

the purpose/ studying the load flow performance of the network. 
The necessary data of the sample system for a completely three- 
phase lines have been provided in Appendix B. The pu series 
impedance and the shunt admittance of the multi-phase line 
connected to an appropriate transformer at each end of the 
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line have been provided in Table 4.2. Two values of the 
leaka g ' roue buncos (i. e. , 8% and 5%) of the transformers have 
been taken to obtain each set of parameters. For, the present 
load flow study, the p.u. equivalent single-phase series 
impedance and the shunt impedance corresponding to 8% leakage 
Impedance of the transformer have been used. The load flow 
investigation has been carried out for case a and case b. 
Firstly, the lino no. 2 (double-circuit three-phase line) is 
converted to a multi-phase line and thereafter line nos. 3 
and 6 are also converted to multi-phase lines to investigate 
the Impact of conversions of double-circuit three-phase 
lines to multi-phase transmission lines. 

4.5 RESULTS AND DISCUSSIONS 

The load flow study has been carried out on a 
sample system shown in Fig.4,lwhere line no.l is a double- 
circuit three-phase line, to investigate the impact of con- 
version from double circuit three-phase line to a single- 
circuit six-phase or to a twelve-phase line under voltage con- 
dition A and B for case a and case b. The bus voltages and 
the line flows for a six-phase conversion have been provided 
in Table 4.3. and 4.4. for case a, and those for a twelve- 
phase conversion in Table 4.5 and 4.6. The bus voltages for 
case b of a six-phase and a twelve-phase conversions have been 
provided in Table 4.7, while the line flows and efficiencies 
in Table 4.8. 
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TABLE 4.3 


for 


* 3us ^tages of the Sample Network (Fig.4.1) 

Uioe a (A double Circuit line (no.l) is converted 
to a single-circuit six-phase line) 


BUS 

NO. 

Type 

of 

Bus 

Generation 

Load 

Voltage 

magnitude 

Angle in 
degree 

Voltage 

condition 

1 

Slack 

3 31. 4 3+j 11 3. 81 

O.O+jO.O 

1.07 

0.0 


2 

P/Q 

O.O+jO.O 

190+j 40 

0.9864 

-10.25 

I 

3 

P/Q 

O.O+jO.O 

130+j 25 

0.9729 

-12.5041 


1 

Slack 

325. 79+j68„41 

O.O+jO.O 

1.07 . 

0.0 


2 

P ,Q 

O.O+jO.O 

190+j 40.0 

1.0329 

-6.8129 

A 

3 

P/Q 

O.O+jO.O 

130+j 25 

1.0006 

-10.5428 


1 

Slack 

3 33. 44+j 162.77 

O.O+jO.O 

1.07 

0.0 


2 

P/Q 

O.O+jO.O 

190+j 40 

0.9229 

-18.393 

B 

3 

P/Q 

O.O+jO.O 

130+j 2 5 

0.9288 

-16.7839 


1 

Slack 

4 70+j 151.97 

O.O+jO.O 

1.07 

0.0 


2 

P/Q 

O.O+jO.O 

330+j 70 

0.998 

-10.9046 

II 

3 

P/Q 

O.O+jO.O 

130+j 25 

0.9786 

-12.7642 



* In the last column, the voltage condition, refers to the 

following * 


I - a double circuit three-phase line 

II - a double circuit three-phase line is converted to 

a single-circuit multi -phase line under voltage 
conditions A with increased lo.ad. 



TABLE 4.4 
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Line Flows and Transmission Efficiencies of the Sample Sys- 
tem shown in Fig. 4,1 for Case a 

(Line No.l is converted to a single-circuit 6-Phase 

line) 


Voltage 

Condition 

Line 

No. 

From 

Bus 

To 

Bus 

Line 

P ij 

(MW) 

flow 

Q ij . 

MVAR 

E f f iciency 

% 

I 

1 

1 

2 

215.91 

72.02 

96.68 


2 

2 

3 

18.76 

1.46 

99.37 


3 

1 

3 

115.52 

41.78 

96.39 

A 

1 

1 

2 

226.76 

41.62 

98.89 


2 

2 

3 

34.26 

9.29 

98.87 


3 

1 

3 

99.02 

26.79 

97.06 

— — *-* — 

— 

— 

— 




B 

1 

1 

2 

184.25 

94.71 

96.75 


2 

3 

2 

11.77 

1.18 

99.57 


3 

1 

3 

. 149.188 

68.06 

95.02 

11 

1 

1 

2 

352.81 

112.92 

98.17 


2 

2 

3 

16.37 

4.59 

99.39 


3 

1 

3 

117.98 

39.04 

96.38 
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TABLE 4 . 5 

Bus Voltages of the Sample System Shown in Fig. 4.1 for Case a 

(Line No.l is converted to a single circuit 

12 -phase line) 


BUS 

NO. 

" Type of 
BUS 

Generation 

Load 

mag. 

Angle 

in 

degree 

Voltage 

Condi- 

tion 

1 

Slack 

324 « 11 + J 13*69 

O.O+jO.O 

1.07 

0.0 


2 

P,Q 

O.O+jO.O 

l90+j40 

1.0573 

-3.5699 

A 

3 

P,Q 

Q.G+JD.O 

1 30+j 25 

1.0152 

-8.7769 


1 

Slack 

327.S5+J97.84 

O.O+jO.O 

1.07 

0.0 


2 

P,Q 

O.O+jO.O 

I90+j,40.0 

1.0073 

-9.8558 

B 

3 

P#Q 

O.O+jO.O 

1 30+j 25.0 

0.9846 

-12.1965 


— 

— 

— 





1 

Slack 

466.38-fj69.31 

O.O+jO.O 

1.07 

0.0 

. ' i 

2 

P,Q 

O.O+jO.O 

3 30+j 70.0 

1.0428 

-5,5717 

■ 11 

3 

P/Q 

O.O+jO.O 

1 30+j 25 

1 .0066 

-9.8657 

1 


TABLE 4.6 

Line Plows and Transmission Efficiencies of the Sample Shown in Pig. 4.1 for Case 

( £ ine No.l is converted, to a Single— circuii. lz— phase line) 
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TABLE 4.7 


Bus Voltages of the Sample System (Pig. 4.1> for Case b 


BUS 

NO. 

Generation 

Load 

Voltage 

Angle Volt- 
in age 

degree Condi- 
tion 

* 

Type 

of 

Multi- 

Phase 

Line 

1 

3 31 . 39+ j 113.32 

O.O+jO.O 

.8935 

0.0 



2 

O.O+jO.O 

190+j 40 

.8362 

-10.14 

A 


3 

1 

O.O+jO.O 

331,39 H 113. 32 

130+j 25 

O.O+jO.O 

0.7910 

1.2415 

-16.175 

0.0 


six- 

phase 

line 

2 

O.O+jO.O 

190+j 40 

1.1401 

-12.81 

B 


3 

O.O+jO.O 

130+j 25 

1.1434 

-11.75 



1 

331. 39 +j 11 3. 32 

O.O+jO.O 

.7396 

0.0 



2 

O.O+jO.O 

190+j 40 

.7028 

-7.7814 

A 


3 

O.O+jO.O 

130+j 25 

.6156 

-21.0299 

12- 

— Phase 

1 

331.39+JH3.32 

O.O+jO.O 

1.0103 

0.0 


line 

2 

O.O+jO.O 

190* j40 

0.9398 

-11.27 

B 


3 

O.O+JO.O 

130+j 25 

0.9129 

-14.02 




* Type of multi -phase line refers to a multi-phase transmiss 
ion lino to which a double circuit three-phase line (No.l) 


has boon converted* 
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TABLE 4.8 


Line 

flows < 

:and Tj 

ransmission Efficiencies of 
(Fig. 4.1) for Case b 

the Sample 

Network 

Line 

No. 

From 

Bus 

To 

Bus 

Line 

P 

(MW) 

flows 

Q 

(MVAR) 

Effici 

n 

% 

. Voltage 
condition 

Type of 
multi- 
phase 
line 

1 

1 

2 

229.83 

73.0 

98.29 



2 

2 

3 

35.198 

12.976 

98.16 

A 


3 

1 


100.84 

40.37 

95.37 


six- 

— 

* — 



— 



- phase 
line 

1 

1 

2 

183,23 

66.29 

97.85 



2 

3 

2 

11.57 

3.196 

99.72 

B 


3 

1 

3 

147.32 

47.1 

96.67 



1 

1 

2 

245.386 

68.69 

98.68 



2 

2 

3 

51.97 

26.09 

95.82 

A 


3 

1 

3 

85.79 

44.64 

93,71 


12- 



— 


. — 

— 



Phase 








line 

1 

1 

2 

216.094 

72.44 

88.1 



2 

2 

3 

21.15 

6,51 

99.1 

B 


3 

1 

3 

114.46 

40.95 

96.0 
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From Table 4.3/ it is observed, that/ when the slack 
bus voltage magnitude is fixed at 1.07 p.u., the reactive 
power at the slack bus are different for the different confi- 
gurations. For example/ the required generation for an 
original (i.e. completely) three-phase system at the slack 
bus is 113.81 MVAR while it reduces to 68.41 MVAR for voltage 
condition A and increases to 162.77 for voltage condition B 
and 151.97 for configuration II. It is also observed that- 
the voltage and angle improves under configuration A and 
they deteriorates under configuration B in comparision to 
configuration 1, while in configuration II the voltages and 
angles are comparable to those under configuration I. When 
line No.l is converted to a twelve-phase line, it is observed 
from Table 4.5 that in all configurations, the voltage and 
angle improves in comparision to those in configuration I, 
bu t the voltage and angle in condition A & II have improved 
more than those of condition B, a similar observation as 
from Table 4.3, and also, the reactive power required under 
configuration A & II are less than those under configuration! 
B } however, for this twelve-phase conversion, the reactive 
power in all co n f igura t ions are less than that of configu- 
ration I. 

Prom Table 4.7, where, the load flow study has been 
carried out by assuming the reactive power (though normally 
not true) of the slack bus (generation) fixed at 113.32 MVAR, 
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it is observed that the bu.s vol tsges under voltage condition 
B for six -phase conversion are higher while those under voltage 
condition A arc lower than the voltages under configuration I. 
For twelve -phase conversion, the voltages are even lower than 
six-phase configuration under voltage condition A> and higher 
chan those for six-phase conversion under voltage condition B 
hut it is still lower than those in configuration I. 

Further, the load flow investigation has been carried 
out on a sample system network shown in Fig. 4.2 where three 
linos, namely lino nos, 2,3 and 6 have been considered as 
double-circuit three-phase lines, to investigate the system 
under different configuration especially with a view to 
examine the impact of increasing number of conversion of 
double-circuit three-phase lines to multi-phase lines for 
various cases. 

The bus voltages and the line flows of the sample 
network under configuration I have been provided in Tables 
4.9 and 4.10. The slack bus voltage has been assumed as 1.04 
for case a. The slack bus generation for six-phase and twelve- 
phase conversions under different voltage conditions and sys- 
tem configurations have been provided in Table 4,11. The bus 
voltages and the line flows for six-phase conversion have been 
provided in Tables 4,12 and 4.13 for case a, while those for 
twelve-phase conversion in Tables 4.14 and 4.15 respectively. 
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Thu bus voltages and the line flows for the sample 
network (Fig. 4 .2) have been provided for six-phase conversion 
in Tables 4.16 and 4.17 respectively for case b where the 
reactive power of the slack bus has been fixed at 27.03 MVAR, 
and those for twelve-phase conversion in Tables 4.18 and 
4.19 respectively. 

It in observed from Tables 4.9 and 4.11 that the 
slack bus reactive power are different for different voltage 
conditions, namely voltage conditions A and B, and different 
system configurations i.e., configuration I, configuration II 
and configuration III for six-phase and twelve-phase conversions. 
For example, the slack bus reactive power is 27.03 MVAR for 
a completely three-phase system, 25.44 mvar for system configu- 
ration I, 148.22 for configuration II and 246.28 MVAR for 
configuration III for six-phase conversion under voltage con- 
dition A, and the corresponding values for different system 
configurations under voltage condition B are 30.82 MVAR, ■ ■ 
45.02 MVAR and 47.3 MVAR. For twelve -phase conversion, the 
respective values arc 99.55 MVAR, 444.92 MVAR and 805.52 MVAR 
under voltage condition A and .15 MVAR, 60.84 MVAR and 107,69 
under voltage condition B. It is observed that the required 
reactive power for twelve-phase conversion is higher than 
those for six-phase conversion. Moreover, the active powers 
have also changed substantially with respect to that of the 
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TABLE 4.9 


Bus Vol fcages and Generations of the Sample System shown in 

Fig. 4. 2 ■ " ' 


BUS 

No. 

Type 

of 

BUS 

Gent 'ration 

Load 

Volt. mag. 

Angle in 
deqree 

1 

Slack 

71. 64-i-j 27.03 

O.O+jO.O 

1.04 

0.0 

2 

P, v 

16 3.0+j 6.7 

O.O+jO.O 

1.025 

9.27 

3 

P, V 

85.0+J10.9 

O.O+jO.O 

1.025 

4.66 

4 

P/Q 

O.O+jO.O 

O.O+jO.O 

1.0258 

-2.22‘ 

5 

P#Q 

O.O+jO.O 

125+j50.0 

.9957 

-3.98 

6 

P/Q 

O.O+jO.O 

90+j 30 

1.0127 

-3.687 

7 

P/Q 

O.O+jO.O 

O.O+jO.O 

1.0258 

3.719 

8 

P/ Q 

O.O+jO.O 

100+j 35 

1.0159 

0.7273 

9 

P/Q 

O.O+jO.O 

O.O+jO.O 

1.0323 

1.96 


TABLE 4.10 


Line 

Flows and 
Fig. 4.2 

Trans 
for a 

. Efficiencies of a Sample System shown in 
Completely Three-Phase Lines. 

Line 

No. 

From 

BUS 

TO 

Bus 

Line 

P ij 

(MW) 

flows 

Q ij 

(MVAR) 

Efficiency 

% 

1 

4 

5 

40.928 

22.875 

99.37 

2 

4 

6 

30.698 

1.0328 

99.458 

3 

7 

5 

86.61 

8.365 

97.35 

4 

9 

6 

60.808 

18.078 

97.78 

5 

7 

8 

76.37 

.7737 

99.377 

6 

9 

8 

24.18 

3.084 

99.63 
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TABLE 4.11 


Slack Bus Generations for Different Configurations of the Sample 

System Shown in Fig. 4. 2 for Case a 


Type of 
Line 

Voltage 

Condition 

Slack bus Generation for different configu- 
rations 


config. I 

config. II 

config. Ill 

6 -phase 

A 

71.14_j2S.44 

69.85-jl48.22 

72.57-j246.28 

line 

B 

71. 72+j 30.82 

72 ,12+j 45.02 

72. 21 +j 47. 3 

12 -phase 

A 

70. 82 -j 99.55 

76.32_j444.92 

107.19-j 805.52 

Line 

B 

71.32+j .15 

69.88-j60.84 

70.32_jl07.69 


Configuration I - when line no. 2 is converted to a single-circuit 

multi-phase transmission line 

Configuration XX- when line non. 2 and 3 are converted to single 

circuit multi— phase transmission lines 
Configuration III - when line nos. 2,3 and 6 are converted to 

single circuit multiphase transmission line. 


TABLE 4.12 
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Bus Voltages of the Sample System Shown in Pig.' 4.2 for 

Case a 

(conversion of double-circuit 3-phase lines to single circuit 

six-phase lines) 


Voltage 

Bus 

Voltage magnitude and Angle of the buses 

condition 

NO. 

config, I 
Voltage Angle 

config. II 
Voltage Angle 

config. Ill 
Voltage Angle 

A 

1 

1.04 

0.0 

1.04 

0.0 

1.04 0.0 


2 

1.074 

8.38 

1.22 

3.67 

1.3466 2.2267 


3 

1.077 

4.28 

1.199 

1.55 

1.3445 -0.1321 


4 

1.0548 

—2 .14 

1.1228 

-1.97 

1.1771 -1.95 


5 

1.0337 

-3.775 

1.174 

-3.238 

1.2698 -3.3148 


6 

1.0611 

-3.11 

1.146 

-3.22 

1.2265 -3.2484 


7 

1.074 

3.32 

1.2196 

-0.253 

1.3456 -0.9964 


8 

1.067 

0.64 

1.2060 

-1.95 

1.354 -2.45 


9 

1.084 

1.84 

1.2055 

-0.42 

1.3498 -1.7047 

B 

1 

1.04 

0.0 

1.04 

0.0 

1.04 0.0 


2 

1.0174 

9.02 

.9928 

12.06 

.9836 12.07 


3 

1.015 

4.055 

0.9917 

6.01 

.9924 6.25 


4 

1.0237 

-2.22 

1.0159 

-2.25 

1.0146 -2.253 


5 

.9917 

-4.15 

0.9774 

-4.565 

.9739 -4.63 


6 

1.0016 

-4.77 

.9833 

-4.075 

,9827 -3.95 


7 

1.0183 

3.38 

.9938 

6,134 

,9848 6.03 


8 

1.007 

0.218 

.9823 

2.47 

.972 2.24 


9 

1.0225 

1.304 

.9994 

3.13 

1.0 3.377 
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TABLE 4.14 

Bus Voltages of Sample System Shown in Pig. 4.2 for Case a 

(conversion of double circuit 3-phase lines to single circuit 

12 -phase lines) 


Volt. 

Bus 

Voltage magnitude and - 

Angle of 

the bus vchtaqes 

Cond. 

do. 

con fig. I 
mag. Angle 

con fig 
mag. 

.11 

Angle 

config. 

mag. 

III 

Angle 


1 

1.04 0.0 

1 « 04r 

0.0 

1.04 

0.0 


2 

1.1301 7.52 

1.5252 

0.0933 

1.9264 

-2.464 


3 

1.1323 3.95 

1.4526 

-0.505 

1.9397 

-3,87 


4 

1 .0958 -2.0*31 

1.2871 

-1.882 

1.4873 

-2.28 

A 

5 

1 .0825 -3.527 

1.4754 

-3.307 

1.8238 

-4.29 


6 

1.109 -2.573 

1.3247 

-2.57 

1.5714 

-3.04 


7 

1.13 2.95 

1.524 

-2.42 

1.925 

-4.03 


8 

1.124 .5802 

1.494 

-3.27 

1.962 

-4.84 


9 

1.1388 1.7344 

1.4574 

-1.85 

1.943 

-4.626 


1 

1.04 0.0 

1.04 

0.0 

1.04 

0.0 


2 

1.052 8.67 

1.1288 

5.98 

1.197 

4.83 


3 

1.054 4.29 

1.1201 

2.72 

1.1948 

1.732 


4 

1.04 -2.17 

1.0744 

-2.06 

1.1 

-2.03 

B 

5 

1.0159 -8.9 

1.0856 

-3.57 

1.1314 

-3.57 


6 

1.04 -3.58 

1.0874 

-3.66 

1.1317 

-3.62 


7 

1.0525 3.4 

1.1287 

1.4 

1 .1966 

0.753 


8 

9 

1.0443 0.56 

1.118 

-0.87' 

1.1963 

-1.31 


1.0612 1.736 

1.1267 

0.4578 

1.2 

-0.2575 



156 


in 


Hh 

a 

§ 

& 


a? 


0 

tn 

rd 

o 

u • 

0 4> 

*r4 

o 

CNI U 
* -H 

a 


Cn 

f'N 


s 

q 

P 

if) 

►h 

tf) 

<u 

**H 

a, 

a 

<0 

CO 

o 

w 

■S 

(J 

0 

-rl 

o 

-H 

m 

m 

M 


w 


<D 

i%H| 

g* 

<0 


o 

Jp 

m 

0 ***** 
a m 

•H 0 
Hf C 
•Hi 
Cl) H 
W 

tt) 
W 
0< rtf 

1 ,C 

<“0 pL( 

•P csl 

“H tH 

3 

B 

•H 

0 


5 


M 

EH 


Cl) 

<HI 

o 

'd 

l H 

0 


'd c 
o 


i 

1 

o 

H 

(h 


•H 

tn 

u 

S 

a 

o 

o 


pi 


‘a 

*H 

H-l 

S 

,° 

0 


& 

S 

•H 

0 

-H 

44 

44 

0 

r g 

c 

cd 

W) 

* * 

o 

H 

4*4 

0 

• 3 , 


cn 

•H 

44 

a 

o 

a 




%\ 

a % ] 

k 


a, | 


P* ^ 


% 

a fe 

A 


a 


g i 


i ! 


I M 

! • 

cn 

■H 

,'W 

I (H 

0 

a 


VP 1 


3 

°S 

«w* F 


ft 


CD 


PI 


0 

55 


■P T3 

0 8 

> a 


H 1 


H< 

ON 

H< 

CO 

sh 

tH 

t> 

UN 


o 

• 

• 

9 

• 

• 

ft 

CM 

V0 

VO 

o 

ON 

ON 

VO 

ON 

On 

VO 

ON 

ON 

00 




CM 

87 

H 

o 

CO 


tH 

st 

» 

9 

ft' 

• 

• 

ft 

tH 

V0 

in 

O 

ro 

ON 

H 


co 


co 

UN 

VO 

UN 

tH 

ro 

tH 

CO 

VO 

ro 

O 

ON 

CM 


co 

co 

00 

vo 

tH 

CO 

* 

• 

'ft 

« 

'ft 

ft 

CM 

sh 

tH 

LO 

tH 

ON 

H< 

VO 

o 

H< 

VO 

CO 



tH 




UN 

O 

UO 


O 


O 

ON 

H 

O 

UN 

CO 

9 

* 

• 

• 

• 

ft 

o 

CO 

ON 


ON 

ON 

t> 

ON 

ON 

ON 

ON 

ON 


CM 

CO 


ro 

oh 

VO 

co 


o 

O 

ON 

• 

4 

ft 

ft 

• 

ft 

H 

tH 

o 

o 

O 

ON 

O 

UN 


VO 

H< 

o 

no 

CM 





ON 

CM 

CO 

CM 

VO 

[> 

!> 

m 

H 

ro 

ON 

VO 

• 

» 

• 

• 

ft 

ft 

tH 


ON 

CO 

ro 

VO 

CM 

LO 

O 

ro 

10 

Hi 



tH 




ID 

CM 

ro 


tH 


V0 

O 

V0 

CM 

UN 

o 

• 

• 

» 

• 

• 

9 

ON 

ON 

l> 

CO 

ON 

ON 

ON 

ON 

ON 

ON 

ON 

ON 

H< 



VO 

CO 


VO 

tH 



tH 


• 

• 

00 

ft 

• 

UN 

CM 

o 

o 

H* 

!> 

o 

1 

H 

co 

tH 

I 

CM 


tH 

« 

1 


• 

ON 

CM 

VO 

ON 

H* 


H< 

CO 

r- 

l> 

CM 

CM 

• 

• 

« 

• 

9 

ft 


ro 

ON 


ro 

O 

ro 

ro 

00 

UN 

t> 

CM 

tH 

CM 

ro 


UN 

VO 


LO 

00 


CM 

CO 

VO 

• 

ft 

♦ 

ON 

ON 

CO 

ON 

ON 

ON 

VO 

CO 

OJ 

CO 

tH 

Hi 

• 


9 

H* 

CM 

CM 

LO 

VO 

tH 

o 

Mh 

I> 

CO 

ON 

tH 

’ ft 

• 

•ft 

H< 

UN 

CM 

CO 

CO 

ON 

CM 

VO 


VO 

VO 

l> 

• 

'• 

ft 

ON 

ON 

00 

ON 

ON 

ON 


VO 

H 1 

tH 

ON 

CM 

ON 

H 1 

VO 

• 

• 

• 


o 

VO 

CM 

ro 



00 

on as 


sh 

CO 


co un cm 

co co a\ 


o 


LO H 1 
O si’ 1 


ON ON j> 
ON ON ON 


VO CM CSI 

ID rl ON 

* • • 
VO ON ^ 

vH *H 


vH 

VO 

* 

O 

H< 


CM 

O co 

ft * 

O VO 
00 00 


r4 (SI CO ^ Lf) VO 


CQ 


60.64 18.16' 97.9 55.23 7.543 98.17 15.48 11.71 97.8 

76.17 3.409 99.41 70.65 .27 99.51 70.82 17.336 99.57 

24.35 3.39 99.65 29.77 . 6.787 99.7 29.51 25.62 99.87 
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TABLE 4.16 

Bus Voltages of the Sample System Shown in Fig. 4.2 for case b 

(conversion o i. double-circuit 3— phase lines to single circuit six- 

phase lines) 


Volt. 

cond. 

BUS 

No. 

TTonfig 

magnitude 

Magnitudes 
.1 ' 

Angle 

and Angles of Bus 
config if ‘ J ““ "" 
magnitude Angle 

Voltages 

config. Ill 
magnitude Angle 


1 

.9775 

0.0 

.8651 

0.0 

.81 

0.0 


2 

.9399 

11.57 

.8227 

10.66 

.7719 

12.43 


3 

.9481 

6.126 

.8166 

5.83 

.7742 

5.61 


4 

.9625 

-2.51 

.8484 

-3.22 

.7924 

-3.68 

A 

5 

• 9 1 9 6 

-4.42 

• 8191 

-4.815 

.762 

-5.32 


6 

. 955 3 

-3.625 

.8328 

-5.25 

.7778 

-6.16 


7 

.9418 

4.96 

.8269 

2.056 

.7778 

2.663 


8 

.9326 

1.4659 

.807 

-1.645 

.7713 

-1.39 


9 

.956 3 

2.97 

.8267 

1.6 

.785 

0.917 


1 

1 . 0444 

0.0 

1.0575 

0.0 

1.0592 

0.0 


2 

1.0278 

8.79 

1.0417 

10.75 

1.039 

10.6 


3 

1.0254 

3.92 

1.0402 

5.26 

1.0459 

5.4 


4 

1 .0303 

-2.197 

1.0435 

-2.143 

1.0452 

-2.1355 

B 

5 

1.0001 

-4.102 

1.0131 

-4.354 

1.0138 

-4.39 


6 

1,0107 

-4.7166 

1.0243 

-3.88 

1.028 

-3.75 


7 

1.0285 

3.2685 

1.0423 

5.37 

1.0396 

5.196 


8 

1.0177 

.1652 

1.0321 

2.042 

1.0284 

1.794 


9 

1.0328 

1.2286 

1.0474 

2.64 

1.0531 

2.81 
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20.91 2.08 99.67 12,904 2.976 99.69 ’11.54 2.195 29.7 
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Bus Voltages of the Sample System (Fig. 4.2) for Case b 

(conversion o£ double circuit 3-phase lines to single circuit 12- 

phase lines). 


Volt. 

cond. 

Bus 

NO. 


Magnitudes and Angles of 

Bus Voltages 


config. X 

Ha gn 1 tude Angle 

config. 

Magnitude 

II 

Angle 

config* 

Magnitude 

III 

Angle 


1 

.9162 

0.0 

.7391 

0.0 

.6705 

0.0 


2 

.8444 

15.02 

.6763 

12.56 

.6061 

16.9 


3 

.864 

8.22 

.6609 

8.257 

.6153 

6.43 


4 

. 9004 

-2.8706 

.7202 

-4.44 

.6505 

-5.42 

A 

5 

.8 395 

-4.95 

.6877 

-5.62 

.6152 

-5.95 


6 

.8969 

-3.5318 

.7088 

-6.19 

.6404 

-7.886 


7 

.8481 

6.83 

.6868 

-0.1064 

.6223 

1.244 


8 

.84 

2.57 

.6537 

-4.53 

.6223 

-3.38 


9 

.8733 

4.425 

.6748 

1.845 

.6319 

-0.908 


1 

1.0065 

o 

• 

o 

.9409 

0.0 

.9021 

0.0 

2 

.9801 

10.32 

.9083 

10.32 

.873 

11.16 

3 

.984 

5.237 

.9051 

5.24 

.8752 

5.34 

4 

.992 

-2.37 

. 9254 

-2.71 

.886 

-2.95 

B 5 

.9554 

-4 .24 

.8962 

-4.54 

.8558 

-4.95 

6 

.9814 

-3.88 

.9077 

-4.79 

.8594 

-5.25 

7 

.9814 

4.238 

.9107 

3.25 

.876 

3.5 

8 

.9718 

.9716 

.8956 

-0.2385 

.8674 

-0.28 

9 

.9918 

2.31 

.914 

1.79 

.8843 

1.71 


TABLE 4.19 
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Line Flows and Trans, Efficiencies of Sample System Shown in 

Fig. 4, 2 for Case b 

(conversion of: double circuit 3-phase lines to Single Circuit 12- 

phase lines) 


Volt. Line 



Line flows 

and efficiencies 



Cond.No. 

con fig. I 


config. II 


config. ±11 



P 

Q 

ri 

P 

Q 

r\ 

P 

Q 

V 


(MW) 

(MVAR) 

% 

(MW) 

(MVAR) 

% 

(MW) 

(MVAR) 

% 

1 

39.35 

53.28 

98.37 

15.04 

21.32 

98.25 

7.44 

22.21 

97.68 

2 

32.85 

30.32 

99.87 

56.55 

0.47 

99.63 

64.0 

2.67 

99.54 

A 3 

89.57 

14.95 

96.01 

111.56 

53.53 

98.73 

119.8 

32.32 

98.36 

4 

58.9 

34.84 

96.6 

35.029 

26.91 

96.15 

27.16 

14.77 

97.14 

5 

7 3.96 

1.76 

99.13 

51.39 

23.87 

98.81 

43.03 

6.202 

99.05 

6 

26.6 3 

18.133 

99.18 

49.92 

6.16 

98.63 

57.65 

7.42 

99.68 

1 

40.92 

29*74 

99.22 

34.46 

20.67 

99.33 

34.36 

21.1 

99.27 

2 

30.79 

6.05 

99.71 

37.06 

2.56 

99.56 

37.21 

1.802 

99.54 

3 

86.86 

9.757 

97.1 

92.48 

20.18 

98.2 

92.74 

14.0 

98.04 

B 4 

60.715 

22.14 

97.58 

54.57 

20.94 

97.43 

54.48 

15.39 

97.28 

5 

76.21 

.885 

99.32 

70.42 

6.72 

99.25 

70.22 

1.128 

99.22 

6 

24.36 

6.79 

99.56 

30.32 

4.79 

99.48 

30.46 

1.12 

99.72 







t 
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completely three-phase system, it moans that slack bus gene- 
rator should generate required values of reactive powers if 
possible or the generator should be replaced by another one 
if necessary to generate required power to keep the slack bus 
voltage constant for all conditions. The original generator 
if employed will limit the number of conversions/ however the 
reactive power of the generator will have to be adjusted. 

The comp axis ion of Tables 4.9/ 4.12 and 4.13 reveals 
that the bus voltages increase under voltage condition A and 
thy decrease under voltage condition B for six-phase uprating 
while they increase under both conditions for twelve-phase, with 
the increasing number of conversions of double-circuit three- 
phase lines to multi-phase lines for case a. 

When Table 4.9 is compared to the results of the 
load-flow study for case b in Tables 4.16 and 4.18, it is 
observed that the bus voltages decrease under voltage condition 
A for six-phase conversion and under voltage conditions A & B 
for twelve-phase conversion for the three types of system con- 
figurations. But under voltage condition B, the bus voltages 
increase for configuration I and II while for system confi- 
guration III, som e bus voltages have decreased in compari- 
sion to those for case 1L; however, they are higher than those 
of the original three-phase systems. 
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Prom Tables 4.10 and 4.13, the power flowing through 
lines 2 and 3 are (40.9+j22.8) and (30.69+jl.0328) in the 
original completely three-phase system while they are (39.2+j 
12 . 27 ), (31.94+J40.75) for configuration I, (25.85+j82.12) , 
(43.99+j 80.4) for configuration II and (26.26+jl43.23) , 
(46.3-Hjl38.15) for configuration III, and therefore, it nay be 
likely that original double circuit three-phase line may not 
be able for the required conversion for case a under voltage 
condition A. But the line flows in lines 2 and 3 have not 
changed much. 

If the two options for increasing the power transmi- 
ssion capability of a double -circuit three-phase line as 
mentioned in section 4,3 is considered, it is observed from 
equations (4,11) and (4.12), or (4.14) and (4.15) the equi- 
valent single-phase parameters for double circuit three-phase 
line with Increased transmission voltage will be the same as 
those of a single-circuit six-phase line if the leakage impe- 
dances of the transformers in each option is assumed to be 
the same ; since the transmission of more power by increasing 
the voltage of a double-circuit three-phase line will also 
require a three-phase/ three-phase transformer at the each end 
of the three-phase line similar to the up rating of a double- 
circuit three-phase line to a single-circuit six-phase line. 
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4.6 CONCLUSION 

The equivalent single-phase parameters for xr-circuit 
representation of a multi-phase line integrated with a three- 
phase wye/n-phase star connected transformer at each end of 
the line under different phase to ground voltages have been 
obtained to analyse the balanced network on a single-phase 
basis. The inpact of converting double -circuit three-phase 
to single— circuit multi— phase lines has been studied, it is 
observed ■flaw the results of the load flow investigation 
that, 

1. For the same slack bus reactive power, the bus voltages 
increase when the phase to ground' voltage of the single- 
circuit six— phase line is the same as that of the double— cir- 
cuit three— phase line while the bus voltages decrease when 
the phase to ground voltage of the line is f 3 times that of 
the three-phase line, and 

2. For the same slack bus voltage, the bus voltages decrease 
when the phase to ground voltage of the multi-phase line is 
the same as that of the double-circuit three-phase line while 
they increase when the phase to ground voltage is increased 
to y*3 times that of the three-phase line. 



CHAPTER 5 


TRANSIENT STABILITY ANALYSIS OP SIX:j?HASE SYSTEMS 
5.1 INTRODUCTION 

In the previous chapter, the load flow investi- 
gation has been carried out to study' the comparative perfor- 
mance of different configurations of" three-phaffi and six-phase 
transmission systems. This chapter presents the transient 
stability Investigation of the mixed three-phase and six-phase 
as well as comple tely six-phase systems in order to compare 
their performances with that of three-phase systems. Since, 
in the transient stability analysis, the model of a synchro* 
nous machine is required and the various models of three-phase 
machines, depending upon the accuracy desired in the analysis 
are available in the literature [55,69 ] . 

The section 5.2 starts with the detailed mathe- 
matical representation of six-phase synchronous machine; 
section 5 . 2.1 presents the direct and quadrature - axis transient 
and subtransient inductances, while section 5.2.2 deals with 
the different time constants for a six-phase synchronous 
machine. In section 5.3, a case study has been taken for 
transient stability Investigation based upon the sample network 
shown in Pig. 5.1, The system representation for this study 
has been discussed in section 5.3.2. In section 5.3.3, the 
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different system configurations of the sample network has 
been presented. The results obtained from the investigation 
have been presented and discussed in section 5.3.4. Finally, 
the chapter concludes in section 5.4. 

5,2 REPRESENTAT I° N 0F MULTI-PHASE (SIX-PHASE) 

Tho synchronous machine is represented by the 
different typos of model, depending upon the desiredi accuracy 
and the requirement of the analysis. These models require 
different typos of inductances, such as direct— and quadrature- 
axis transient and sub transient inductances, and also various 
types of time constants of the machine. The inductances and 
the time constants of a multi-phase (six-phase) machine have 
been calculated in the following subsections* 

5.2.1 Direct and Quadrature Axis Transient and Subtrarisient 

Inductances 

These inductances of a multi-phase (six- phase) 
synchronous machine have been derived with the help of the 
flux linkages and voltage equations in d q o component in 

equations (2,38) and (2.40). 

Since transformer voltages and are much 
smaller than the speed and and therefore, neglect- 

ing these voltages, the equations (2.38 ) and (2.40) in the 
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matrix form can be written as follows. 


y = Li - V* 3 

x q q q 

*!■ = V 3 I W :L a 

*G = - 1f3L mG i q 

*d “ V ' 3 l «d i a 

\ ^ 3 ^mQ^q 


L „i„ 
mF P 

L mG i G 


-f3 


L mD i D 

L mQ ± Q 


+ ^fV + L FD i D 
+ L GG. i G- " L GQ i d 
+ L PD 1 P + L DD 1 D 

’ l gq 1 g + l qq 1 q 


and 


v a 

ss 

- 

• r i , - 
a d 

Wm 

Xq 

v q 

22$ 

- r i + 
a q 

“fa 

v_ 

F 

XS 

- 



V G 

$25 

0 

0 

0 

U 

1 

11 

-"o 

V D 

5= 

0 

= - r D i D 

- pT d 

v a 

5S 

0 

= ~ Vq 

- P5 b 


Prom equations (5.1) and (5.8)# 'w e 9 et # 


V = - r i + »L a i. tf3« L i p + f 3“ L^i 

q a q a d 

?rom equations (5.3) and (5.5), we have 


El 


i D. 


, i_ 

*A 1 


VW"Vto +( r 3 L mD L FD" ^D L mF );L c 
-W/, pD +f D L FF + (f 3 HnP 1 ^^ 3 L FF L mD );L d 


(5.1) 

(5.2) 

(5.3) 

(5.4) 

(5.5) 

(5.6) 

(5.7) 

(5.8) 

(5.9) 

(5.10) 

(5.11) 

(5.12) 

(5.13) 


(5.14) 
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where a 1 » ^ L pp L DD -L pD ) 


Substituting i p and from equation (5.14) in (5.13), the 

following equation is obtained. 

V q r a i q +(J ^ L d“IT (L FF L rr0 + \p L DDT' 2 \p L m]7 D PD^ ^ i d 


uf 


F 


A 1 


(L 


i3 Of 


F^DD^mD^FD^ + A 1 ^ L mD L FF“ Ij mF L FD^ 


m. 


Prom equations (5.2) and (5.7), we have. 


v d - 


r a i d “ <d ^ L q i q “ ^ 3 L mG i G L mQ i Q^ 


(5.16) 


From equations (5.2) and (5.4) 


to -*» 


1 

A 2 


¥ G L QQ + \ h QQ + ^ 3L QQ L mG + ^ 3 WW 1 q 


¥ L +¥ L +(f 3L_ L ^+lf3 L Lw,)i 
G Gfl fl GG: GQ mG mQ GG J q 


(5.17) 


where a 2 « (L^L^-L^) 

Substituting i^ and i from equation (5.17) in equation (5.16), 


the following expression is obtained. 

s ** r a ± d" “ t h q~ a! (L QQ L mG +L mQ L GG +2L mG L GQ L mQ ) ^ x q 


v 


V3 w¥ V3w¥ _ 

. ’ *G , T T i T T firT, ._t, +L _L__) (5.18) 

+ ' A 2 L mG L QQ +L miQ G£T A 2 mG GQ mQ GG 


From equations (5.15) and (5.18), the subtransient inductances 


are 
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a 

ii 

■j* 

- TT 

(5.19) 

and 

cF 

u 

3 2 2 

l A' ,, 2 ll{L QQ L m G +L mQ L GG +2L mG L QQ L rnQ J ■ 

(5.20) 

Further* when the damper windings D and Q are not co ns ids red* 

we approach 

the transient state* and hence 

transient indue- 

ances are 

L d = 

2 

L d “ 3 L mP^ L PF 

(5.21) 

1 

and L ** 

q 

L q ~ 3 L mc/ L GG 

(5.22) 


where d and q refer to direct and quadrature axis quantities 

respectively. 


The expressions (5.21) and (5.22) give the direct and quadra- 
ture axis transient inductance of the machine; while expressions 
(5.19) and (5.20), the direct and quadrature axis subtransient 
inductances of the six-phase machine. 

5,2.2 Different Time Constants of the Six-Phase Machine 

The different time constants of the six-phase 
synchronous machine have been obtained in this section. 

These time constants are required for the different types of 

synchronous machine model. 
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Direct- axis transient open-circuit time constat T* : 

do 

In this case, a sudden voltage is applied to 
the field winding keeping armature open circuited without 
damper winding. 


Prom equations (5.3) and (5.9), 

voltage equation become 


L FF P *‘F +r F ;1 ‘F ^ F 


Thus, 


T „ « 

do 


I»pp/ 3 


F 


the field 


(5.23) 

(5.24) 


Direct-axis transient short-circuit time constant Tj : 

In this case the condition is similar to above 
except that the stator windings are short-circuited. 


From equations (5.3) and (5$, 

ft + + Vr = - v f 

From equations (5.8) and (5.1), we have, 

L,i, + f 3 L _ i„ = 0 
d d mF F 

Substituting values of Pi d from (5.26) in (5.25), 
[ L pp -( V 3 L mp ) 2 /L d ] Pi p + r p i p = - v p 

From equation (5.27 ) , 

. V ^-( f = T ' 

T , * 7 — C" h A do 


(5.25) 


(5.26) 


(5.27) 


(5.28) 
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Similarly# the quadrature— «axis opsn circuit transient time 

constants are given by 





(5.29) 


Direct Axis Subtransient open-circuit time constant T, s 

do 

The condition is similar to above except for the 

presence of damper windings. 


From equations (5.3) and (5.9)# we have* 

L FF P ip + L FD Pi D + r F i F = “ V F (5.30) 

Similarly, from equations (5.5) and (5.11) 

LpjjPip = 0 

At t *0 + # IL* 0# hence from equation (5.5)* 

^FD^F + '^DD'So ~ ® (5.32) 


Substituting ip from equation! (5.32) in equation (5.30), and 
Pip from equation (5.30) in (5.31)* the following differential 

equation is obtained. 

2 2 

PiiQ +C(r D L pF +rpL DD )/(L DI) Lp F -Lp D )] i D = - v F L Fl/ ^D^FF^FD^ 

(5.33) 


and from equation (5.33)# 

2 

M ^DD^F F ** L FD 

"^dO r D^FF + 


(5.34) 
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Direct Axis short circuit subtransient time constant • 

From equations (5;9) and (5,11) after substi- 
tuting the values of and ¥ from equations (5.3) and 
(5.5), the field and the damper d-axis winding voltages are 

^ ■ I imF P ^d + L FF P ^F + L FD P '*D + r F^F = ” V p (5.35) 

^" 3 L mD Pi d + L FD Pi F + L DD Pi D + Vd ~ 0 (5.36) 

At times t « o + , '? F =0 as the voltage is suddenly applied. 


Neglecting the stator winding resistance from the equation 
(5.8), and assuming = 0 since stator windings are short 

circuited, we get from equations (5.1) and (5.3), 


L d i d + L mF i F + HnD^) 0 

^ 3 L mF i d + L FF 1 F + L FD ± D = ° 


From equations (5.37) and (5.38), 


’V 

d, 

1 

L FF^ 3 L mD^ 3 L mF L FD 

ip 

(Lp p L d -V 3L mp V 3 ^ nE i5 r ” 

~^ 3 L mF^ 3L mD + L FD L d 


From equation (5.36) 


Pi- « - rp-[f3 I^ 0 Pi d +Ita D pl i3 +r E ;L D ^ 

t U FD 


(5.37) 

(5.38) 


(5.39) 


(5.40) 
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Substituting the equation (5.40) in the equation (5.35) and. 
the value of ip from (5.39), the following differential equa- 
tion is obtained. 

[ (L FF L DD~ L FD > (L FF L d ~^" 3 L mP ) + (L PF^ 3 L mD^ 3 L mP L FD J 
-LpplfS L^) ] Pi D + [ r D L FF (L Fp L d -(\T3L mF ) ) “ r F L FD ^ 3L mF^3L rnD 

-Lp p L d ) 3 i D * V F L FD ( L FF L d ~ ^ 3 L mF^ ( 5 . 4 i) 

From equation ( 5 . 41 ), . 

T " „ L FF^ L d ( '^FF I t)b‘‘* L FD^ ~ L FF ^ 3 L mD^ + 2 ^ 3 L rnD^ 3 L mF L FD~ ^ 3 L mF* ^ 

r D L FF (L FF L d ” ( ^ 3 L mF^ * ” r F L FD ( ^ 3 L mF^ 3L mD ' L FF L d* 

(5.42) 


Assuming r p bo be negligible in equation (5.34), 


H 

« * T r 

T d _ L d " 


T 


do 


(5.43) 


Similarly, the quadrature-axis open circuit and short-circuit 
subtransient time constants can be obtained which are as follows 


T 


i i 


qo 


l QQ l gs ~ L op- 

r Q L GG + r G L QQ 


and. 


I* 


T 


»* 


L a " 
T ,o 

q 


(5.44) 


(5.45) 
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Further, til© voltage behind the stou transient reactances and. 


their rate of change have been obtained as follows. 


From equation (5.15) , 


E - The voltage behind q-axis subtransient reactance 

q 




Af’’ L mF L DD”^ 3 L mD L FD* + ' 


Al " ^ 3L mD L FF" ^ 3L mF L FD^ 


(5.46) 


Differentiating equation (5.46) w.r.t. time, and substituting 


P X J' F and P V I' D from equations (5.9) and (5.11), we get. 


^ 3wIj mP ^ 1 " L rnD L FD //L mF L DD^ ^r p + V ^ 3WL mD i D^ 1 ~ L mF'4’D //L inD L FF^ 
PE'* = >5 ! + T 


LfF L DD“ L FD 


l ff l, dd" l fd 


r F L DD 


r I *v i t 

D FF 


(5.47) 


Substituting for - L FD from equation (5.34) in (5.47) , 


we get 


PE * 1 

q 


(E +E >K1 


K2 E. 


(5.48) 


*» x do< 

T 1 . (1 + ) 

d ° T-. 

D 


M x n 

T ac < 1+ V-> 

■ L do 


where. 


^ * 1 " L mD L FD / ' L mF L DD 


K2 * l“ L mF L FD^ L mD L FF 


T D “ W r D 
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Similarly, from equation (5.18), 


E 


I * 


V 30% If 30^ 

>2 — (jj t, +h l ) + Q. 

-2 mGGQ rnQ! GQ ; 


L GG L QQ”"’ I i3Q 


T T ^ L mG L GQ +L mQ L GG^ 

L gg l qq- l gq 


(5.49) 


Differentiating the above equation w.r.t time and substituting 
the values of and from (5.10) and (5.12), we have. 


0 


pE ^ ^ 3 ML mG [ 1 +L mQ L QG //L mG Ll aQ \ I G J 3(JL raQ [ 1 +L mG L G£/ L mQ L GG ] I G 

” J 1Tr vr 1,T,r 111 ' 2 v /„ " ' , _ _ I 2™ 


M 


( Xj cx} Ij OQ“ L GQ^ //L QQ r G 


^ l qq l gg“ l gq ,/:c q l gg 


(5.50) 


Introducing the equation (5.44) in the equation (5.50), we 

have 


PE 


r« 


K3 E, 


T 


T 


1 (1+ -B®> 


qo 


Q 


K4 E 




t (l + — r^-) 

qo m 

qo 


(5.51) 


where, 

K 3 * 1 + Ij mQ L GQ // ' L mG Ij QQ 
K4 » 1 + ] 4nG L G</ L mQ L GGJ 

T Q * L GQ /r Q 

Equations (5.48) and (5.51) give the rate of change of voltage 
behind q and d-axis subtransient reactances. 



175 


when the d— and q-axis damper windings are not considered, 
the transient state is approached, and the equations (5.47) 
and (5.50) reduce to 


PE 31 

q 




E ) 

q 


(5.52) 


and 


PE d " 



(5.53) 


5.3 A CASE STUDY 

The transient stability investigation has been 
carried out on a sample system as shown in Fig. 5.1 for the 
different system configurations in order to compare the per- 
formance of a mixed three-phase and six- phase system with 
that of a completely three-phase system. Also, a completely 
siX-phase sys tern has been compared with the conventional 
three-phase system, 

5.3.1 Sample Network 

A nine bus three-machine sanple network [ 55 ] 
has been taken for transient stability investigation. The 
line and the transformer data are provided in Appendix B : j 
the synchronous machine data in Table 5.1. Prefault condi- 
tions for the different configurations with the same loading 
conditions are given in chapter 4. The load flow data for 



m 


six-phase line 


FIG. 5.1 Nine-bus three machine sample system 
for the transient stability Investigation 
Line nos. 2,3 and 6 are six phase lines 
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the increased loading has bean shown in Table 5*2* The 
pre-fault voltage for the Increased loading for a mixed system, 
and for a completely six-phase system with the same loading 
and with increased loading have been provided in Table 5.3, 

5.3.2 System Representation 
1) Machine Representation : . 

Synchronous machines have been represented by the classical 
model whore ‘the changes in flux linkages and the saliency 
are neglected. Thus, the synchronous machine is represented 
by a constant voltage source behind the direct-axis tran- 
sient reactance for the purpose of this study. 

ii) Network Representation : 

The network performance equations employing single-phase 
equivalent representation of various elements are described 

as 

YV = I (5.54) 

The solution of equation (5.54) is obtained by Gauss- 
Siadol iterative technique. 

iii) Load Representation: 

The load is modelled as a constant impedances from bus to 

ground. 
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iv) Swing equation : 

The equations for electro-mechanical motion of the machine 
employing the simplified model are 

dd, 

dt * W i (t> ~ 2nf (5.55) 


d0) i = jk£ 
dt 


- PJt)) 


(5.56) 


where, i*l, . . . ,m; m = 3, the number of machines in operation 
the effect of damping and governor action has been neglected 
in this model. 


5.3.3 Sample Network with Different Configurations : 

The different configurations under transient 
stability investigations are : 1) a completely three- phase 
system. Line No. 2,3 and 6 have been assumed as double-circuit 
three-phase lines; 2) a mixed three- phase and six-phase 
systems where the six-phase lines are obtained by converting 
double— circuit three-phase lines; thiee cases have been considered, 
for example t a) when only line no. 2 is converted, b) when two 
lines (nos. 2 and 3) are converted and c) when all the double- 
circuit three-phase lines (nos. 2,3, and 6) are converted ; 

3) a completely six-phase system where the single-circuit as 
well as the double circuit three-phase lines are converted to 
single circuit six-phase lines. 
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When a three-phase line is converted to a six- 
phase line/ two voltages have been considered for the study 
of mixed three-phase and six-phase systems with the same 
loading as that of the completely three-phase systemsj in one 
case the line to ground voltage of the six-phase line is f 3 
times that of the double circuit three-phase line, and in the 
second case, the line to ground voltage of the two lines (i.e., 
three-phase and the six-phase) is assumed to be equal. If the 
power system is assumed to meet the increased demand/ say by V" 3/ 
then with this increased system loading/ the line to line 
voltages of the two systems have been assumed to be equal. 

For a completely six-phase system, the transient stability 
has been Investigated with the same loading as the original 
three-phase system, and also with the increased loading by 
assuming the same phase to ground voltage of the two systems. 
The system has been assumed to be balanced and therefore, the 
study has been conducted based on one line diagram. 

As the conversion requires one three-phase/ six- 
phase transformer at each end of the six-phase line, the 
single-phase equivalent of a six-phase line integrated with 
two transformers has been calculated with 8% leakage impedance 
of the transformer from equation (4.12) and (4.13). 



180 


When the system becomes completely six-phase, it 
requires six— phase transformers in place of three-phase trans- 
formers, in this case, the leakage impedances of the six- 
phase transformers have been assumed to be the same as those 
of three-phase transformers. The six-phase synchronous 
machines have been assumed to have the same induc tan ces (i.e. 
inductive coefficients) as those of three-phase machines. 


TABLE 5.1 

Machine Data on 100 MVA Base of the Sample System 

Shown in Fig. 5.1 


Machine 

NO. 

connecting 

Bus 

Inertia constant 
(Sec.) 

Direct-Axis Trans, 
reactance 

1 

1 

23.64 

(1) 

.0608 




(3) 

.1216 

2 

2 

6.4 

(1) 

.1198 




( 3 ) 

.2396 

3 

3 

3.01 

( 1 ) 

.1813 




(3) 

.3626 

In the last column 

(1) and (3) refer to 

the 

direct-axis 


transient reactance for three-phase and 6 -phase machines 

respectively. 
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TABLE 5.2 

Assumed Bus Voltages for the Load Flow Analysis of the 
Sample System Shown in Fig. 5.1 for increased Loading 


Bus Assumed bus Generation Load 


No. 

voltage 

MW 

MVAJR 

MW 

MVAR 

1 

l,04+j0.0 

0.0 

0.0 

0.0 

0.0 

2 

l.OO+jO.O 

163.0 

6.7 

0.0 

0.0 

3 

l.OO+jO.O 

85.0 

-10.9 

0.0 

0.0 

4 

l.OO+jO.O 

0.0 

0.0 

0.0 

0.0 

5 

l.OO+jO.O 

0.0 

0.0 

216,5 

86.6 

6 

l.OO+jO.O 

0.0 

0.0 

156.0 

52.0 

7 

l.OO+jO.O 

0.0 

0.0 

0.0 

0,0 

8 

l.OO+jO.O 

o 

♦ 

O 

0.0 

173.0 

60.6 


l.OO+jO.O 


9 


0.0 


0.0 


0.0 


0.0 
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TABLE 5.3 


Bus Voltages of the Sartple System Shown in Pig. 5.1 


Bus _ B us Voltage _______ 

No. (T) when line nos 2,3, & (2) for a completely six- 

6 are converted to phase 

six— phase lines with same with increased 

lo ading loading 


1 

1.04 

0.0 

1 . 04 

0.0 

1.04 0.0 

2 

1.078 

-8.174 

1.1563 

4.15 

1.0107 -4.223 

3 

1.0746 

-11.605 

1.1582 

1.6943 

1.0198 -6.6475 

4 

1 . 0464 

-9.212 

1.0698 

-1.0496 

1.0237 -4.7036 

5 

1.0369 

-15.863 

1.0768 

-2.3654 

.9905 -9.192 

6 

1.0434 

-13.867 

1.1183 

-2.1237 

1.0136 -.30.225 

7 

1.0783 

-13.2016 

1.1553 

1.9645 

1.0099 —7.0835 

8 

1.0687 

-16.471 

1.1536 

.4756 

.9975 -9.87 

9 

1.0815 

-14.0615 

1.1612 

.6331 

1.0232 -8.0152 


Slack bus generation (I) = (302.47, 12.74) 

Slack bus generation (Il) = (70.76, -107.0) 

with same loading 

Slack bus generation with increased loading = (303.14,71.16} 
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5.3.4 Result, and Discussions 

A s y mine tx leal (n— phase to ground) fault is 
assumed to occur at the buses 4, 7 and 9, one at a time for 
the different power system configurations under various post- 
fault conditions. Graphs (Machine angle difference vs. time) 
have been plotted as shown in Figs. 5.2 - 5.7 for 3 cases : 

I) a completely three-phase system II) a mixed three-phase and 
six-phase system under voltage conditions A; III) a complete- 
ly six-phase system under voltage condition B; the same system 
loading has been considered. The machine angle difference vs. 
time for the fault at the buses 4,7 and 9 are observed to 
compare the performance affected by the conversion of double- 
circuit three-phase lines to six-phase lines, and also to 
study the impact of converting a three-phase system to a comp- 
letely six-phase system. It is observed from the graphs that 
there is irrpor o vemen t in the stability with the six-phase 
system. 

The critical clearing times for the faults at 
the buses have been shown in Tables 5. 4. -5. 11. It is observed 
from the tables that when a double circuit three-phase line 
is converted to a single circuit six-phase line with the same 
line to line (adjacent phase) voltage,, the critical clearing 
time increases as the number of conversions increases 
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FIG. 5-6 Plot of ( 62 - 5 ,) Vs time for the fault near bus 
The fault is cleared in 0-14 Sec. by opening line 




near bus 
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irrespective of the fauJlt locations and the post-fault con- 
figurations, But when the power is transmitted at the sane 
phase to ground voltage as that of the three-phase line, it 
is observed from Table 5,5 that, when line no. 2 is converted 
to six-phase line, the critical clearing time for different 
post— fault conditions are the same. For example, when the 
fault near the bus 4 is cleared by tripping line 2, which 
is a six— phase line, the critical clearing times for the two 
configurations are ,24 sec. - ,25 sec. 

When line nos. 3 and 6 are converted to six- 
phase lines, it is observed that the critical clearing time 
is less when the line no.l is tripped, and the critical 
clearing time under configurations II and III are the same, 
i.e., ,21 _ ,22 Sec. It neans that the conversion, of line 
no. 3 reduces the stability and line no. 6 does not affect it 
when the fault is at bus no. 4, Therefore, it can be concluded 
that, whether the conversion to six-phase system reduces the 
stability or not depends upon the location of that line. From 
Tables 5.7 and 5.9, which provide the clearing times for the 
fault at buses 7 and 9 respectively, the critical clearing 
time increases and decreases as the power is transmitted at 
the same adjacent phase voltages ox at the same phase to 
ground voltage for the same loading. Further, the critical 
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TABLE 5.4 

f ° r ^ ^ r6e -Phase to Ground 
Fault near the different Buses of a Completely Three. 

Phase Sample Network 


Faulted bus Post fault condition 


cct (in Sec.) 


a) Fault removed 

b) Line 1 is tripped 

c) Line 2 is tripped 


.25 - .26 
.22 - .23 
.24 _ .25 


a) Fault removed 

b) Line No. 3 is tripped 

c) Line No. 5 is tripped 


.16. - .17 
.09 - .10 
.14 - .15 


a) Fault removed 

b) Line No. 4 is tripped 

c) Line No. 6 is tripped 


.19 - .20 
.16 - .17 
.17 - .18 
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TABLE 5*5 


a /iT f ° r the n-Phase to Ground fault 
at the bus No. 4 of the Sample System for C ase a 


System Post fault conditions 
con fig. 


cct (in Sec. ) 
V oltage condition 
A B 


I 


a) 

fault removed 

.26 - .27 

.25 - .26 

b) 

line 1 is tripped 

.23 - .24 

.22 r .23 

c) 

line 2 is tripped 

.24 _ .25 

.24 - .25 


XI 


a) 

fault removed 

.29 - .30 

.23 - ..24 

b) 

line 1 is tripped 

. .27 - .28 

.20 - .21 

c) 

line 2 is tripped 

.29 - .30 

.21 - .22 


III 


a) 

fault removed 

.34 - .35 

.22 - .23 

b) 

line 1 is tripped 

.32 - .33 

.19 - ,20 

C) 

line 2 is tripped 

.33 - . 34 

o 21 - .22 


Voltage configuration! : 

A refers to the mixed system when line to line (i.e 

adjacent line) voltage of six-phase line and the double 
circuit three-phase lines are the same, 

B refers to the mixed system when the six— phase line and 
the double circuit three-phase lines have the same 
jftase to ground voltage. 


193 


TABIjE 5,6 


«»a- to Ground 
*» WO. 4 of the Sample System for Case b 


System 

config, Post fault condition 


cct(in Sec.) 
Voltage condition 


a) fault removed 

b) line 1 is tripped 

c) line 2 is tripped 

a) fault removed 

b) line 1 is tripped 

c) line 2 is tripped 

a) fault removed 


.21 

- .22 

.25 

- .26 

.19 ; 

- .20 

.22 

- .23 

.20 . 

- .21 

.24 

- .25 

*18 - ,19 

.24 

- .25 

.14 - .15 

.21 

_ .22 

.17 - 

• *38 

.23 

- .24 

.16 _ ,17 

-24 

- .25 

.12 - 

■ *13 

.21 

- .22 

.15 - 

.16 

.23 

- .24 
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TABLE 5.7 


Critical Clearing Time for the n-Phase to Ground Fault 
at Bus No. 7 of the Sample System for Case a 


config. P° s 't fault condition cct (in Sec.) 

Voltage condition 
A B 


I 

a) 

fault removed 

.17 - .18 

.16 - .17 


b) 

line 3 is tripped 

.12 - .13 

.07 - .08 


c) 

line 5 is tripped 

.15 - .16 

.14 - .15 

II 

a) 

fault removed 

.20 - .21 

.15 - .16 


b) 

line 3 is tripped 

.16 - .17 

.05 - .06 


c) 

line 5 is tripped 

.19 - .20 

.09 - .10 


III 

a) 

fault removed 

.22 - .23 

.15 - .16 

t 


b) 

line 3 is tripped 

.19 - .20 

0.0 - 0.0 


c) 

line 5 is tripped 

.21 - .22 

.08 - .09 
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TABU3 5.8 


Critical Clearing Time for the n —Phase to Ground Fault 
at Bus No. 7 of the Sample System for Case b 


con fig. Post fault condition 


cat (in Sec.) 
)ltage ^condition 
A B 


a) fault removed 


.15 - .16 .17 - .18 
.05 - .06 .07 - .08 


b) line 3 is tripped .05 - .06 .07 - .08 

c) line No. 5 is tripped .11 - .12 .14 - .1$ 


a) fault removed 


.14 - .15 


b) line no. 3 is tripped 0.0 - 


.16 - .17 


.08 - .09 


c) line no. 5 is tripped .12 - .13 .11 - .12 


a) fault removed 


.12 - ,13 .16 - .17 


b) line no. 3 is tripped 0.0 

c) line no. 5 is tripped .09 - .10 


.06 - .07 


11 


12 
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TABLE 5.9 


Critical Clearing Time for the n— Phase to Ground. Fault 
at Bus No. 9 of the Sajple System for Case a. 


config. Post fault “nation 


cct (in Sec.) 
Vol tage condi tion 
X b“ 


I a) fault rem oved .19 - .20 .19 - ,2o 

b) line 4 is tripped .17 - *18 .16 - .17 

c) line 6 is tripped .18 - .19 *17 - .18 


II a) fault removed *21 - .22 .18 - .19 

b) line 4 is tripped .20 - .21 .13 - .14 

c) line 6 is tripped .20 - .21 .17 - .18 


III a) fault removed 

b) line 4 is tripped 

c) line 6 is tripped 


.22 - .23 .18 - .19 
.22 - .23 .11 - .12 
.21 - .22 .17 - .18 
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TABLE 5.10 


Critical Clearing Time for the n-Phase to Ground Fault 
a Bus No. 9 of the Sample System for Case b 


System 

config. 


Post fault condition 


cct (in Sec.) for 
different voltage condi. 


A 


B 


I 


a) fault removed 

.18 - .19 

.19 - .20 

b) line no. 4 is tripped 

.13 - .14 

' .16 - .17 

c) line 6 is tripped 

,16 - .17 

.17 - .18 


XI 


a) 

fault removed 


.16 - .17 

.19 - .20 

b) 

line no. 4 is 

tripped 

.11 - .12 

.15 - .16 

c) 

line no. 6 is 

tripped 

.14 - .15 

.17 - .18 


III 

a) 

fault removed 


.15 - .16 


b) 

line no. 4 is 

tripped 

• 

0 

1 

» 

h* 


c) 

line no. 6 is 

tripped 

.13 - .14 


.18 - .19 
.13 - .14 
.17 - .18 
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TABLE 5.11 

Critical Clearing Time for n-Phase to Ground Fault with 
Increased System Loading for Case a 


Faulted Post fault condition 

Bus 


cct(in Sec.) 

System configuration 

when line For a completely 

nos , 2 , 3&6 6 -phase system 

are 6— phase (volt, config. B) 

lines same increased 

(Volt, config.) loading loading 


A 

1) 

fault removed .77 - .78 

.33 - .34 

1.54-1.55 

2) 

line no.l is tripped .47 - .48 

.32 - .33 

1.19-1.2 

3) 

line no. 2 is tripped ,60 - .61 

.32 - .33 

1.24-1.25 

1) 

fault removed .27 - .28 

.20 — . 21 

.26-. 27 

2) 

line no. 3 is tripped .26 - .27 

.18 - .19 

.24-. 2 5 

3) 

line no. 5 is tripped .28 - .29 

.18 - .19 

.24-. 25 


1) fault removed * 28 - * 29 

2) line no. 4 is tripped .27 - .28 
line no. 6 is tripped .27 - .28 


.21 - .22 .26-. 27 
.19 - .20 .25-. 26 
.20 - .21 .25-. 26 
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clearing times for the fault at the different buses with the 
help of prefault conditions obtained for case b are provided in 
Tables 5*6, 5.8 and 5.10, it is observed from the Tables 
that the critical clearing times are less under both trans- 
mission voltages/ however , the critical clearing times are 
more in case when the six— phase lines have the same phase to 
ground voltage as the three-phase line. 

In Table 5.11, the critical clearing times for 
the increased loading i.o, , ^"3 times that of the completely 
three— phase system, have been tabulated for the different 
cases : 1) when linos 2, 3 and 6 are converted to six-phase 
lines and the power is transmitted at the same line to line 
(adjacent phi so) voltage of the two systems and (2) a comple- 
tely six-phase system; case 2 has been considered for two 
loading conditions with phase to ground voltage equal to that 
of three-phase system. 

It is observed from the comparision of Tables 
5.5, 5.7, 5.9 and 5.11 that the critical clearing times for 
the mixed sys tom with increased loading are higher with those 
having the same loading. For a completely six-phase system 
with the same loading as well as with the increased loading, 
the critical clearing time for the two loadings are more than 
that for Idle completely three-phase system. It is further 
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observed that the critical clearing times for a completely 
six-phase system are higher with increased loading than with 
the same loading as the corrpletely three-phase system. 

5.4 CONCLUSION 

The direct- and the quadrature- axis transient 
and sub transient inductances as well as the different time 
constants of a six-phase synchronous machine have been cal- 
culated which will be useful in analyzing the six— phase 
machine as well as the power system network wher ever it is 
incorporated. The transient stability has been investigated 
for the different power system configurations epecially with 
a view to study the impact of converting double circuit 
three-phase lines to single circuit six-phase lines. The 
results reveal the following observations. 

1. If the phase voltage i.e. line to neutral voltage of the 
six-phase line is equal to V" 3 times that of the double 
circuit three-phase line, the stability increases with the 
same loading as well as with increased loading. 

2. if the uprating of the double-circuit three-phase lines 
are carried out on the same loading with equal phase 
(i.e. line to neutral) voltage, the stability reduces. 
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3. If a three-phase system is converted to a completely six- 
phaae system , the stability of the six-phase system 
increases with the »ame loading as well as with increased 

loading, when the phase voltages of the two systems are 

equal. 

4. The impact of converting a double-circuit three-phase 
lino to a single— circuit six-phase line on the stability 
depends on the location of the converted line from the 

fault. 

5. If the power is transmitted through six-phase line having 
the adjacent voltages equal to those of double-circuit 
three-phase lines, the stability of the mixed system is 
more with increased loading than with the same lo adi ng 

as that of the completely three-phase system, and also, 
a completely six-phase system is more stable with increased 
loading than with the same loading when both the systems 
have the same phase to ground voltage. 



CHAPTER 6 


DYNAMIC STABILITY ANALYSIS OF SIX-PHASE SYSTEMS 

6.1 INTRODUCTION 

In the previous chapter the transient stability 
investigation has been carried, out on a sanple network for 
its different configurations in order to study the impact of 
converting double circuit three-phase lines to single-circuit 
six— phase lines on the transient stability. This chapter 
presents an elaborate study of dynamic stability of the 
multi-phase (six-phase) system. For the dynamic stability 
study, various mathematical models of three-phase synchronous 
machines arc available [ 55,56,69]. A mathematical model of 
the multi-machine power network has been provided in section 

6.2 for the present study. A linearized model of the m-port 
reduced network has been given in section 6.2.1, and in 
section 6*2.2 a linearized model of a power system network 
has been provided incorporating the machine details. A 
case study for dynamic stability investigation is conducted 
in section 6.3. A sample system for the purpose of study 
has been shown in section 6.3.1* The various system confi- 
gurations of the sample network have been discussed in 
section 6.3.2. Further, the results of the case study have been 
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presented and discussed in section 6.3.3. Finally, this 
chapter concludes in section 6.4 based on the discussions of 

the results of investigations. 

6.2 DYNAMIC MODEL OF THE MULTIMACHINE POWER SYSTEM 

Various models, such as two-axis model, one-axis 
model, classical model etc, for dynamic stability study of 
multi-machine system are available in the li ter ature[ 55,70,71], 
The consideration of a particular model for the study depends 
upon the availability of data and also to what extent the 
accuracy is desired in the analysis. Here, for the present 
study, the two-axis model and the classical model for a syn- 
chronous machine has been considered. The voltage and the 
current relations of machines are presented. The linearized 
model for the power system network has been given. 

6.2.1 Linearized Model of Power System Network 

The linearized model of a reduced power system 
network (Fig. 6.1) can be expressed by a set of equations 

as given below (Appendix C). 
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( 6 . 1 ) 


6.2.2 Linearized Models of Synchronous Machines 

The two models for synchronous machines are presen. 

ted. 

1. Twrw,xi' k Mo del Of a S ynchrpn ou£jMachine 

in the two axis model, the transients effects are 
while the sub transient effects are 


taken into consideration 
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neglected [55 J . The transient effects are dominated by the 
rotor circuit. The rotor has two circuits; one main field 
winding on its d—axi s and the other equivalent q— axis formed 

by the solid rotor* 


The two -axis model of each machine is given by 


Vi ^di 88 " E di - < X qi " V V 


( 6 . 2 ) 


t' . e' = E wn . -E . -1- (x,. - x' .) I,, 
do! qi FD1 qi dx dx dx 


(6.3) 


M *) = T , - ( I , . E* . + I . E* . ) - D. o). (6,4) 

“i i mi v di di qi qi ’ x x 


i i a 


(6.5) 


where 1 * 1,2, m is the number of machines. 

A linear model for synchronous machine has been 
obtained by linearizing equations (6.2)-(6.5), which are as 

follows. 

(6.6) 


qoi A ^df ” (x qi - V AI qi" AEdi 


* I 


T AE = AE„. -AE . + (x,. " x di^ Al di (6»7) 

l do± qi FDx qi dl ai 


M< A*S , = AT. - of i -«dio AE di * E di° iIdl 


' l ± " w i ~ mi 


, B ‘ a I . + I a io AE qi^ 

+ a qio qi qi ° H 


( 6 . 8 ) 
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A A, =aW 

°i i 


(6.9) 


It is observed from equations (6.6) — (6.9) that 
each machine represents a fourth order system with state 
variables a E^^a E^, a^K and a6 ^ ^but the unknowns are six, 
namely aI^ and aI^ in addition to the state variables.- 
These two unknown parameters of each machine can be found 
from equation (6.1). 


In equation (6.1), a is increments in internal 


generator voltages, 


2. A Classical Model of a Synchronous Machine 


The classical model of the synchronous machine 


is given by 


' S i - T mi - E i V 'Vi 


( 6 . 10 ) 


6 i = W i ” W o (6.11) 

Thus, the linearized classical model of the synchronous machine 
from equations (6.10) and (6.11) are 


H ±A W i » AV - E io Al qi - Djl AW i 
A ® j _ 


( 6 . 12 ) 


(6.13) 
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where i * 1.2,,.., r; r is the number of machines represented 

by classical model. 


Thus# the machine represented by the classical model will 
have only and a 6 r as state variables, in equation (6.12). 

E io known and AI^.^ f s unknown variable that is eliminated 

with the help of equation (6.1). 


If the machines in the power system network are 

represented by both two-axis model and the classical model, 

then in order to eliminate aI of the classically represented 

qr 

machine, a relation between aI and the currents of the 

qr 

machines represented by two-axis model is obtained. Also, the 
vol tages aV^ andAV^ of the machines represented by two-axis 

model are obtained ih terms of aI . and aI .. of these machines, and 

qi di 

voltages E r of the machines represented classically [55] . 

The dynamical modelling for the three-machine 
power system network has been given below. 


A Particular Case : 

A linear model of a power system network with 
three synchronous machines has been developed. 

Taking linear equations (6.12) - (6.13) for a 
classical model of machine 1 and equations (6.6) - (6.9) for 
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a two-axis model of machines 2 and 3, the resulting equations 

in the matrix form will be 
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( 6 . 14 ) 


From equation (6*1) with V replaced by E and for a three' 
machine system (using 6 2 i = '“ <5 i2 / 6 31 = ~ 6 13^ # we have 
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where t 

X = f A« n A E ' A E ' A« A E ' . A E ' A« A A S. ., ] 
-to L 1 q2 d2 2 q3 d3 3 12 13 J 

. t 

a«n = k T ml A E PD2 a T m2 A E FD3 A T m3 ] 

The system matrix [ A^j is obtained by substituting the current 
variables from equation (6.16) in (6.14), and by combining 


the first two 

terms 

of equation 

(6.14) 

which 

then 

gives 

the 

coefficient matrix 

[*J o£ 

the , 

state variables X 

— m 

. The 


matrix 1 

:A J is 

given 

below 








a n 

a l2 

a i3 

0 

a 15 

a 16 

0 

a l8 

*• — t 

a 19 


0 

a 22 

a 2 3 

0 

a 25 

a 26 

0 

a 28 

a 29 


0 

a 32 

a 33 

0 

a 35 

a 36 

0 

a 38 

S 39 


0 

a 42 

a 43 

a 44 

a 45 

a 46 

0 

a 48 

a 49 


0 

a 52 

a 53 

0 

a 55 

a 56 

0 

a 58 

a 59 


0 

a 62 

a 63 

0 

a 65 

a 66 

0 

a 68 

a 69 

: 

0 

a ?2 

a 73 

0 

a 75 

a 76 

a 77 

a 78 

a 79 


1 

0 

0 -1 

0 

0 

0 

0 

0 


1 

0 

0 

0 

0 

0 

-1 

0 

0 


(6.18) 
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where 

D 1 

a ll = ~ 

E 1 

a l2 ~ - y ll * a 13 

E 1 

■ a 15 * hJ’ y l3' a X6 

E X 

a 18 “ “ y l5 

E 1 

a l9 M 1 F 16 

a 22 = [ -1 + (x d2 - x 2 ) B 22 ] 

do2 

a 23 = --T [ (x d2 - x^) G 22 3 

(io 2 

a 25 85 “T 1 ^ x d2 x 2* y 43 

• T do2 

a 26 = ~7 < X d2 "* *2* y 33 

do2 

a 28 = (x d2 “ x 2^ y 45 

i do2 

a 2 9 = -4~ <* d2 “ x 2> y 46 

a o2 


a 12 

“ ” a 15 
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a 32 

a 33 

a 35 

a 36 

a 38 

a 39 

a 42 

a 43 

a 44 

a 45 

a 46 

a 48 


1 

' (x q2 

“ x 2^ G 22 

T qo2 

X 

' E 1 + 

(x q2 * X 2^ B 22 3 

- t 

T 0 
qo2 

1 

<X q2 

- X 2> y 33 

T qo2 

1 

<X q2 

“ X 2^ y 34 

T qo2 

1 

Cx q2 

- x 2> y 35 

T qo2 

1 

<x q2 

1 

_ "V" \ T» 

T qo2 

X 2 > y 36 


~ M 2 

[ X q20 + 

E q20 G 22 + E d20 B 22 

1 

" M 2 

^ I d2o *“ 

E q2o B 22 + E d20 G 22^ 

^2 



M 2 



1 

£ E q20 Y 33 

+ E d20 y 43] 


“ " Pq20 y 34 + E d20 y 44] 

= " t E q2o y 35 + E d2o Y 45 J 
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t 

a 49 ~ " ^ E g20 Y 36 + ^20^46^ 

a 52 = “~T (x d3 - x 3 ) y fil 
x do3 

a 53 = 7* (x d3 " X 3* y 51 
do 3 

a 55 " “/ ^“ 1 + (x d3 ~ X 3^ B 33^ 

do 3 

a 56 = r* ~ (x d3 " X 3* 
do 3 

a 58 = T 7 (x d3 ~ X 3* y 65 
X do3 

a 59 = -~T~ (x d3 - X 3> y 66 
do3 

a 62 = “ ^ X q 3 - X 3> y 51 

qo3 

a 63 "7 (x q3 ~ X 3* y 52 

qo3 

a 65 = “ “77 (X q3 " X 3* G 33 
qo3 

a 66 5=1 77 ^ 1 “ ^ X q3 “ X 3^ B 33^ 

qo3 

a 68 = ~ 7 1 (x q3 *" X 3* y 55 
qo3 
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a 69 = " - 7 — (x q3 - x 3> ^56 

qo3 

a 72 L E q30 Y 51 + E d30 Y 61 3 

a 73 = ~ M^C E q30 Y 52 + E d30 Y 51 3 

a 75 “ “ I q30 + E q30 G 33 + E d30 B 33 3 

1 r 1 1 "i 

a 76 ^dSO ~ E q30 B 33 + E d30 G 33 3 

°3 

a 77 "“Mg 

a 78 = “ t E q30 Y 55 + E d30 Y 65 3 

a 79 = “ C E q30 Y 56 + E d30 Y 66 3 

6.3 CASE STUDY FOR DYNAMIC STABILITY INVESTIGATION 

A case study has been conducted for the dynamic 


stability investigation for the different configurations of 
a sample network in order to study the performance of six- 
phase system in comparision to three-phase systems. 
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6.3.1 A Sample System 

A sample network as shown in Fig, 4.2 has been 
taken for the present investigation at the initial operating 
point obtained from the load flow study. The loads have 
been simulated as constant impedances based upon the initial 
operating conditions. Linearized machine equations have been 
used. The classical model for generator 1 and the two-axis 
model for generator nos. 2 and 3 have been used. The data 
for the 3 synchronous machines are given in Table 6.1. The 
initial operating conditions are obtained with the help of 
the data obtained from the load flow study . provided in 
Section 4.4. 

6.3.2 Different Configurations of the Sample System under 
Investigation 

The following configurations of the sample system 
have been considered for the dynamic stability investigation.! 

1. A compl e tely three-phase system having line nos. 2.3 and 
6 as double -circuit three-phase lines 

2. A mixed system when line no. 2 is converted to a six- 
phase line 

3. A mixed system when both the lines no. 2 and 3 are con- 
verted to six-phase lines. 
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TABLE 6.1 


Three-Machine System Data 


Quantity 

unit 

Generator 1 
(classical) 

Generator 2 
(two-axis) 

Generator 3 
(two-axis) 

H (MW.S/100MVA) Second 

23,64 

6.4 

3.0 

Mj =2H 0) B 

pu 

17824.14 

4825.4863 

2269.4865 

Mj =2H/0 b 

sec/elect 

radian 

.1254 

.03395 

.01597 

X d 

pu 

0.146 

0.8958 

1. 3125 

* 

X d 

PU 

0.0608 

0.1198 

0.1813 

X q 

pu 

0.0969 

0.8645 

1.2578 

X q 

i 

pu 

0.0969 

0.1969 

0.25 

T 

Sec 

0.0 

0.535 

0.6 

qo 

i. 

T 

pu 

0.0 

201.69 

226.19 

qo 

t' 

qo 

sec/elect 

radian 

0.0 

.001419 

.00159 

i 

T, 

do 

sec 

8.96 

6.0 

5.89 

1 

T 

ao 

pu 

3377.8404 

2260.9467 

2220.4777 

i 

T ■ 
ao 

sec/elect 

radian 

.02376 

.01591 

.01562 
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4. A mixed three- phase and six— phase system when all the 
double circuit three-phase lines i.e., line nos. 2,3 
and 6 are converted to single circuit six-phase lines. 

The investigations for configurations (2), (3) 
and (4) have been carried out under two voltage conditions, viz 
A when the phase to ground voltage of a six-phase line is 
V 3 times that of the double-circuit three-phase lines and 
B when the phase to ground voltage of the six-phase line 
and the double circuit three-phase line are the same. 

The dynamic stability investigation has been carried 
out for the above configurations for case a and case b as 
described in Section 4.3, The above investigations have been 
carried out for the same loading of the various configurations 
of the sample network. Further, a mixed three-phase and six- 
phase system have been considered with the increased loading 
i.e., 1.732 times that of the initial loading, when all the 
double circuit three-phase lines 2,3 and 6 have been converted 
to single circuit six-phase lines -under voltage condition A 
for case a (i.e. when phase to ground voltage of a six-phase 
system is V"3 times that of a three-phase system). 

Further, the dynamic stability investigation has 
been conducted on a completely six-phase system with the 
same and with the increased loadings having the same phase to 
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ground voltage as that of the completely three-phase system 
for case a. In this case obviously, three-phase/six-phase 
transformers will not be required, rather six-phase trans- 
formers will be required in place of three-phase transformers* 
The inertia constant and the damping coefficient of six-phase 
synchronous machines have been taken to be equal to those of 
three-phase synchronous machines. The per phase reactances 
and the time constants of six-phase machines have also been 
assumed to be equal to those of three-phase machines. 



The dynamic stability of the different configurat- 
ions has been determined with the help of eigenvalue technique. 
When controlling variables u^ — 0 in equation (6.17), we have 



(6.19) 


The eigenvalues of the system matrix [ A m ] are determined for 
two damping coefficients of the synchronous machines, namely 
Di = 0 or 1.0 pu for i = 1,2,3. 

6.3.3 Results and Discussions 

The eigenvalues of system matrix [A m ]for the 
different configurations of the sample network for two 
damping coefficients have been obtained and are provided in 
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Table 6.2 - 6.12. Table 6.2 shows the eigenvalues of system 
matrix [A^] corresponding to the original completely three- 
phase sample network for two damping coefficients of the machines; 
Table 6.3 to 6.6 shows the eigenvalues for three configurations 
under voltage conditions A and B for two damping coefficients 
respectively for case ajand the Tables 6.7 to 6.10/ those for 
case b. In table 6.11, the eigenvalues of system matrix [ A^ ] 
corresponding to voltage conditions A of case a are tabulated 
for configuration HI with the increased loading for two 
damping coefficients of the machines. The results for the 
completely six-phase system with the same as well as increased 
loadings to that of the original three-phase system are pro- 
vided in Table 6*12. 

The initial value responses have been obtained 
using DVERK subroutine and are shown in Figs. 6.2-6. 7. Figs. 
6.2-6. 3 show the responses (6 l2 and 6 l3 vs «t) for configuration I 
and configuration III under voltage condition A for case a; 
figs. 6. 4-6. 5/ those under voltage condition B for case b ; 
figs. 6.6— 6.7/ those for the completely three-phase system 
and the completely six-phase system for case a under voltage 
condition B. 

From the comparision of Tables 6.2 - 6.6/ it 
is observed that # as the number of conversion from double- 
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TABLE 6,2 

Eigenvalues of System Matrix [a^] for a Completely 
Three-Phase System 


Eigenvalues Dampi n g Coefficients 

Df = D 2 = D 3 - "o D 1 *= D“ =T1 d 3 =1,0 


X l* X 2 

-,926+jl3,0 5 7 3 

1. 004+j 13. 0621 


-,l994+j8,664 

-0.2 34+j 8. 665 

X 5 

-6.28 

- 6.275 

X 6 

-3.908 

-3.91 

X 1 

-.1058 

-0.07501+j ,0487 

X Q 

0.0 


X 9 

-.1727 

-0.1715 
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TABLE 6.3 

Eigenvalues of System Matrix [ A ] under voltage Con- 
ditions A for case a 

(Neglecting danping coefficients of machines) 


Eigenvalues 


System configura 

tions 



configuration I 

configurations 

II configuration III 

X 1' X 2 


1 .0828+] 3.857 

—1.4564 + 3 ' 15. 9 39 

-1.8756+jl8.4579 

X 3' X 4 


-0.24067+j9.162 

-0.4126+jll.031 

-0.52+3*12.27 

h 5 


-5.9948 

-5.2408 

-4.601 

Ns 


-3.76146 

-3.3166 

-2.8197 

X 7 


-0.1188 

-0.1521 

-0.1714 

X 8 


0*0 

0.0 

0.0 

X 9 


-0.1907 

-0.2361 

-0.2936 

System configuration refers to 

following meanings ; 

the mixed system having the 

config.. 

I 

: when line no. 

2 is a six-phase 

line 


config. II : when Ijoth the lines nos. 2 and 3 are six-phase 
lines 

config. Ill : when all the double circuit three-phase lines 
viz line nos. 2,3 and 6 are six-phase lines. 
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TABLE 6.4 

E Igenvalues of System Matrix [ A^J under Voltage 
Conditions B for case a 

(neglecting danping coefficients of machines) 


Eigenvalues 

s 

configuration i 

ystem Configuratic 
configuration II 

ms 

configuration III 

X 1* K 2 

-,883+jl2.8097 

— . 8l48+jl2.6178 

-0.715+j 12.06 

X 3' X 4 

-0.l813+j8.4783 

-0.1572+3 7.974 

-0.1563+j7.892 

X 5 

-6.3166 

—6.405 

-6.139 

X 

-3.8631 

-3.8105 

-3.815 

6 




Kj 

-0.1016 

-0.081 

-0.07825 

X 8 

0.0 

0.0 

0.0 

X 9 

-0.1682 

-0.1594 

-0.1494 
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TABLE. 6.5 

Eigenvalues of System Matrix [ A^] under Voltage Condi- 
tions A for case a 

(uniform darrping coefficients 1.0 pu of machines) 


Eigenvalues 


system configurations 

"configuration I 

configuration II 

configuration III 

A-i #7^2 

-1.16+13.862 

-1.5337+jl5.946 

-1.9522+J18.455 

^3*^4 

-0.2754+j9.152 

-•0,44864^ 11,0318 

— 0. 5569+j 12.272 

X 5 

-5.9896 

-5.235 

-4.5964 

X 6 

-3.764 

-3.3207 

-2.8255 

X 7' X 8 

-0.0812+j.03796 

-0 .06 39 1 1 -0 . 128 6 

-0.1589*— 0.0513 

X 9 

-0.1895 

-0.23593 

-0.2931 
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TABLE 6.6 

Eigenvalues of System Matrix [ A^ ] under Voltage Conditions B 

for case a 

(assuming uniform damping (1.0 pu) of the machines) 


Eigenvalues Syst em conf igurations ; 

CONFYgURATION ' “ Y ~ CONFIGURATION II CONFIGURATION III 


^1 # 

X 2 

— . 9596+j 12 . 814 

-0.8916+j 12.522 

.7922+jl2.064 


X 4 

-0.2l66+j8.4785 

-0.1918+j 7.9744 

-0.l902+j7.891 

^5 


-6.312 

-6.402 

-6.135 

X 6 


-3.866 

-3.8136 

-3,8182 

# 

X 8 

, -.0728+3.0525 

-,0633+j .060172 

-.0623+3.060547 

X 9 


-0.1669 

. -0.1576 

-0.1471 
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TABLE 6.7 


Eigenvalues of the System Matrix [ A^] under Voltage Condi- 
tion A for case b 

(neglecting danping coefficients of machines) 


Eigenvalues 


System configurations 

configuration I 

configuration II 

configuration III 

k 1* x 2 

-,72o2+jl2.l36 

-0.4569+jl0.488 

-0.4198+3*10.083 

^3 '^4 

-0.l734+j8.0088 

-0.2039+3*7.1977 

-0.267+j6.688 


-6.73 

-7.336 

- 7.672 

*6 

—4 .1226 

-4.445 

-4.37 

X 7 

-0.07554 

-0.032 

+0.01 

^8 

0.0 

0.0 

0.0 

X 9 

-0.1457 

-0.09858 

-0.08415 
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TABLE 6.8 


Eigenvalues of System Matrix [a ] under Voltage Condidi- 

tion B for m case b 

(neglecting damping coefficients of machines) 


Eigenvalues 


system configurations 


~ configuration" *T 

configuration II 

configuration III 


*-0.912+12.9 387 

-0.9495^13.1187 

— 0.852+jl2.703 

X 3 #X 4 

-0.1865+j8.654 

-0.1724+J8.32 

-0.l717+j8. 287 

*5 

-6.2578 

-6.134 

-5.852 

X 6 

-3.834 

-3.694 

-3.685 

X 7 

-.105 

-0.098 

-0.0977 

X 8 

0.0 

0 .0 

0.0 


-0.1717 

-0.176 

-0.16837 
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TABLE 6,9 

Eigenvalues of System matrix [ A^ ] Under Voltage Condition A 

for case b 

(assuming uniform damping coefficients 1,0 pu of the 

machines.) 


Eigenvalues 

configuration I 

System configurations 
configuration II configuration ifl 


-0.7976+jl2. 14046 

— 0.5329+j 10.491 

-0.495+jl0.0867 


—0.2069+3*8.009 

-0.2369+37.1979 

-0.2996+3*6.6887 

X 5 

-6.7265 

-7.3346 

-7.671 

X 6 

-4.1254 

-4.4475 

-4.378 

Xr ] ,K 8 

-0.0608+j .06285 

—0.0397+3 .0832 

-0.02115+j.0922 

X 9 

-0.144 

-0.09827 

-0.0809 
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TABLE 6.10 

Eigenvalues of System Matrix [ A^ ] under Voltage Condition B 

for case b 

(assuming uniform damping coefficients (l.Opu) of machines) 


Eigenvalues system configuration 



configuration I 

configuration II 

configuration III 

^1*^2 

— 0.9887+jl2,9433 

— 1 .026+j 13. 12 33 

-0.9293+] 12. 707 

x 3 #x 4 

— 0.2219 +j 8 *654 

-0.2074+j8.324 

-0.2062+j8.287 

X 5 

-6.253 

-6.129 

-5.8577 

*6 

-3.837 

-3.698 

-3.688 

00 

*% 

!> 

<c 

-0.0744+j .05064 

-0.0714+]' .05223 

-0.07l4+j.05l78 


-0.1705 

-0.1745 

-0.1664 
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TABLE 6.11 

Eigenvalues of System Matrix C ]] with the increased System 
Loading under Voltage Condition A . for case a 


Eigenvalues 

A mixed three-phase and six-phase system with 
increased loading 

(Line no* 2,3 and 6 are six-phase lines) 


D 1 - D 2 = D 3 = 1 

D 1 - °2 = D 3 = 0 

K l' X 2 

-1. 2 3018+j 14.0282 

-1.15239+3*14.02258 

\ 3 ,\ 4 

-0,35373+j9.522 

-0.31647+3*9.5208 

>5 

-6.17694 

-6.1825 

>6 

-3.9176 

-3.9112 

x 7 

-0.19403 

-0.1954 

K 8 

-0.1225 

-0.1754 

X g 

-0,08883 

-0.0 
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TABLE 6.12 

Eigenvalues of System Matrix [ A^] for a completely Six- 
Phase System under voltage condition B for case a 


E i gen- with same loading with increased loading 


values 

D ± = i * 0 

d 1 =’D 2 =:D 3 =0 *0 

Di= D 2=D 3 =1.0 

d 1 =d 2 =d 3 =o.o 

X 1' X 2 

-1.9954+j 13.82- 

1.9l73+jl3.8l 

-1.7175+j 12.12 

1.638+5.12 . 115 

X 3'^4 

— 0.60 36 +j 8. 61 

— 0.5669+j8.61 

-0.5496+j7.897 

-0. 51 32+j 7.896 


-3.3828 

-3.3889 

-4.05702 

-4.064 


-2.391 

-2 . 391 

-2.7284 

-2.7263 

X 7 

-0.05026 

—0.0 

-0.05473 

-0.0 

k q 

-0.3765 

-0.37754 

-0.30698 

-0.21794 

K 

-0.21657 

-0.22356 

-0.204318 

-0.3096 



234 

circuit three-phase lines to single-circuit six-phase lines 
under voltage condition A (i.e. , both six-phase- and three- 
phase system having same line to line voltage) are increased, 
the eigenvalues moves towards left. But when the same conver- 
sion is carried out under voltage condition B (i.e. both six- 
phase and three-phase system having same phase to ground 
voltage),, the eigenvalues shift towards right. It thus 
reveals that the conversion to the six-phase system increases 
the dynamic stability as the number of conversion increases 
under voltage configuration A while it decreases under voltage 
configuration B with the loading as that of the original 
three-phase system. Further, when Tables 6.7 - 6.1.0 are 
observed, where the results for case b with the same loading 
as the previous one are tabulated, it shows that the eigen- 
values are shifting right thereby the system being less 
dynamically stable irrespective of the voltage conditions A 
and B with the fixed loading. 

When the results of Table 6,11 <xre compared to 
those of Tables 6.2, and to those of Tables 6.3 and 6.5 under 
configuration III, the. eigenvalues of the system matrix [a^] 
corresponding to the increased loading system under configu- 
ration III for voltage condition A have shifted towards right 
as corrpared to fixed loading and, hence, the power system 
is less dynamically stable with increased loading than with 
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the same loading. If, however, the eigenvalues corresponding 
to configuration III of the system with increased loading is 
compared to those of the completely three-phase system with 
the original loading, the system with the increased loading 
is more dynamically stable than the original three-phase 
system. 


When the results of Table 6.12 are compared with 
those of Tables 6.2, 6.4 and 6.6, it is observed that a 
completely six-phase system with the same loading and also 
with the increased loading is more dynamically stable than a 
completely three-phase system. It i s also observed that, with res 
pect to system configuration under voltage condition B, a 
corrpletely six— phase system with the same loading or with 
increased loading is dynamically more stable than a completely 
"three— phase system.. Further, comparisions of Table 6,12 
with tables 6.3 and 6.5 show that the completely six— phase 
system under voltage condition B is dynamically more stable 
than a mixed system under configuration III for Voltage 
condition A for the same loading, while the former with 
increased loading is dynamically less stable than the' latter. 


From the initial value responses (a..6^ 2 , a 6 13 ~- t) 
shown in Figs. 6.2 - 6.3, it is observed that the dynamic 



I - Line no. 2 is a six-phase line 
III - Line nos. 2.3 and 6 are six-phase lines 



FIG. 6-2 Initial value Responses (hS 17 yt) of mixed-phase 
system under voltage condition A for case a 


and 6 are six- phase lines 



o 




Line no. 2 is a six- phase line 

Line nos 2*3 and 6 are six- phase lines. 



FIG- 6.4 Initial value Responses (A6j 2 ~t ) of mixed- phase 
systems under voltage condition B for case b . 




completely three-phase system, 
completely six -phase system- 



FIG-6-7 InVUal value Responses ( A ) of completely three-phase 
and six- phase systems . 
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stability of the system has considerably improved with the 
increasing conversions of double circuit three-phase lines 
to single circuit six-phase lines under voltage condition A 
for case a. But the dynamic stability of the mixed system 
under voltage condition B for case b as observed from the 
initial responses (figs. 6.4 - 6.5) has not been affected, 
it has rather slightly reduced with increasing conversions. 
•Further, it is observed from Figs. 6.6 and 6.7* that a comp- 
letely six-phase system is better than a completely three- 
phase system from the stability point of view. 

6.4 CONCLUSION 

A dynamic stability investigation has been carried 
out for different system configurations in order to 

observe the relative performances of a completely three-phase 
system* a mixed three-phase and six— phase system, and a comp- 
letely six-phase system. It is observed that a mixed three- 
phase and six-phase system improves the stability with the 
increasing number of conversions from- double circuit three- 
phase lines to single-circuit six-phase lines with the same 
line to line (i.e. adjacent phase) voltage, but the system 
becomes less stable with the increasing conversions taking 
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place at the samephase to ground voltage* Further* it is 
observed that a completely six-phase system with the same or 
increased system loading is dynamically more stable than a 
completely three-phase system* 



CHAPTER 7 


CONCLUSIONS 

The mathematical modelling for the various multi- 
phase elements of power system network has been developed and 
based upon these models* various studies such as load flow* 
transient and dynamic stability studies on a sample network 
for its different configurations have been conducted with a 
view to investigate the performance of a mixed phase as well 
as a completely multi-phase i.e.* completely six-phase system 
in comparision to that of a conventional three-phase system. 

The inductances of a six-phase machine in the ' 
phasor form have been obtained from the fundamentals.. A 
detailed dynamic model of a six-phase machine has been deve- 
loped with the help of an orthogonal transformation. Further* 
a linearized model of the machine has been developed in terns 
of c ur rent variables. The dynamic stability of a six- phase 
machine has been investigated on a sample system* employing 
eigenv al ue technique*. the results of the six-phase synchronous 
machine have been compared to those of the three-phase machine 

Further* the mathematical model of a multi-phase 
(n— phase) synchronous machine in the phasor form has been 
developed and the corresponding model of a six-phase^ as well 
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as twelve phase machine has been derived. For the purpose of 
n -phase transmission, three-phase/n-phase transformers are 
necessary to be compatible with the system; three— phase/n— phase 
transformer models for different winding connections have been 
developed with off-nominal turn-ratios on primary as well as 
on secondary side. 

The three-phase equivalent of a multi-phase (n-phase) 
line connected through two three-phase/n-phase transformers 
has been derived in terms of three-phase impedance matrix 
and the three-phase ABCD parameters for 71 circuit representa- 
tion, and further based upon the three-phase equivalent re-' 
presentation, a single phase equivalent has been obtained. 

The n-phase equivalent of a three-phase line integrated with 
two three-phase/n-phase transformers has also been derived. 

Further, two types of problems have been formulated 
and investigated for the load flow analysis based upon one 
line diagram. The equivalent single-phase representation has 
been derived for the six-phase and twelve-phase transmission 
lines with three -phase/ six-phase or three-phase/ twelve— phase 
transformers at both the ends of multi-phase lines, which 
also helps in obtaining the parameters of single-phase equi- 
valent in pu on the same base as that of three-phase line. 

The load flow investigation has been carried out on a sample 
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network for its different configurations obtained by conver- 
ting double— circuit three-phase lines to multi— phase (six— phase 
or twelve— phase) lines with the phase to ground voltage equal 
to t or Y" 3 times as that of the double— circuit three-phase lines. 

The direct and the quadrature axis transient and 
subtransient inductances of a six— phase machine have been 
obtained from the fundamentals, and the different time cons- 
tants for a six— phase machine have also been calculated. The 
transient stability of a sample network has been investigated 
to study the impact of converting double— circuit three-phase 
lines to single-circuit six-phase lines for the two cases as 
outlined in the load flow study. A completely six-phase 
system with the same voltage has also been investigated for 
transient stability studies. 

The dynamic stability of a multi— machine system 
has been carried out on a sample network by conducting investi- 
gation for the different configurations of the sample network. 
The relative performances of mixed-phase systems, completely 
six-phase system and a conventional three-phase system have 

also been studied. 

Based upon the case studies for mixed-phase and 
six-phase (multi-phase) systems, the performances regarding 
high-phase order have been observed which are discussed briefly 

in the following paragraphs. 
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The six-phase machine is dynamically more stable 
than the three-phase machine when the same power is delivered 
by both the machines * When the six -phase machine delivers 
double the power 9 it shows a little improvement in stability 
as compared to that of a three-phase machine* 

The magnitude of bus voltages of mixed— phase 
systems having the same slack bus voltage increases or 
decreases with the increasing number of conversions of double- 
circuit three-phase lines with equal line to line adjacent 

phase) voltage or equal phase to ground voltage. When the 
same reactive power at the slack bus is considered, the bus 
voltages of the mixed-phase system for the same system loading 
increase or decrease with the same phase to ground voltage or 
with the same line to line (i.e. adjacent phase) voltage, as 
that of the double-circuit three-phase lines, with the increas- 
ing conversions of double-circuit three-phase lines to single- 
circuit six-phase lines. 

The transient stability of mixed-phase system is 
increased with the same as well as with increased loading to 
that of the completely three-phase system, when the double- 
circuit three-phase lines are converted to six-phase lines 
with the same line to line (i.e., adjacent phase) voltages, 

takes place with the same line to 


tout if the conversion 



neutral voltage with the same loading, the stability reduces. 
Also, the stability increases when a three-phase system is 
converted to a completely six-phase system with equal line 
to neutral voltage for the same loading as well as for increased 
loading. It has further been. observed that the impact of the 
conversion on the stability depends upon the - location of the 
converted line from the fault. 

The dynamic stability of multi -machine power 
system shows that a mixed-phase system is dynamically more 
stable than the three-phase system with increasing conversions, 
when the line to line (adjacent phase) voltages of the six-phase 
and the double-circuit three-phase lines are the same; the 
stability decreases with equal line to neutral voltage. When 
mixed-phase system is considered with increased loading, the 
system is more stable than a completely three-phase system. 

For a completely six-phase system with equal and increased 
loading, it shows an improvement in stability as compared to 
that of a three-phase system. 

Based upon the investigation and results, we now 
give the future scope of the work in this area. 



249 


Future Scope of Work 

A three— phase/multi-phase transformer is one of 
the important elements of the mixed-phase systems; so its 
mathematical modelling, an important aspect in the analysis. 

A simple model of the transformer representing it by a series 
leakage impedance along with tappings has been employed in 
the analysis. It is however desirable to develop its model 
including details, such as non-linearities caused by satu- 
ration, phase shifting effects, etc. 

The multi-phase lines have been assumed throughout 
the study as transposed lines, though actually these are . 
un transposed. It is difficult to achieve the complete trans- 
position of the multi -phase lines, it would however be realistic 
to analyse the multi-phase with mutual coupling by taking the 
lines as untransposed. Further, it would be more realistic * 
if the performance of mixed-phase and the completely multi- 
systems (i.e. completely six-phase or twelve-phase systems) 
would be examined having based on the practical data obtained 
from the designs of multi-phase lines for various conductor 
configurations and tower geometries. 

A classical model of the synchronous machine for 
the purpose of transient stability study has been employed 
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where the synchronous machine is represented as a constant 
voltage source behind the direct-axis transient reactance, 
and the loads by constant impedances* It is .desirable to 
employ the machine models incorporating Such details as 
the effect of field flux decay, saturation, effects of exci- 
tation control and governing systems etc. Also it would be 
more realistic if the loads would be represented more accu- 
rately, such as dependence of loads on frequency, etc. 

Although experimental multi-phase lines have been 
constructed and the field testing and the simulation studies 
have been reported, yet there is a need to extend experi- 
mentation in order to gain the operational experience, and 
to be acquainted further with the newer problems as well. 

The dynamic stability investigation of a six-phase 
synchronous machine has been conducted by taking its parameters 
based upon those of three-phase machine. It is, however,/ 
desirable to carry out the investigation based, upon the 
actual parameters from its practical design and tests data* 
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APPENDIX A 


LINE PARAMETERS FOR A THREE— BUS SINGLE -MACHINE SAMPLE 

NETWORK (FIG. 4.1) 


TABLE A.l 


Line Parameters and Various Capabilities of Different 
T YP es of- 138 kV Lines of Sample System Shown in Fig. 4.1 


Type of 
line 

Series Imp. 
Q /phase/ 
mile 

Shunt admi. 
jus/ph/mile 

SIL 

(MW) 

Thermal capac- 
ity (MW) 

Single- 

circuit 

line 

0.1185+ 
jO. 789 

5.387 

49.76 

185.48 

Double- 

circuit 

line 

0.05925+ 
jO. 353 

12.077 

111.39 

370.96 

Six- 

Phase 

Line 

0.1185+ 
jO. 744 

5.7248 

316.96 

642.53 


Lenth of each line = 50 miles. 


TABLE A. 2 

p.u. Line Parameter of Different Types of 138 kV lines on 100MVA 
Base of a Sanple System in Fig. 4.1 


Type of line Series Impedance half shunt admittance 

S t 1 / 2 w 


Single circuit ,0311+j .2071 . 0.0+j.0257 

three-phase line 


Double-circuit ,0155+j .09268 0.0+j.0575 

three-phase line 


Single-circuit .0311+j.l953 0.0+j.0273 

six-phase line 


APPENDIX B 

PARAMETERS OF A NINE-BUS THRE E —MACHINE SAMPLE NETWORK 

(Fig. 4.2) 

TABLE B.l 

p.u. Line Parameters of the Sample Network 


Line 

No. 

Buses 
From To 

m (*T 

Types 
of line 

Series Impedance 

Z 

pq 

Half line changing 
admittance 1/2 y 

1 

4 

5 

1 

.Ol+j ,085 

0.0+j 0.088 

2 

4 

6 

2 

,0l7+j 0.092 

0.0+j 0.079 

3 

5 

7 

2 

.032+30.161 

0.0+j 0.153 

4 

6 

9 

1 

,039+j 0.17 

O.O+jO.179 

5 

7 

8 

1 

,0085+j .072 

0.0-i-j 0.0745 

6 

8 

9 

2 

,0119+j .1008 

0.0+j 0.1045 


* 


Types of line in column 3 have the following meanings 

1 — Single-circuit three-phase line 

2 - Double-circuit three-phase line 


TABLE B. 2 

Pu Transformer Parameters of a Sample System Shown in 


Fig. 4.2 

No. 

connecting 

Bus 

Reactance 

off-nominal 

From 

To 


ratio 

1 

1 

4 

.0576 

1.0 

2 

2 

7 

.0625 

l.o 

3 

3 

9 

.0586 

1.0 


APPENDIX C 


LINEARIZED MODEL FOR MULTIMACHINE SYSTEM WITH CONSTANT 

IMPEDANCE LOADS 


Let and 1^ be the voltages and currents of 
synchronous machines where i = l,2,...,mf m is the number of 
synchronous machines in the network. The complex vectors 
V and I can be expressed as. 



(C. 2) 


where the axis q ± is taken as the phasor reference 
in each case. Let us consider a power system network with 
m machine buses and r load buses. Since the loads are re- 
presented as constant impedances, the network has only m 
nodes connected with the active sources i.e. synchronous 
machines. The reduced network is shown in Fig. C.l For this 

m 
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network, the node voltages and currents are ejqpressed in 
phasor notations and these are V^, V 2 ,...,V m and I I# I 2/ ..,I m 

respectively. At steady— state, these currents and voltages 
can be represented by phasor w.r.t. a common (network) frame 
of reference and threir relations are given by, 

A _ A 

I = YV (C.3) 

where. 



and Y is the short-circuit admittance matrix of the reduced 
network. The network in the transient state can be described 
by equation (C.3). Similar to those describing its steady- 
state behaviour under the following assumptions : 

a) the system angular speed does not depart appreci- 
ably from the rated (i.e. synchronous) speed i.e., 

= 03 

R 

b) the transformer voltages are negligible in compari- 
sion to the speed voltages. 
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Convers ion of Machine Co-ordinates to System Reference 

Let us consider a voltage at any node i. 

It is desired to obtain ^ from v a j Dc j_ by applying Park* s 
transformation. ^ is expressed in phasor notation as. 

"^i ^ rorn ( c «l) using the rotor magnetic axis of machine i as, 
reference axis. It can also be expressed to the system 
reference as using the transformation 

A. _ J 6 i 

V i = V ± e 1 (C.4) 

where 6^ is the angle between this rotor magnetic axis and 
a synchronously rotating reference frame. The rotor angle 
of this machine is given by 

e i = “r* + I + 6 i (C * 5) 

Equation (C.4) can be generalized to include all the nodes, , 
and the generalized transformation is 
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V Q1 +3V D1 

Let * = v Q2+d v D2 

_W 5V Dm_ 

Then from equations (C.l) and (C.4) -(C.7) , we have. 


(C.7) 


A 

V = T V (C.8) 

Thus T is a transformation which transforms the 
d and q quantities of all machines to the synchronously 

rotating system frame of reference. T is an orthogonal 

-1 * 

transformation i,e,, T = T . 

Prom equation no. (C.8), 

V = T* V and similarly I = T I ; I = T* I (C.9) 

Relation between Machine Currents and Voltages 

From equation (C.3), using (C.8) and (C.9), 

T I = Y T V 


i.e. I = M V 

M = (T"* 1 Y T) 


(C.10) 


where 


(C.ll) 
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matrix Y 

of the network is 

of the form 

~~ J©n 
Y ll e 

-*®12 

Y 12 e 

Y, e JSlm 
lm 

jG 2l 

Y 2l e 

-j©22 

Y 22 e 

Y e jG2m 
Y 2m S 

* 

* * • * 

• • • 9 

• 

« 

j 9 - 
Y e ^ 

j9 m2 

Y m2 e ••• 

Y mm eJSrrm 


(C.12) 


and 


M = 


■11 


'21 


, j9 ll 

j ( e 12 “ 6 i2^ 

j (©. 

Y 12 e 

* * Y im e 

j (©21* 6 21 > 

J '©22 

j (©. 


lm 


2m °2iirr 


•22 


'2m 


g J (8 m2-« m 2) 

ml m2 mm 


(C.13) 


Thus equation (C.lo) gives the desired relation between the 
terminal voltages and currents of the machines. 

Linearized Model for the Network 

Linearizing equation (C. 10), we have. 


A i s H q aV + a M V o 


(C.14) 
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where M q is evaluated at the initial angles i = 1,2, . .m , 

a nd V o is the initial value of the vector V. 

Let 6^ = 6^ o + A then the matrix M in equation (C.13) 
becomes 


M= 


Y ll e 


j© 


11 


j ( e 1 2’" 6 i2o“ A<5 12^ J {G lm" 6 lmo“ A6 ln 

*12 e ** *lm 


Y mir e 


J ( G ml' 6 ml0" Ae rd \ r J ( 0 m2 ‘■ 6 m2o’" A5 m2 5 J G 


Y m2 e 


. . Y e 

mm 


mm 


The general term m. . of the matrix M is Y . . e 

XJ Xj 


(c.15) 

j ( © . . —6 . . ~a<5 . . } 
ij ijo ij 


— — j ^ e ii“^iio^ ^ 

Thus m. . of the matrix M is Y,.e J J e 


ij 


'ij 


±j (C.16) 


Substituting cos a 6. • = 1 and sin a 6. . = A 6 . . , (C.16) becomes, 

-kJ 


j( 9 ii~ 6 iio ) 

m, , = Y, , e 1J 1JO (1-j AS,,) 


ij U 


i3 


(C.17) 


Therefore, 


— - . v ^®ij”*^ijo^ 

Amy = -J YyS 


(C.l 8) 
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Hence, a MV 

o 


Thus, the linearized equation (C. 14 ) 


a 1 

i 


Y ll e 

♦ * ♦ 

^2 

ss 

j(9 21 -S 210 ) 

Y 2l e 

• • « 

* * 

AT 

m 


• *■ • •«. 

^ ie 3 ( ®nl- 6 nlo 5 

• mm 

• • m 


S Y ,. e (e lk- S lW) > 
k=l ik 


*2V e 

k=l /Jc 


m j ( q 4-6 ) 

2 Y e mk mk ° ; 
k=l mk 


V k0 A6 lk 


v k0 A 6 2k 


V A 6 

kO °mk 


(C.19) 


becomes. 


Y e 3 ' ^ e im _6 lm0^ 
lm 


J* 

L 

Y e J(e 2 m - S 2 m0 > 

2m 


av 2 

j© 

Ye 1 ™ 

mm 


AV 

m 

, 

i 

— — j 






— — : 


? 7 Y e j(9 lk-W 

kO lk 


A 6 lk 

3 

Z V Y e J t8 2k _a 2kO ) 

k=l k ° Y 2k e 


A6 2k 


m j (© , -6 , J) , 

2 V, Y v e mk mk0 
i=l k0 mk 


• m m m 

A 6 mk 


(C.20) 
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The set of equations (C.2o) gives the conplete 
description of the system. 

Now the above equation (c.2o) has been considered 
for the sample system with three synchronous machines. 

Separating real and imaginary parts of equation 
(C.2o) we have* 


AI 

qi 


y ll 

y 12 

Y 13 

y 14 

y 15 

y 16 


AE q2. 

Al dl 


y 21 

y 22 

Y 2 3 

y 24 

y 25 

y 26 


AE d2 

AI q2 


G 22 

” B 22 

y 33 

y 34 

y 35 

y 36 


ae' . 
q3 

• AI d2 


B 22 

G 22 

Y 4 3 

y 44 

y 45 

y 46 


AE d3 

AI q3 


y 51 

y 52 

G 33 

~ B 33 

y 55 

y 56 : 


A6 12 

AI d3 


y 61 

y 62 

B 33 

°33 

y 65 

y 66 , 

; 

r 

A6 13 


y U = Y 22 = Y 12 cos < e i2- 6 l2 0 ); y 2l = “ y i2 = Y 12 sini (e i2“ 6 12o> 

Y l3 = Y 24 = Y 13 COS (e i3- 6 130>- Y 23 = ~ y i4 * Y 13 Sin 

Y 15 = Y 12 [ E d20 COS (e i2- 6 l2 0 )+E q20 sin ^12’’ 6 12o ) 3 

y 16 ~ Y 13^ E d30 COS (9 13“ 6 130 ) + E q30 Sin ( ®13” 6 130^ 

y 25 * Y 12 C E d20 Sin (e i2“ 6 120^ *“ E q20 GOS ^ e l2*“ a i2o^ 3 



269 


y 26 

= y 13 

r e * 

L d30 

sin 

(6 13“ 6 130 5 

- E ' 

q30 

cos 

(©13-6130) ] 

y 33 

= y 44 

= ^23 

cos 

^ e 23” 6 230 ) 

* y 43 = 

“ y 34 

= y 23 sin (©. 

y 35 

= 

Y^2 sin (© 

12 +5 120 ) “ 

i 

E Y 

d30 23 

cos 

(©2 3 -6 2 30 ) 



- 1 

- E 

q30 Y 23 

sin (© 23 

- 5 

230* 

y 36 

" Y 2 3 l 

r e 

L d30 

cos 

(e 23“ 6 2 30* 

1 

+ E 

g30 

sin 

(9 23 - 6 230>J 

y 45 


^ 0 12 

+6 120 ) - E 

d30 Y 23 Sin '®23- 6 230 ) 





+ Vo Y 23 cos 

(© 2 

3“ 6 230* 

y 46 

= Y 23 

^ E d30 

sin 

(9 23“* 6 230 

) - E* 

q30 

cos 

^ 0 23 _6 23O* ^ 


y 51 = y 62 " Y 23 cos <S 23 + 6 23o ,; y 61 * ~ y 52 = Y 23 sln ( ®23 +6 23 0 > 
y 55 = Y 23 t E d2 0 oos (S 23 +5 230 ,+E q2o sin (e 23 + S 230 > ] 
y 56 “ - E 1 Y 13 sin <®1 3 + 6 13o> “ Y 23 E d20 OOS C ®23 + 6 230 > 

_Y 23 E q20 Sln ^®23 + S 230^ 

y 65 = Y 23t E d20 sln ( ®23 + 'W ~ E g2 0 oos ( ®23 + 6 230 ) ] 
y 66 = E 1 Y 13 c ° s <®13 + W " E d2o Y 23 sin <®23 + W 

+ E q20 Y 23 COS ( ®23 + S 230 > 
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